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THERMO-MECHANICAL  AND  THERMAL 


BEHAVIOR  OF  HIGH-TEMPERATURE  STRUCTURAL  MATERIALS 


PREFACE 


Technical  ceramics  because  of  their  chemical  inertness,  high  melting 
point,  good  wear  resistance,  excellent  mechanical  stability  at  high  tempera¬ 
ture  and  other  unique  properties,  represent  a  class  of  materials  eminently 
suited  for  many  critical  engineering  applications.  Unfortunately,  because 
of  their  brittleness  and  unfavorable  combination  of  pertinent  material 
properties,  technical  ceramics  generally  are  highly  susceptible  to  catas¬ 
trophic  failure  in  non-uniform  thermal  environments,  which  give  rise  to 
thermal  stresses  of  high  magnitude. 

Thermal  stress  failure  analysis  of  structural  materials  represents 
a  multi-disciplinary  problem  which  involves  the  principles  of  heat  trans¬ 
fer,  mechanics  and  materials  engineering.  Over  the  last  few  decades  much 
general  understanding  of  the  nature  of  thermal  stress  failure  of  brittle 
materials  has  been  generated.  However,  due  to  the  multi-disciplinary 
nature  of  the  problem,  the  ability  to  predict  thermal  stress  failure  quan¬ 
titatively  for  design  or  other  purposes  has  lagged  behind  the  progress  made 
in  other  engineering  fields.  The  objective  of  the  present  program  is  to 
improve  the  qualitative  and  quantitative  understanding  of  the  nature  of 
thermal  stress  failure  of  brittle  structural  materials,  including  the  ex¬ 
perimental  as  well  as  theoretocal  variables.  In  order  to  achieve  this 
objective,  the  participating  investigators  and  scope  of  the  program  are 
organized  such  that  full  advantage  is  taken  of  the  com!  ined  inputs  from  a 
number  of  engineering  disciplines.  In  a  similar  spirit  a  number  of  studies 
were  conducted  in  cooperation  with  investigators  at  other  institutions. 

The  effort  of  this  program  consists  of  four  main  themes,  including: 
experimental  thermal  shock  testing,  measurement  of  thermophysical  pro¬ 
perties  relevant  to  thermal  stress  failure,  the  analysis  of  thermal  stress 
failure  and  the  dissemination  of  information  on  thermal  stresses  in 
the  form  of  review  articles,  conferences,  etc. 

Studies  completed  within  the  period  covered  by  this  report  are  pre¬ 
sented  as  individual  chapters  in  the  main  body  of  this  report.  The  title 
of  these  chapters  within  their  main  theme  together  with  a  brief  comment 
are  as  follows: 


i 


A.  Experimental  Thermal  Shock  Testing 


Chapter  I:  "Effect  of  Bath  and  Specimen  Temperature  on  the  Thermal 
Stress  Resistance  of  Brittle  Ceramics  Subjected  to 
Thermal  Quenching" 

This  study  shows  that  the  heat  transfer  coefficient  of  fluid  quenching 
media  can  be  highly  variable.  Such  variability  must  be  taken  into  account 
in  the  quantitative  interpretation  of  experimental  data  for  the  relative 
thermal  stress  resistance  of  different  materials. 

Chapter  II:  "Improvement  of  Thermal  Shock  Resistance  of  Brittle 

Structural  Ceramics  by  a  Dispersed  Phase  of  Zirconia" 

This  study  shows  that  the  thermal  shock  behavior  of  structural  ceramics 
such  as  alumina  and  silicon  nitride  can  be  improved  significantly  as  the 
result  of  toughening  by  the  zirconia  dispersed  phase. 

B.  Measurement  of  Thermophysical  Properties  Relevant  to  Thermal 
Stress  Failure 

Chapter  III:  "Effect  of  Phase  Composition  and  Microstructure  on  the 
Thermal  Diffusivity  of  Silicon  Nitride" 

The  results  of  this  study  indicate  that  the  heat  transfer  properties 
of  hot-pressed  and  reaction-sintered  silicon  nitride  are  strongly  affected 
by  the  a/a+B  ratio,  grain  boundary  phase,  structure  and  size  of  pore  as 
well  as  grain  morphology. 

Chapter  IV:  "Orientation  Effects  on  the  Thermal  Diffusivity  of 
Hot-Pressed  Silicon  Nitride" 

The  thermal  diffusivity  of  silicon  nitride  was  found  to  be  greater  in 
perpendicular  than  parallel  to  the  hot-pressing  direction.  This  effect  re¬ 
sults  from  the  preferred  crystallogrophic  orientation  of  the  B-Si^N^  grains. 

Chapter  V:  "Observations  on  the  Sub-Critical  Growth  and  Healing 
of  Microcracks  in  Brittle  Ceramics" 

The  stability  of  microcracks  in  brittle  structural  materials  was 
monitored  by  measurements  of  the  thermal  diffusivity.  The  experimental 
observations  of  crack-growth  and  healing  provided  evidence  that  the  micro¬ 
crack  configurations  depend  strongly  on  temperature  levels  as  well  as  time. 

C.  Analysis  of  Thermal  Stresses 

Chapter  VI:  "Thermal  Stresses  in  Partially  Absorbing  Flat  Plate 

Symmetrically  Heated  by  Thermal  Radiation  and  Cooled 
by  Convection" 


Chapter  VII: 


"Analysis  of  Thermal  Stress  Resistance  of  Partially 
Absorbing  Ceramic  Plate  Subjected  to  Asymmetric 
Radiation,  II,  Convective  Cooling  at  Front  Surface" 


The  above  two  studies  are  in  the  general  area  of  thermal  stresses 
generated  by  a  spatially  non-uniform  internal  heat  generation,  relevant 
to  IR-windows,  solar  energy  collectors  and  nuclear  fusion  reactors.  The 
results  of  these  studies  establish  the  role  of  the  material  properties  and 
heat  transfer  variables  which  affect  thermal  stress  resistance  under  the 
particular  heating  and  cooling  conditions  considered. 

Chapter  VIII:  "Effect  of  Nature  of  Concavity  of  Temperature  Distribution 
on  Position  and  Sign  of  Maximum  Thermal  Stress" 

In  a  previously  reported  study  of  the  thermal  stresses  in  a  semi¬ 
absorbing  ceramic  plate,  heated  at  the  front  by  radiation  and  cooled  at 
the  back  by  convection,  it  was  found  that  the  maximum  tensile  stresses 
occurred  in  the  front  surface,  which  is  at  the  highest  temperature.  This 
finding  contrasts  with  the  usual  result  found  for  many  other  instances  of 
thermal  stresses,  in  which  generally,  the  regions  at  highest  temperature  are 
in  a  state  of  compression.  As  demonstrated  in  this  chapter,  the  occurrence 
of  the  maximum  tensile  stresses  in  the  hottest  region  of  the  plate  is 
due  to  the  concave  downward  nature  of  the  curvature  of  the  temperature  dis¬ 
tribution  within  the  plate  due  to  the  internal  heat  generation.  In  the 
case  of  heat  transfer  at  the  surface  as  is  generally  the  case  in  the  vast 
majority  of  instances  of  transient  heat  transfer,  the  temperature  distribution 
is  concave  upward,  which  results  in  a  compressive  stress  state  at  the  position 
of  highest  temperature. 

D.  Dissemination  of  Information  on  Thermal  Stresses 

Chapter  IX:  "Effect  of  Crack  Interaction  on  the  Fracture  Initiation 
and  Crack  Propagation  in  Brittle  Ceramics  Subjected  to 
Severe  Thermal  Shock" 

It  is  shown  that  the  interaction  between  cracks  generally  is  beneficial 
to  thermal  stress  resistance,  with  the  exception  of  co-planar  cracks  which 
can  coalesce  to  result  in  complete  structural  disintegration. 

Chapter  X:  "Stresses  Due  to  Thermal  Trapping  in  Semi-Absorbing 
Materials  Subjected  to  Intense  Radiation" 

This  study  reviews  the  principal  theoretical  results  of  the  thermal 
stresses  in  semi-absorbing  ceramics  subjected  to  thermal  radiation. 

Chapter  XI:  "Effect  of  Spatial  Variation  of  Thermal  Conductivity  on 

Magnitude  of  Tensile  Thermal  Stresses  in  Brittle  Materials 
Subjected  to  Convective  Heating" 

This  article  reviews  the  results  reported  previously,  which  shows  that 
a  spatial  variation  of  thermal  conductivity  can  significantly  lower  the 
magnitude  of  the  tensile  thermal  stresses  in  convective  heating. 

Within  the  general  theme  of  the  dissemination  of  information  on  thermal 
stresses, a  conference  was  organized  on  the  topic,  "Thermal  Stresses  in  Material 
and  Structures  in  Severe  Thermal  Environments,"  held  at  Virginia  Polytechnic 
Institute  and  State  University,  March  19-21,  1980.  The  proceedings  of  this 
conference  which  contain  a  total  of  forty-three  papers  on  thermal  stresses  in 
materials  and  structures  are  available  from  Plenum  Press. 


For  the  convenience  of  the  reader.  Table  1  lists  the  status  of 
publications  presented  in  earlier  reports.  For  most  of  these  studies, 
reprints  are  available  and  will  be  sent  on  request. 
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CHAPTER  I 


EFFECT  OF  BATH  AND  SPECIMEN  TEMPERATURE 
ON  THE  THERMAL  STRESS  RESISTANCE  OF  BRITTLE 
CERAMICS  SUBJECTED  TO  THERMAL  QUENCHING 


by 


J.  P.  Singh,  Y.  Tree  and  D.  P.  H.  Hasselman 


Department  of  Materials  Engineering 
Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  VA  24061  USA 


' 


ABSTRACT 


The  effect  of  specimen  and  bath  temperature  on  the  failure  of 
brittle  ceramics  in  a  thermal  quench  experiment  was  studied  by  quench¬ 
ing  glass  and  alumina  rods  in  water  and  silicone  oil  baths  at  different 
temperatures.  The  results  were  discussed  in  terms  of  the  variation  of 
heat  transfer  coefficient  of  the  quenching  media  and  the  change  in 
material  properties  as  a  function  of  temperature.  It  was  found  that 
the  usual  assumption  of  constant  heat  transfer  coefficient  and  material 
properties  may  lead  to  considerable  errors  in  the  quantitative  inter¬ 
pretation  of  the  results  of  thermal  quench  experiments.  Effective  values 
for  the  film  coefficient  of  heat  transfer  for  water  and  oil  baths  were 
estimated  as  a  function  of  film  temperature  from  thermal  quench  data. 
Recommendations  were  made  for  the  selection  of  quenching  media  and  for 
the  procedure  to  be  followed  in  reporting  the  results. 
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1  .  INTRODUCTION 

Because  of  their  favorable  properties  at  high  temperature  as  well  as 
corrosive  and  erosive  environments,  ceramic  materials  are  finding  ever 
wider  applications  in  industrial  and  many  other  fields.  Because  of  their 
inherent  brittleness  and  combination  of  other  pertinent  properties, 
ceramic  materials,  however,  exhibit  poor  resistance  to  catastrophic  failure 
under  thermal  conditions  which  generate  thermal  stresses.  Therefore,  for 
the  purpose  of  reliable  engineering  design,  it  is  important  that  the 
variables  which  control  thermal  stress  failure  of  brittle  ceramics  are 
well  understood.  For  a  given  thermal  environment,  the  selection  of 
materials  with  optimum  thermal  stress  resistance  can  be  based  on  "thermal 
stress  resistance  parameters"  obtained  from  theory1  ^  or  on  comparative 
testing. 

One  such  method  of  testing  thermal  stress  resistance  of  ceramic 

materials  consists  of  quenching  appropriate  specimens  from  high  temperature 

7-9 

into  a  fluid  medium  at  lower  temperature.  Typical  specimens  are  rods 

with  circular  or  square  cross-section.  Quenching  media  can  include  water, 
oils  or  fluidized  beds.  Although,  this  test  method  enjoys  considerable 
popularity  because  of  its  relative  simplicity,  the  experimental  results 

g _ x  o 

obtained  frequently  do  not  agree  well  with  predictions  based  on  theory. 

In  part,  such  discrepancies  possibly  can  be  attributed  to  differences  in 
the  susceptibility  of  ceramic  materials  to  subcritical  crack  growth 
which  was  shown  to  be  significant  in  the  water  quencihng  of  a  glass.11 

Other  possible  variables  to  which  a  discrepancy  between  theory  and 
experiment  can  be  attributed,  can  be  ascertained  by  the  theoretical 
expression  for  the  critical  temperature  difference  (AT  )  required  to 
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induce  thermal  stress  failure  in  a  ceramic  specimen  subjected  to  a  thermal 
quench: 

AT,  =  A  St(l-v)/aE  +  B  S  (l-v)K/aEah  (1) 

where  A  and  B  are  constants  which  depend  on  the  geometry  of  the  specimen, 

a  is  the  characteristic  dimension  (i.e.,  radius  of  a  circular  cylinder), 

S  is  the  tensile  strength,  E  is  Young's  modulus,  v  is  Poisson's  ratio,  a 

is  the  coefficient  of  thermal  expansion,  K  is  the  thermal  conductivity  and 

h  is  the  heat  transfer  coefficient. 

Eq .  1  indicates  that  ATc  depends  on  a  number  of  material  properties, 

the  dimensions  of  the  specimen  and  the  heat  transfer  coefficient,  h.  The 

properties  such  as  the  coefficient  of  thermal  expansion.  Young's  modulus, 

Poisson's  ratio  and  thermal  conductivity  as  well  as  their  temperature 

dependence  are  readily  measured.  An  assessment  of  the  tensile  strength  is 

more  difficult  due  to  the  statistical  nature  of  brittle  fracture.  This 

latter  phenomenon  requires  the  need  of  establishing  fracture  probabili- 

lities  based  on  stress  distributions,  specimen  dimensions,  etc.,  by  such 
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theories  advanced  by  Weibull  or  others.  The  need  for  such  an  approach 

13 

was  established  in  a  previous  study  of  the  thermal  fatigue  of  a  glass. 

The  constants  A  and  B  are  readily  determined  by  analytical  or  numerical 
methods. ^  ^  Specimen  dimensions  can  be  obtained  to  high  precision.  This 
leaves  the  heat  transfer  coefficient  as  the  last  remaining  significant 
variable,  which  could  be  responsible  for  the  discrepancy  between  theory 
and  experiment.  Generally,  in  the  comparison  and  interpretation  of 
quenching  data,  the  heat  transfer  characteristics  of  the  fluid  medium 
are  assumed  to  be  invariant. 
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To  establish  the  validity  of  this  latter  assumption,  it  should  be 
noted  that  the  transfer  of  heat  between  a  solid  and  fluid  is  governed  by 
a  number  of  variables.  Firstly,  the  heat  transfer  coefficient  depends  on 
the  fluid  properties  such  as  the  coefficient  of  volume  thermal  expan¬ 
sion,  the  density,  the  specific  heat,  the  thermal  conductivity  as  well 
as  the.  viscosity.  For  purposes  of  the  objectives  of  this  study,  it  is 
important  to  note  that  these  fluid  properties  especially  the  viscosity 
can  show  a  strong  temperature  dependence.  Secondly,  the  rate  of  heat 
transfer  depends  on  the  geometry  and  dimensions  of  the  solid  and  the 

difference  in  temperature  between  the  solid  surface  and  fluid. ^  Thirdly, 

* 

for  a  specimen  falling  through  the  fluid  at  its  terminal  velocity,  the 

rate  of  heat  transfer  also  is  a  function  of  the  relative  difference  in 
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density  of  the  fluid  and  solid.  For  a  water  bath  ,  the  heat  transfer 
is  also  affected  very  strongly  by  the  effect  of  nucleate  boiling 
which  enhances  the  heat  transfer  coefficient  and  the  effect  of  film 
formation  at  very  high  temperature  levels  which  tends  to  suppress  the 
transfer  of  heat. 

Calculations  of  heat  transfer  coefficients  generally  are  based  on 
the  properties  of  the  fluid  corresponding  to  the  "film"  temperature 
taken  as  the  average  of  the  temperature  of  the  fluid  and  the  surface  of 
the  solid.  It  should  be  noted  then  that  ceramic  specimens  with  differences 
in  values  of  the  critical  quenching  temperature  (AT^)  will  show  differences 

*  Terminal  velocity  is  the  constant  velocity  attained  by  the 
specimen  during  its  fall  through  the  fluid  bath  when  the  downward 
force  of  its  weight  is  balanced  by  the  upward  buoyancy  and  shear 
forces . 
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in  the  film  temperature  even  for  identical  bath  temperatures.  In  view 
of  the  temperature  dependence  of  the  fluid  properties,  especially  the 
viscosity,  differences  in  values  of  AT^  will  lead  to  differences  in  the 
values  of  the  heat  transfer  coefficient.  At  least  in  part,  then,  the 
aforementioned  discrepancies  between  predicted  and  measured  thermal 
stress  resistance  in  quenching  experiments,  is  because  the  heat  transfer 
coefficient  is  a  function  of  temperature  rather  than  a  constant  as  usually 
assumed. 

The  validity  of  this  hypothesis  can  be  verified  by  measuring  AT^ 
for  a  range  of  film  temperatures  which  can  be  achieved  by  varying  the 
bath  temperature.  For  a  given  material  and  bath  temperature,  the  film 
temperature  can  be  varied  by  changing  the  size  of  the  specimen.  The 
purpose  of  the  present  paper  is  to  report  the  results  of  such  a  study. 

2 .  EXPERIMENTAL 

2.1.  Materials  and  Quenching  Media 

Specimens  consisted  of  circular  rods  6.25  cm  long  and  diam.  of 
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0.53  cm  of  a  soda-lime-silica  glass**  also  used  in  previous  studies 
and  a  polycrystalline  aluminum  oxide***  with  a  range  of  diameters.  The 
values  for  the  material  properties  of  the  glass  are  presented  in  Table  1. 
Similarly  Table  2  lists  the  flexural  strength  of  the  aluminum  oxide  as 
a  function  of  diameter  size  and  other  relevant  properties. 

The  quenching  media  consisted  of  water  and  two  silicone  oils^  with 
nominal  viscosities  of  5  and  100 centistokes  at  25°C.  The  alumina  specimens 
were  tested  in  both  the  water  and  silicone  oils.  The  glass  rods  were  tested 
in  the  silicone  oils  only,  in  view  of  the  considerable  stress  corrosion 

**  R-6,  Owens  Corning  Fiberglass  Corp.,  Toledo,  Ohio 
***  998  Alumina,  McDanel  Refractory  Porcelain  Co.,  Beaver  Falls,  PA. 
t  Type  200,  Dow-Corning  Corp.,  Midland,  Michigan. 
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TABLE  1.  PROPERTIES  OF  Na20-Ca0-Si02  GLASS  11 

Flexural  Strength  at  Liquid  Nitrogen  Temperature  (MPa)  194 
Young's  Modulus  of  Elasticity  (GPa)  69 

Coefficient  of  Thermal  Expansion  (°C.  1)  9.3x10 

Poisson's  Ratio  0.25 

Thermal  Conductivity  (Cal. cm.  loC  1 . S .  1)  2.5x10 
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TABLE  2.  PROPERTIES  OF  POLYCRYSTALLINE  ALUMINA^ 


Flexural 

Strength  (MPa) 

Strength  Relative  to  Diameter 

Diam  (cm) 

of  0.318  cm 

0.318 

327  (±25) 

1.0 

0.475 

311  (±64) 

0.95 

0.635 

271  (±51) 

0.83 

0.953 

251  (±20) 

0.77 

Young's  modulus  of  elasticity  (GPa)  393 

Coefficient  of  Thermal  expansion  (°C  7.4x10 

Poisson's  Ratio  0.265 

Thermal  Conductivity  (Cal, cm.  l0C.  ^S. 

200°C  0.0437 

300°C  0.Q312 

350° C  0.0291 


ft  The  data  on  flexural  strength  and  the  thermal  conductivity  were  obtained  in 
our  laboratory  and  the  other  material  properties  were  supplied  by 
the  manufacturer. 
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effects  encountered  in  water.  * 

2.2  Quenching  Procedure 

The  quenching  media  were  contained  in  a  stainless  steel  beaker 
with  diameter  of  18  cm  by  20  cm  deep.  The  steel  beaker  was  held  in 
a  thermostatically  controlled  glass  tank  30  cm  in  diam.  by  30  cm  deep. 
Prior  to  quenching,  the  ends  of  the  glass  rods  were  wrapped  with  glass 
insulating  tape  to  prevent  fracture  due  to  unknown  stresses  at  the  rod 
ends.  Because  of  the  very  high  density  of  cracks  in  the  alumina,  this 
procedure  was  not  required  for  these  specimens. 

Sets  of  five  specimens  were  heated  slowly  to  a  pre-selected  tem¬ 
perature  in  a  small  electrically  heated  laboratory  furnace  and  held  for 
10-15  minutes  to  assure  thermal  equilibrium,  after  which  they  were  dropped 
into  the  quenching  bath.  This  procedure  was  repeated  for  a  range  of 
quenching  temperature  difference  (AT) ,  defined  as  the  difference  in 
initial  temperature  of  the  specimen  and  the  bath  temperature.  For  the 
glass,  the  critical  quenching  temperature  difference  (ATc)  was  defined 
as  the  value  of  At  at  which  three  or  more  specimens  had  visible  cracks. 

For  the  alumina  specimens  in  which  the  cracks  were  not  easily  detectable, 
the  critical  temperature  difference  was  obtained  by  measurements  of  the 
strength  in  4-point  bending  with  loading  and  support  spans  of  1.6  and 
5  cm,  respectively,  and  a  crosshead  speed  of  0.05  cm.  min  Prior  to 
the  strength  test,  all  specimens  quenched  in  water  were  dried  at  ,vll0oC 
for  2-3  hours.  The  critical  temperature  difference  was  reported  in  terms 
of  a  temperature  range.  The  lower  limit  of  this  range  corresponded  to 
the  value  of  AT  at  which  one  specimen  in  five  exhibited  a  significant 
strength  loss.  The  upper  limit  corresponded  to  the  value  of  AT  at  which 
four  or  more  specimens  exhibited  a  strength  loss. 
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3.  RESULTS  AND  DISCUSSION 

Figs,  la  and  lb  show  the  value  of  AT^  as  a  function  of  bath  tem¬ 
perature  for  the  glass  specimens  quenched  into  the  silicone  oils  with 
viscosity  of  100  and  5  cst,  respectively  (at  25°C) . 

The  decrease  in  ATc  with  increasing  bath  temperature  for  the  oil 
with  viscosity  of  100  cst  above  30°C  is  thought  to  be  due  primarily  to 
the  decrease  in  viscosity  with  increasing  temperature.  Such  a  decrease 
in  viscosity  increases  the  heat  transfer  coefficient  which  in  turn  de¬ 
creases  AT^.  On  comparison  of  figs,  la  and  lb,  this  conclusion  is  sup¬ 
ported  by  the  fact  that  for  a  given  bath  temperature  the  oil  with  the 
higher  viscosity  also  yields  the  higher  value  of  ATc_  At  least  for 
the  100  cst  oil,  the  variation  of  AT£  with  temperature  cannot  be  attri¬ 
buted  to  the  temperature  dependence  of  the  relevant  material  properties 
of  the  glass,  since  the  initial  specimen  temperature  (i.e.,  bath  tempera¬ 
ture  +  AT^)  is  constant  over  the  total  range  of  bath  temperature. 

Of  interest  in  fig.  lb  is  the  decrease  in  AT^  with  decreasing 
bath  temperature  below  30°C.  This  behavior  is  directly  opposite  to  the 
effect  noted  for  temperatures  in  excess  of  30°C.  It  is  thought  that  this 
apparent  anomaly  can  be  attributed  to  at  least  two  effects.  Firstly, 
a  calculation  of  the  heat  transfer  coefficient  from  the  fluid  properties 
near  room  temperature  indicates  that  the  temperature  dependence  of  the 
heat  transfer  coefficient  decreases  with  decreasing  fluid  temperature.  For 
this  reason,  the  effect  of  bath  temperature  on  ATc  is  expected  to  be  less 
pronounced  at  the  lower  than  at  the  higher  bath  temperatures.  Secondly, 
it  should  be  noted  that  with  decreasing  bath  temperature  the  initial 
specimen  temperature  also  decreases.  This  implies  that  the  temperature 


CRITICAL  QUENCHING  TEMP  OIFF,  &TC  (°C)  CRITICAL  QUENCH 


9 


350 


340 


330 


320 


310 


300 


_l _ J _ I _ I _ 1 _ L- 


0  10  20  30  40  50  60  70  60 


BATH  TEMPERATURE  ( °C ) 


ig.  1.  Effect  of  bath  temperature  on  ATC  for 
glass  rods  quenched  in  silicone  oil: 
(a)  100  centistokes  (b)  5  centistokes 
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dependence  of  the  properties  of  the  glass  should  be  taken  into  account. 

L8 

In  this  respect,  the  rather  high  temperature  dependence  of  Young's  modulus 
and  strength  of  a  soda-lime-silica  glass  over  this  range  of  initial  specimen 
temperature,  tends  to  reduce  A with  decreasing  bath  temperature.  For 
a  range  of  bath  temperature  of  10  to  30°C, the  strength  increases  due  to 
increasing  initial  specimen  temperature  prior  to  quenching  beyond  which 
the  strength  remains  constant  due  to  approximately  constant  initial  specimen 
temperature.  The  anomaly  in  ATc  for  the  low  bath  temperature  shown  in 
fig.  lb,  results  from  the  combined  effects  of  the  temperature  dependence 
of  the  heat  transfer  coefficient,  Young's  modulus  and  failure  stress,  with 
the  temperature  dependence  of  these  latter  two  properties,  contributing  the 
major  effect.  In  view  of  the  relatively  low  temperature  dependence  of  the 
coefficient  of  thermal  expansion  and  the  thermal  conductivity,  these 
two  properties  are  thought  to  have  little  influence  on  the  anomaly. 

Support  for  this  latter  conclusion  is  given  by  the  experimental 

data  shown  in  fig.  2  for  AT  for  the  alumina  rods  with  diam.  of  0.63  cm 

c 

quenched  in  the  same  silicon  oil  (5  cts)  as  the  glass  data  shown  in 
fig.  lb.  For  these  alumina  specimens  the  anomaly  at  the  low  value  of 
bath  temperature  is  absent.  The  value  of  AT^  decreases  monotonically 
with  increasing  bath  temperature  probably  due  primarily  to  the  decrease  in 
viscosity  and  corresponding  increase  in  heat  transfer  coefficient  with 
increasing  bath  temperature.  At  the  initial  specimen  temperature  of  450°C 
the  aluminum  oxide  is  not  expected  to  exhibit  a  significant  change  in 
fracture  stress  during  the  15  min.  anneal.  This  suggests  that  indeed  the 
anomaly  shown  in  fig.  lb  is  due  to  the  temperature  dependence  of  the 
strength  and  Young's  modulus  of  the  glass  over  the  range  of  initial 


CRITICAL  TEMPERATURE  DIFFERENCE,  ATc(°C) 


Effect  of  bath  temperature  on  ATC  for  alumi¬ 
na  rods  (d  =  0.63  cm.)  quenched  in  silicone 
oil  (5  centistokes)  . 
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specimen  temperatures  involved. 

Figure  3  shows  the  effective  heat  transfer  coefficient  for  the 
two  silicone  oils  as  a  function  of  the  film  temperature  calculated  by 
means  of  eq .  1  from  the  experimental  data  for  AT^  and  the  relevant 
material  properties.  In  these  calculations  the  strength  at  liquid 
nitrogen  temperature  (~196°C)  was  taken  to  be  194  MPa  as  measured  by 
Badaliance  et  al .  ^  and  its  temperature  dependence  was  assumed  to 
be  similar  to  that  in  reference  18.  The  strong  temperature  dependence  is 
evident.  Clearly  these  data  are  relevant  only  to  the  present  quenching 
experiments  and  probably  are  not  appropriate  to  other  heat  transfer 
conditions.  Nevertheless,  they  may  serve  as  a  guide  for  the  design  of 
quenching  experiments  of  other  materials. 

For  water  as  the  quenching  medium,  the  experimental  data  for  AT£ 

differed  significantly  from  those  obtained  for  silicone  oil.  Figure  4 

shows  the  dependence  of  ATc  on  the  water  bath  temperature  for  the  alumina 

rods  with  diameter  of  0.318  cm.  In  contrast  to  the  data  obtained  for  the 

silicone  oil,  ATc  for  the  water  bath  shows  an  upward  curvature  with  a 

minimum  near  60°C.  At  the  water  temperature  of  90°C,  only  a  small  fraction 

of  the  specimens  could  be  fractured  by  the  thermal  stress  generated  at 

quenching  temperature  differences  from  225  to  425°C.  Figure  5  shows  the 

distribution  of  the  percentage  of  specimens  fractured  during  the  quench 

as  a  function  of  AT.  No  specimens  could  be  fractured  for  values  of  AT  > 

425°C.  For  these  reasons,  for  the  bath  temperature  of  90°C,  a  value  for 

AT^  could  not  be  defined.  Comparison  of  the  data  in  Figs.  2  and  4  indicates 

that  AT  for  the  silicone  oil  exceeds  AT  for  the  water  bath,  due  to  large 
c  c 

differences  in  heat  transfer  coefficients  as  will  be  shown  later. 

Figure  6a  shows  the  size  dependence  of  AT^  for  the  alumina  specimens 
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g.  3.  Effective  heat  transfer  coefficient  for  sili¬ 
cone  oil  baths  inferred  from  the  thermal 
fracture  data  of  figures  la  and  lb. 


Effect  of  bath  temperature  on  ATC  for  alumina 
rods  (d  =  0.318  cm.)  quenched  in  water. 
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Fig.  5.  Probability  of  failure  of  alumina  rods  (d  - 
0.318  cm.)  quenched  in  water  at  90°C. 
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quenched  into  the  water  bath  at  25  and  90°C.  For  purpose  of  interpre¬ 
tation  of  the  data,  it  is  convenient  to  eliminate  the  effect  of  rod 
diameter  on  the  flexural  strength.  This  was  accomplished  as  shown  in 
fig.  6b  by  dividing  the  observed  values  of  AT^  by  the  ratio  of  the 
flexural  strength  relative  to  the  strength  for  the  rod  diameter  of 
0.318  cm,  listed  in  Table  2.  For  this  reason,  AT^  in  fig.  6b  is 
referred  to  as  "normalized"  critical  quenching  temperature  difference, 
ATc'.  Of  interest  to  note  is  that  AT  ^  decreases  with  increasing  value 
of  1/d,  in  contrast  to  an  increase  predicted  by  eq .  1.  Since  the 
derivation  of  eq,  1  was  based  on  theory,  this  discrepancy  cannot  be  at¬ 
tributed  to  a  lack  of  fundamental  understanding  of  thermo-mechanical 
principles. 

At  least  a  qualitative  explanation  for  the  data  shown  in  figs.  4  and 
6a  can  be  based  on  the  general  temperature  dependence  of  the  heat  trans¬ 
fer  coefficient  for  a  hot  surface  in  water.  This  is  shown  schematically 
in  fig.  7  as  a  function  of  the  difference  in  temperature  of  the  hot 
surface  and  the  water  temperature.  At  low  temperature  differences,  the 
heat  transfer  coefficient  increases  with  increasing  temperature  as  the 
direct  result  of  the  nature  of  convective  heat  transfer  coupled  with 
nucleate  boiling. ^  At  the  higher  temperature  differences,  the  rate  of 
heat  transfer  is  limited  by  the  formation  of  steam-film  on  the  hot  sur¬ 
face.  Because  of  the  existence  of  such  a  film,  the  value  of  the  heat 
transfer  coef  ficient  will  go  through  a  maximum  to  be  followed  by  a  de¬ 
crease  with  further  increases  in  temperature  difference.  At  sufficiently 
high  temperatures  where  radiation  from  the  hot  surface  becomes  signi¬ 
ficant,  an  increase  in  "effective"  heat  transfer  coefficient  can  be  ob- 


]  7 


served.  Increases  in  water  bath  temperature  towards  the  boiling  for  a 
given  hot  surface  temperature  will  enhance  the  formation  of  the  steam 
film  accompanied  by  a  corresponding  decrease  in  heat  transfer  coefficient. 

Examination  of  fig.  A  in  view  of  the  results  of  fig.  7  suggests 
that  the  minimum  in  the  curve  for  A Tc  vs  water-bath  temperature  corresponds 
to  the  maximum  in  the  curve  of  heat  transfer  coefficient  as  a  function 
temperature  difference.  The  data  in  fig.  6a  for  25°C  can  be  explained  in 
terms  of  the  rapid  rise  in  heat  transfer  coefficient  with  temperature 
difference  in  the  range  of  nucleate  boiling.  An  increase  in  AT^  automa¬ 
tically  increases  the  heat  transfer  coefficient  which  in  turn  limits  the 
increase  in  ATc  which  can  be  achieved.  In  fact,  the  plot  in  fig.  6b 
suggests  that  the  relative  rate  of  increase  in  heat  transfer  coefficient 
exceeds  the  rate  at  which  AT^  should  increase  with  decreasing  specimen 
size.  In  spite  of  these  effects,  ATc  increases  with  decreasing  specimen 
size  because  of  the  corresponding  increase  in  strength  with  decreasing 
specimen  size.  A  similar  explanation  can  be  given  for  the  data  for  the 
bath  temperature  of  90°C,  if  it  is  assumed  that  the  heat  transfer  coef¬ 
ficient  corresponds  to  the  low  values  near  the  minimum  between  the  regions 
of  steam-film  formation  and  radiation  indicated  in  fig.  7. 

Support  for  the  above  explanations  is  found  in  the  form  of  the  "ef¬ 
fective"  heat  transfer  coefficient  for  the  water  bath  calculated  from 
the  observed  values  of  ATc  (figs.  A  and  6a)  with  exception  of  those  for 
90°C  for  the  0.318  cm  diam.  rods,  and  the  property  data  for  the  alumina 
given  in  Table  2.  These  calculated  values  are  shown  in  fig.  8  as  a 
function  of  film  temperature  defined  earlier.  These  calculated  data 
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TEMPERATURE  DIFFERENCE,  AT 


Fig.  7.  Schematic  for  the  variation  of  the  heat 
transfer  coefficient  of  water  as  af¬ 
fected  by  the  solid  surface  tempera¬ 
ture. 


show  the  same  general  trends  as  those  shown  schematically  in  fig.  7. 

It  should  be  noted  that  these  calculated  data  vary  over  an  order  of  magni¬ 
tude.  Again,  it  is  cautioned  that  these  data  are  valid  only  for  the 
particular  testing  conditions  and  materials  of  the  present  study.  Un¬ 
less  verif  ied  ,  their  use  for  other  test  conditions  could  lead  to  misleading 
results . 

The  data  for  the  effective  heat  transfer  coefficient  for  the  water 
bath  shown  in  fig.  8  permit  a  reasonable  explanation  for  the  absence  of 
thermal  stress  failure  in  the  majority  of  specimens  with  diameter  of 
0.318  cm  quenched  into  water  at  90°C.  Because  of  their  small  size, 
as  indicated  by  eq.  1,  these  specimens  require  high  values  of  AT  to  induce 
thermal  stress  failure.  Such  high  values  of  AT  imply  high  film  tempera¬ 
ture  with  corresponding  low  values  of  the  heat  transfer  coefficients. 

The  few  specimens  in  which  failure  could  be  induced  represented  the  weakest 
specimens  of  the  strength  distribution  for  this  size  of  specimen.  It  is 
likely  that  if  the  strength  values  for  these  specimen  had  been  only 
slightly  higher,  failure  would  not  have  occurred  for  any  of  the  specimens. 
On  the  other  hand,  slight  reduction  in  the  strength  would  have  lowered 
ATc  with  a  corresponding  large  increase  in  heat  transfer  coefficient  that 
all  specimen  would  have  failed. 

Manson  and  Smith'*''’  found  for  disc  specimens  peripherally  quenched 
in  water,  that  the  thermal  stresses  were  nearly  constant  with  AT,  ex¬ 
hibiting  only  a  slight  concavity  with  a  maximum  near  ATc  =  120°C.  At 
least  qualitatively  this  suggests  that  the  heat  transfer  coefficients  for 
those  experiments  corresponded  to  the  range  of  film  temperature  shown  in 
fig.  8  which  includes  the  maximum  and  the  regime  over  which  h  decreases 
rapidly  with  increasing  film  temperature. 
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Becher  et  al.  calculated  the  effective  heat  transfer  coefficient 

from  ATc  for  a  number  of  ceramic  materials.  A  plot  of  these  values 
against  AT  resulted  in  a  curve  similar  to  the  one  shown  in  fig.  8. 

However,  if  the  values  of  h  calculated  from  the  present  data  were  plot¬ 
ted  in  a  similar  manner,  the  h  values  for  90°C  bath  temperature  would 
fall  to  the  left  of  the  data  for  the  lower  bath  temperatures.  This,  then, 
no  longer  would  result  in  the  expected  inverted  V-shaped  curve.  For 
this  reason,  it  appears  preferable  to  use  the  concept  of  film  temperature 
(or  perhaps  initial  specimen  temperature)  in  order  to  properly  reflect 
changes  in  heat  transfer  coefficient  due  to  differences  in  bath  tempera¬ 
ture. 

Davidge  et  al.7  found  that  on  increasing  the  water  bath  from  20 

to  95°C  resulted  in  a  decrease  in  AT  from  168  to  88°C  for  an  alumina 

c 

and  from  95  to  45°C  for  a  polycrystalline  magnesium  oxide.  These  ob¬ 
servations  appear  to  contradict  the  present  findings.  However,  as  judged 
by  the  relatively  low  values  of  ATc  for  the  samples  of  Davidge  et  al, 
(presumably  due  to  low  values  of  tensile  strength) ,  the  corresponding 
film  temperatures  for  bath  temperatures  of  20  as  well  as  95°C,  fall  in 
the  region  of  nucleate  boiling.  In  this  region,  the  heat  transfer  coef¬ 
ficient  ^increases  rapidly  with  increasing  film  temperature  so  that  a 
decrease  in  AT^  on  increasing  water  bath  temperature  is  reasonable. 

For  the  alumina  of  this  study,  however,  AT^  is  sufficiently  high  so  that 
the  heat  transfer  coefficient  decreases  with  increasing  value  of  AT. 
Consequently,  ATc  increases  with  increasing  bath  temperature. 
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4 .  CENERAL 

The  data  presented  in  this  paper  are  illustrative  of  the  dif¬ 
ficulties  encountered  in  establishing  the  relative  thermal  stress  re¬ 
sistance  of  ceramic  materials  by  the  quenching  method.  The  interpreta¬ 
tion  of  experimental  data  obtained  by  this  method  requires  a-priori 
quantitative  data  for  the  heat  transfer  characteristics  of  the  fluid 
media  used.  Such  data  are  not  easily  obtained.  Furthermore*  since  the 
specimen  continuously  changes  temperature,  the  heat  transfer  coefficient 
is  expected  to  vary  during  the  duration  of  the  quench.  Under  these 
conditions  a  calculation  of  the  thermal  stresses,  which  takes  into  ac¬ 
count  time-dependent  film  temperatures  and  corresponding  heat  transfer 
coefficients,  may  have  to  rely  on  numerical  methods.  Spatial  variations 
of  the  heat  transfer  coefficients  around  the  specimens  introduced  an 
additional  complexity. 

For  the  above  reasons,  the  heat  transfer  data  shown  in  figs.  3  and  8 

can  only  be  described  as  "effective"  since  they  were  inferred  from  the 

experimental  data.  Although  these  data  show  the  anticipated  dependence 

on  film  temperature,  their  general  applicability  to  other  materials  must 

be  treated  with  caution.  This  is  true,  in  particular,  for  the  water  bath, 

in  which  the  heat  transfer  coefficient  as  affected  by  nucleate  boiling  and 

steam  formation  is  expected  to  be  a  function  of  the  surface  condition 

(roughness,  etc.)  of  the  specimen.  It  is  even  conceivable  that  differences 

in  thermal  stress  resistance  of  different  materials  (i.e.  values  of  AT  ) 

c 

can  be  found  entirely  due  to  differences  in  surface  condition  which  modifies 
the  heat  transfer  characteristics. 
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Even  for  materials  with  surface  conditions  which  yield  identical 
heat  transfer  coefficients  for  a  given  film  temperature,  the  assessment 
of  the  relative  thermal  stress  resistance  based  on  data  for  AT^  is  not 
simple.  Consider  two  different  materials  with  relatively  low  thermal 
stress  resistance  so  that  at  ATc,  the  corresponding  film  tempera¬ 
ture  for  a  water  quench  is  in  the  nucleate  boiling  range.  The  material 
with  the  higher  thermal  stress  resistance  will  require  a  higher  value 
of  ATC.  This  however  will  involve  a  higher  film  temperature  and  corresponding 
higher  value  of  the  heat  transfer  coefficient.  This,  in  turn,  will 
suppress  the  value  of  AT^.  For  this  reason,  differences  in  relative 
thermal  stress  resistance  of  two  materials  subjected  to  a  water  quench 
involving  nucleate  boiling  may  be  difficult  to  assess. 

The  opposite  effect  occurs  for  two  materials  for  which  AT^  cor¬ 
responds  to  values  of  film  temperatures  which  involve  steam-film  forma¬ 
tion.  In  this  range  even  small  differences  in  AT  will  involve  large 
differences  in  heat  transfer  coefficient  and  corresponding  large  dif¬ 
ferences  in  ATc«  This  situation,  for  instance  could  arise  on  comparisons 
of  a  polycrystalline  aluminum  oxide  and  a  highly  thermal  shock  resistant 
material  such  as  high-strength  hot-pressed  silicon  nitride.  It  seems 
reasonable  to  point  out  that  water  despite  its  convenience,  appears  less 
useful  for  measurement  of  thermal  stress  resistance  by  the  quenching 
method . 

Silicone  oils,  on  the  other  hand,  appear  more  suitable  in  view  of 
their  monotonic  increase  of  the  heat  transfer  coefficient  with  increasing 
film  temperature.  For  these  oils  the  relatively  low  value  of  effective 
heat  transfer  coefficient  may  prevent  the  use  of  very  small  specimen  which 
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could  be  Inconvenient  especially  for  newly  developed  materials  for  which 
larger  specimens  may  not  be  available.  It  is  anticipated,  however, 
that  oils  can  be  found  or  developed  with  heat  transfer  characteristics 
higher  than  those  used  in  the  present  study. 

A  final  recommendation  is  in  order  for  purposes  of  facilitating 
quantitative  comparison  of  results  of  different  investigators.  In  re¬ 
porting  such  results,  the  conditions  of  the  quenching  tests  should  be 
stated.  This  includes  the  bath  temperature  as  well  as  the  range  of 
quenching  temperature  differences,  in  order  to  estimate  the  appropriate 
film  temperature.  In  view  of  the  rapid  variation  of  the  heat  transfer 
coefficient  with  film  temperature,  especially  for  water,  the  actual 
temperature  of  the  bath  should  be  stated.  Baths  at  nominal  room  tem¬ 
perature  does  not  appear  adequate,  in  view  of  large  regional  as  well 
as  seasonal  variations  of  this  quantity. 
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Figure  Captions: 


Fig.  1.  Effect  of  bath  temperature  on  ATC  for  glass  rods  quenched 

in  silicone  oil:  (a)  100  centistokes  (b)  5  centistokes. 

Fig.  2.  Effect  of  bath  temperature  on  ATC  for  alumina  rods  (d  =  0.63  cm.) 

quenched  in  silicone  oil  (5  centistokes) 

Fig.  3.  Effective  heat  transfer  coefficient  for  silicone  oil  baths 

inferred  from  the  thermal  fracture  data  of  figures  la  and  lb 

Fig.  A.  Effect  of  bath  temperature  on  ATC  for  alumina  rods  (d  —  0.318  cm.) 

quenched  in  water. 

Fig.  5  Probability  of  failure  of  alumina  rods  (d  =  0.318  cm.)  quenched 

in  water  at  90°C 

Fig.  6a.  Effect  of  specimen  size  on  ATC  for  alumina  rods  quenched  in 

water  at  25  and  90°C 

Fig.  6b.  Normalized  ATC'*  for  alumina  rods  quenched  in  water  at  25  and  90°C. 

Fig.  7.  Schematic  for  the  variation  of  the  heat  transfer  coefficient 

of  water  as  affected  by  the  solid  surface  temperature 

Fig.  8.  Effective  heat  transfer  coefficient  for  water  inferred  from 

the  experimental  ATC  values  in  figs.  4  and  6a. 
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ABSTRACT 

This  paper  discusses  the  effect  of  unstabilized  zirconia  dis¬ 
persions  on  the  thermal  shock  resistance  of  brittle  structural 
ceramics . 

A  general  discussion  is  given  of  the  preferred  direction  of 
modification  of  the  pertinent  material  properties  which  affect 
thermal  shock  resistance.  The  nature  of  the  crystallographic 
phase  transformation  in  the  zirconia  and  its  effect  on  these  Dro- 
perties  is  presented.  The  zirconia  dispersed  phase  can  lead  to 
significant  improvements  in  fracture  toughness  by  transformation  - 
or  microcrack  toughening.  The  volume  change  during  the  zirconia 
phase  transformation  also  lowers  the  effective  coefficient  of  ther¬ 
mal  expansion.  Surface  compressive  stresses  which  result  from  the 
phase  transformation  during  surface  grinding  also  are  beneficial 
in  improving  thermal  shock  resistance.  The  zirconia  phase  also  can 
change  the  unstable  (catastrophic)  mode  of  failure  to  the  more  pre¬ 
ferred  stable  one  with  decrease  in  fracture  stress. 

These  effects  are  illustrated  by  experimental  data  for  com¬ 
posites,  consisting  of  zirconia  dispersions  in  aluminum  oxide, 
silicon  nitride  and  zircon  matrices. 


INTRODUCTION 

The  advantageous  properties  of  ceramic  materials  such  as  their 
excellent  high-temperature  strength,  creep  and  corrosion  resis¬ 
tance,  make  them  ideal  candidates  for  the  replacement  of  such 


materials  as  metal  super-alloys  f>r  turbine  engines  and  many  other 
engineering  applications  involving  extreme  temperatures.  Unfortu¬ 
nately,  ceramic  materials  display  one  unfavorable  characteristic, 
namely  an  extreme  brittleness  or  notch-sensitivity .  This  high  de¬ 
gree  of  brittleness  renders  such  materials  highly  susceptible  to 
catastrophic  failure  under  conditions  of  thermal  shock.  Because  of 
the  absence  of  any  kind  of  plasticity,  the  non-uniform  thermal  ex¬ 
pansions  which  result  from  the  non-uniform  temperatures,  cannot  be 
accommodated  by  non-linear  deformation.  The  low  value  of  fracture 
toughness  of  brittle  ceramics  also  contributes  to  extensive  crack 
propagation  following  thermal  stress  failure.  This  frequently 
renders  any  structure  or  component  made  from  brittle  ceramics  com¬ 
pletely  incapable  of  continued  satisfactory  performance  following 
the  fracture  due  to  thermal  stresses. 

In  order  to  fully  capitalize  on  the  very  favorable  properties 
of  structural  ceramic  materials,  it  is  imperative  that  their  thermal 
shock  resistance  be  improved  significantly.  In  this  respect  it 
should  be  noted  that  ceramic  composites  exhibit  greater  thermal 
shock  resistance  than  single-phase  ceramic  materials  (1,2,3).  More 
recent  results  have  shown  that  the  introduction  of  an  unstabilized 
Zr02  dispersed  phase  in  a  brittle  matrix  can  lead  to  a  significant 
enhancement  of  many  mechanical  properties,  such  as  strength  and 
fracture  toughness,  of  importance  to  the  resistance  to  thermal 
stress  fracture.  This  effect  is  attributed  to  the  unique  tetrago- 
nal-to-monoclinic  phase  transformation  encountered  in  zirconia  (4- 
9). 


The  purpose  of  this  paper  is  to  summarize  the  role  of  this 
phase  transformation  in  enhancing  the  thermal  shock  resistance  of 
brittle  ceramic  materials. 


FUNDAMENTALS 

In  order  to  establish  the  role  of  the  zirconia  dispersions, 
a  brief  overview  will  be  given  of  the  pertinent  material  properties 
which  govern  the  thermal  stress  failure  of  brittle  materials. 


a.  Initiation  of  Thermal  Stress  Failure 

For  structures  or  components  in  which  thermal  stress  failure 
cannot  be  tolerated,  the  material  properties  for  a  given  design 
must  be  selected  such  that  the  maximum  thermal  stresses  do  not 
exceed  the  failure  stress  of  the  materials.  Generally,  brittle 
materials  exhibit  values  of  the  tensile  strength  well  below  the 
compressive  strength.  For  this  reason,  thermal  stress  failure  of 
brittle  materials  usually  occurs  in  tension. 

The  usual  procedure  followed  in  ceramic  technology  for  estab¬ 
lishing  the  role  of  the  pertinent  material  properties  which  affect 
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thermal  stress  failure  of  a  brittle  material  subjected  to  a  given 
thermal  environment,  is  to  derive  in  expression  for  the  maximum 
tensile  thermal  stresses  developed.  By  equating  this  value  of 
stress  to  the  tensile  strength,  an  expression  is  then  obtained  for 
the  maximum  thermal  environment  (i.e.,  temperature  difference,  heat 
flux,  etc.)  to  which  the  material  can  be  subjected  without  incur¬ 
ring  thermal  stress  failure.  As  an  example,  the  maximum  tempera¬ 
ture  difference  to  which  the  ceramic  structure  can.  be  subjected, 
can  be  expressed: 

Ao  j(l-v)K 


where  A  is  a  constant  depending  on  the  geometry  of  the  structure, 
ci f  is  the  tensile  fracture  stress,  v  is  Poisson's  ratio,  K  is  the 
thermal  conductivity,  uc  is  the  coefficient  of  thermal  expansion  and 
E  is  Young's  modulus. 

In  fracture-mechanical  terms,  the  tensile  strength  can  be  ex¬ 
pressed  : 

»f  ■  YKIc/ao1/2  <2> 

where  Y  is  a  constant  related  to  the  geometry  of  the  failure-ini¬ 
tiating  crack,  Kpc  is  for  critical  stress  intensity  factor  (frac¬ 
ture  toughness)  the  mode  I  crack-opening  and  ao  is  a  measure  of  the 
crack  size  or  depth. 


Substitution  of  eq.  2  into  eq.  1  yields: 

AYKt  ( 1-v) K 
.  _  ic 


If  in  the  region  of  the  failure  initiating  flaw,  a  compressive 
stress  (ac)  exists,  AT^ax  can  be  expressed  by: 

AYK  (l-v)K  Ao  (l-v)K 

AT  =— - —  +  - -  (4) 


Crack  Propagation  in  Thermal  Stress  Fields 


In  high-temperature  technology,  frequently  the  situation  is 
encountered  that  even  in  materials  with  the  highest  thermal  stress 
resistance,  thermal  stress  failure  cannot  be  avoided.  In  this 
case,  thermal  stress  resistance  is  not  measured  in  terms  of  the 
maximum  thermal  condition  to  which  the  materials  can  be  subjected 
without  failure.  Instead,  thermal  shock  resistance  is  measured 
by  the  ability  of  the  material  to  render  continued  satisfactory 


service  in  spite  of  having  undergone  thermal  stress  failure.  One 
such  popular  measure  for  this  criterion  is  the  strength  retained 
following  thermal  stress  failure.  The  retained  strength  is  a  mea¬ 
sure  of  the  size  of  the  arrested  cracks  which  resulted  from  the 
failure  process. 


A  simple  fracture-mechanical  analysis  of  the  nature  of  dyna¬ 
mic  crack  propagation  in  a  thermal  stress  field  has  shown  that  the 
size  of  the  crack  (af)  after  propagating  and  arrest  is  inversely 
proportional  to  a  ratio,  which  involves  the  energy  ( Y f )  required 
to  create  unit  area  of  new  crack  surface,  Young's  modulus  (E) , 
the  initial  fracture  stress  (of)  and  the  number  of  cracks  (N) , 
written  as: 


af  =  f(YfE/af2N)  (5) 

From  eq.  2,  the  strength  (aa)  retained  after  thermal  stress 
failure  is: 


o 


a 


YK 


Ic 


(6) 


From  eqs .  2,  5  and  6,  it  can  be  concluded  that  for  a  given 
number  of  cracks,  the  strength  retained  after  failure  is  an  inverse 
function  of  the  initial  strength  or: 

oa  =  f(of_1)  (7) 

In  terms  of  the  objectives  of  this  paper,  eqs.  1,  3  and  7,  indi¬ 
cate  the  problem  faced  by  the  materials  technologist  in  selecting 
or  developing  brittle  materials  subjected  to  high  values  of  thermal 
stress.  As  shown  by  eq.  3,  for  a  given  value  of  Kqc,  thermal 
stress  resistance  as  measured  by  ATmax  can  be  increased  by  de¬ 
creasing  the  flaw  size.  This,  in  effect,  increases  the  initial 
tensile  strength.  Unfortunately,  however,  as  indicated  by  eq .  7, 
if  this  increase  in  strength  is  not  sufficient  and  failure  still 
occurs,  the  only  effect  would  be  to  increase  the  size  of  the  cracks 
after  fracture  and  to  decrease  the  resulting  load-bearing  ability 
of  the  material.  In  fact,  for  ultra-high-strength  structural  cera¬ 
mic  materials,  thermal  stress  failure  can  be  quite  explosive,  ren¬ 
dering  the  structure  or  component  completely  useless  for  further 
service.  In  fact,  the  preferred  approach  would  be  to  increase  the 
original  crack  size.  This  has  the  further  advantage,  that  if  the 
flaw  size  becomes  large  enough,  crack  propagation  will  occur  in  a 
stable  manner.  At  the  same  time  high  strength  can  be  retained  both 
before  and  after  failure,  by  increasing  the  fracture  toughness  as 
much  as  possible.  In  addition,  reductions  in  the  coefficient  of 
thermal  expansion  and  Young's  modulus,  if  possible,  results  in  fur¬ 
ther  improved  thermal  stress  resistance. 


These  conclusions  are  summarized  in  Fig.  1,  which  shows  the 
typical  strength  behavior  following  thermal  shock  of  increasing 


severity.  The  discontinuity  of  uirength  at  ATC  is  the  result  of  the 
dynamic  crack  Propagation.  For  able  crack  proDagation,  no  such 
instantaneous  strength  loss  occurs  with  strength  decreasing  mono- 
tonical ly  for  values  of  AT  >  ATC .  Included  in  Fig.  1,  are  the  ra¬ 
tios  of  the  pertinent  material  variables  which  need  to  be  modified 
to  increase  thermal  stress  resistance.  Further  improvements  in 
thermal  stress  resistance  as  indicated  by  ecj.  4  can  also  be  achiev¬ 
ed  by  creating  a  compressive  stress,  which  opposes  the  tensile  ther¬ 
mal  stress,  in  the  region  of  the  material  where  the  failure  origi¬ 
nates. 


F.xperimenta  1  data  have  indicated  chat  incorporating  a  zirconia 
dispersed  phase  in  a  brittle  matrix  material  can  lead  to  signifi¬ 
cant  improvements  in  thermal  stress  resistance  by  increasing  the 
fracture  toughness,  lowering  the  coefficient  of  thermal  expansion, 
introducing  compressive  stresses  and  by  changing  the  failure  mode 
to  the  preferred  stable  mode  of  ctack  propagation  as  demonstrated 
in  the  following  sections. 


SEVERITY  OF  THERMAL  SHOCK 


Fig.  t.  Variation  of  strength  of  typical  brittle  ceramic  subject¬ 
ed  to  thermal  quench  of  increasing  severity  and  appropriate 
ratios  of  material  properties  relevant  to  thermal  stress 
resistance.  Solid  curve:  high-strength  ceramic;  dotted 
curve:  low  strength  ceramic. 
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MATURE  OF  ZIRCONIA  PHASE  TRAMS FOP MAT I ON 

In  order  to  understand  the  role  of  the  zirconia  dispersions  in 
improving  thermal  stress  resistance,  the  nature  of  the  tetragonal- 
monoclinic  phase  transformation  should  be  examined. 

Figure  2  shows  a  diagram  of  a  zirconia  particle  incorporated 
in  a  matrix  before  (t  =  tetragonal)  and  after  (m  =  monoclinic)  its 
transformation.  The  lower  half  of  Fig.  2  shows  the  thermal  expan¬ 
sion  behavior  of  a  brittle  matrix  containing  zirconia  inclusions. 

On  cooling  from  high  temperatures,  the  inclusions  go  through  the 
transformation  at  a  temperature  Tg,  which  involves  a  3-5%  volume 
expansion.  The  tangential  stresses  in  the  matrix,  caused  hy  this 
expansion,  can  lead  to  microcrack  formation  when  the  particle  dia¬ 
meter  exceeds  a  critical  value  dpc  (12).  On  reheating  the  parti¬ 
cles  can  retransform  to  the  tetragonal  phase  at  a  temperature  T^. 
Generally  Tg  >  Tg  so  that  the  expansion  behavior  shows  a  hysteresis 
The  height  of  the  hysteresis  curve  is  a  function  of  the  volume 
fraction  of  zirconia.  The  temperatures  Tg  and  Tg  are  a  function  of 


Fig.  2.  Zirconia  particle  contained  within  ceramic  matrix,  before 
(t)  and  after  (m)  the  tetragonal-monoclinic  phase  trans¬ 
formation  and  associated  thermal  expansion  behavior  of 
the  composite. 
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the  zirconia  particle  size  and  tin-  constraint  provided  by  the  mat¬ 
rix.  Tp  and  Tp  tenoral  Iv  decrease  with  decreasing  particle  5;  i  ze . 
Thus,  b v  controlling  the  particle  size,  tlie  nature  of  the  micro- 
cracking  in  the  matrix  as  well  a>  the  temperature  range  over  which 
the  transformation  occurs  can  be  controlled. 


EFFECT  OF  PHASE  TRANSFORMATION  IN  ZIRCONIA  ON  MATERIAL  PROPERTIES 
AND  STATE  OF  STRESS 

The  tetragonal-to-monocl inic  phase  transformation  in  the  zir- 
conia  dispersions  has  a  profound  ef  feet  on  a  number  of  material  pro¬ 
perties  which  control  thermal  stress  resistance  as  follows: 

1.  Fracture  toughness  (Rjc)  can  be  modified  by  two  methods  re¬ 
ferred  to  as  micro-crack  toughening  and  transformation  toughening. 

At  sufficiently  large  zirconia  particle  size,  the  hoop  stresses 
can  cause  extensive  micro-cracking  in  the  matrix.  During  failure, 
these  micro-cracks  promote  crack  bifurcation  and  an  enhanced  dis¬ 
sipation  of  energy,  thereby  effectively  increasing  the  fracture 
toughness.  Figure  3  indicates  the  increase  in  fracture  toughness 
which  can  be  achieved  in  this  manner  for  four  different  ceramic 
matrix  materials. 


Fig.  3.  Fracture  toughness  at  room  temperature  of  ceramic  matrix 

materials  without  (white  bars)  and  with  (grey  bars)  various 
volume  percentages  of  a  dispersed  zirconia  phase  indicated 
by  the  number  above  the  grey  bars. 
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Transformation  toughening  can  occur  by  zirconia  dispersions  of 
sufficiently  small  size,  which  as  the  result  of  the  matrix  con¬ 
straint  retain  their  tetragonal  crystal  structure  at  room  tempera¬ 
ture.  However,  during  fracture,  the  high  stresses  in  the  immediate 
region  of  the  crack  tip  cause  a  temporary  removal  of  the  constraint. 
This  permits  the  phase  transformation  to  take  place.  The  particles 
now,  however,  have  on  larger  volume,  so  that  a  passage  of  the  stress 
field  of  the  crack,  a, major  modification  in  the  matrix  stress  field 
is  encountered,  thereby  effectively  increasing  the  fracture  tough¬ 
ness  significantly. 

2.  Compressive  stresses.  The  zirconia  phase  transformation 
can  cause  the  creation  of  compressive  stresses  which  can  oppose 
the  tensile  thermal  stresses.  This  effect  can  be  achieved  by  two 
mechanisms  (16,17).  Steady-state  compressive  stresses  can  be 
created  by  the  stress-induced  phase  transformation.  One  very  ef¬ 
fective  method  by  which  this  can  be  achieved  is  to  cut  or  grind  the 
surface.  The  volume  expansion  associated  with  the  transformation 
leads  to  the  formation  of  high  compressive  stresses  in  the  immediate 
surface  regions.  Figure  4  shows  the  increase  in  bend  strength  and 
fracture  toughness  which  can  be  achieved  by  this  method.  Figure  4 


Fig.  4.  Effect  of  anneal  and  milling  time  on  flexural  strength  and 
fracture  toughness  of  hot-pressed  alumina  with  15  vol.X 
Zr02,  in  which  tetragonal-to-monoclinic  phase  transforma¬ 
tion  was  induced  by  surface  grinding  and  cutting. 
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also  contains  data  tor  the  eifeo.t  of  micro-crack  toughening  on 
strength  and  fracture  toughness,  labelled  "ground  and  annealed." 
Annealing  causes  a  relaxation  of  the  surface  stresses  so  that  any 
remaining  improvements  in  strength  and  toughness  must  be  attributed 
to  effects  other  than  surface  compression.  In  fact,  at  the  very 
short  (milling)  time  for  mixing  the  alumina  and  zirconia  powders, 
the  observed  increase  in  toughness  is  due  to  micro-crack  toughen¬ 
ing,  because  the  zirconia  particle  size  is  still  quite  large.  The 
decrease  in  fracture  toughness  with  increased  milling  time  is  due 
to  the  accompanying  decrease  in  particle  size. 

Thermal  stresses  arise  from  temperature  non-uniformities.  Be¬ 
cause  of  this  non-uniformity  of  temperature,  the  temperature  de¬ 
pendence  of  the  phase  transformation  permits  the  introduction  of 
transient  compressive  thermal  stresses  which  oppose  the  tensile 
thermal  stresses.  During  cooling,  for,  instance,  the  maximum 
thermal  stresses  occur  in  the  surface,  where  also  the  transformation 
will  occur  first.  The  volume  expansion  associated  with  this  trans¬ 
formation  created  a  compressive  stress  in  the  surface  which  opposes 
the  tensile  thermal  stress  due  to  the  greater  degree  of  cooling. 

It  should  be  noted  here  that  this  mechanism  for  improving  thermal 
stress  resistance,  is  effective  only  for  cooling  from  temperature 
immediately  above  or  within  the  range  of  transformation  tempera¬ 
ture.  A  similar  statement  must  be  made  for  compressive  stresses 
created  by  surface  grinding.  Such  stresses  are  useful  only  when 
the  tensile  thermal  stresses  exist  in  the  surface  alone  as  in  the 
case  for  transient  cooling. 

3.  Coefficient  of  thermal  expansion.  As  indicated  in  Fig.  2, 
the  volume  expansion  which  accompanies  the  phase  transformation  on 
cooling,  in  effect,  reduces  the  effective  coefficient  of  thermal 
expansion.  Since  the  range  of  transformation  temperature  is  a 
function  of  the  particle  size  incorporating  a  bi-modal  particle 
size  distribution  results  in  a  widening  of  this  temperature  range. 

In  fact,  a  coefficient  of  thermal  expansion  near  zero  can  be 
created  between  about  80°C  and  room  temperature.  This  effect  is 
indicated  schematically  in  Fig.  5. 

4.  Stable  crack  propagation.  When  a  high  density  of  micro¬ 
cracks  of  sufficient  size  is  introduced  in  a  matrix  such  that  its 
elastic  behavior  is  affected,  crack  propagation  in  a  thermal 
stress  field  will  become  stable.  Strength  loss  with  increasing 
severity  of  thermal  shock  will  be  monotonic  rather  than  discon¬ 
tinuous  for  unstable  crack  propagation  as  shown  in  Fig.  1.  Crack 
sizes  resulting  from  stable  crack  propagation  generally  are  smal¬ 
ler  than  those  which  result  from  unstable  crack  propagation.  For 
this  reason,  although  perhaps  initially  rather  weak,  materials 
which  fail  in  a  stable  manner  under  conditions  of  thermal  shock 
will  exhibit  superior  strength  behavior  following  the  thermal  shock. 
This  fact  is  used  to  advantage  in  refractory  structures  in  which 
thermal  stress  failure  generally  cannot  be  avoided,  but  in  which 
the  structural  materials  nevertheless  must  render  satisfactory 


Fig.  5.  Thermal  expansion  behavior  of  ceramic  matrix  with  zirconia 
inclusions  with  bi-modal  particle  size  distribition. 

service  in  spite  of  being  in  a  partially  fractured  condition. 

The  introduction  of  a  zirconia  dispersed  phase  can  also  be 
used  successfully  to  change  the  undesirable  unstable  failure  mode 
of  a  brittle  matrix  to  the  more  preferred  stable  one.  This,  how¬ 
ever,  will  require  relatively  large  particle  sizes  and  high  volume 
fractions  to  assure  a  significant  decrease  in  Young's  modulus,  a 
prerequisite  for  stable  crack  propagation.  As  an  example,  20  vol.% 
Zr02  particles  with  a  diameter  of  approximately  10  pm  in  an  alumina 
matrix  caused  a  decrease  in  Young's  modulus  and  bend  strength  from 
about  400  GN.irT^  and  500  MN.m-^  to  about  150  GN.m-^  and  170  MN.m~2. 
Of  interest  to  note  is.  that  the  microcracking  in  the  matrix  had 
little  or  no  effect  on  the  strain-at-fracture,  frequently  used  as 
a  measure  for  thermal  stress  resistance.  In  fact,  extensively 
microcracked  materials  can  exhibit  strains-at-f racture  an  order 
of  magnitude  above  that  for  crack-free  materials. 


THERMAL  SHOCK  DATA 

Experimental  data  were  obtained  to  illustrate  many  of  the  a- 
bove  phenomena  (18,19,20).  For  all  materials  studied,  specimens 


0 


were  subjected  to  thermal  shook  by  quenching  from  higher  temperature 
into  a  water  batli  at  room  temperature  followed  by  a  strength  de¬ 
termination  in  4-point  bending. 

Figure  6  shows  the  strength  behavior  of  an  alumina  matrix  with 
about  20  vol.%  Zr02  dispersions  of  various  sizes.  The  100%  alumina 
with  initial  strength  of  approximately  500  lfN.m-2  exhibited  a 
critical  quenching  temperature  difference  of  approximately  200°C. 

The  introduction  of  the  zirconia  particles  with  a  diameter  of  ap¬ 
proximately  1  pm,  as  the  result  of  transformation  toughening,  rais¬ 
ed  the  strength  to  about  650  MN.m~2  and  ATC  to  about  300°C,  a 
substantial  improvement.  Note  that  for  the  latter  two  materials 
as  indicated  by  the  instantaneous  strength  loss  at  ATC  fracture 
occurred  by  unstable  crack  propagation.  This  unstable  mode  of 
failure  can  be  changed  to  the  stable  one  by  incorporation  of  narti- 
cles  with  sizes  of  5  and  10  pm.  Mote  that  the  incorporation  of 
these  larger  particles  not  only  has  improved  ATC,  but  also  that 
the  strength  for  the  composites  for  both  particle  sizes,  at  high 
values  of  AT  exceeds  the  corresponding  values  for  the  100%  alumina 
as  well  as  the  alumina  with  the  1  pm  dispersions. 


AT  of  Thermoshock  ,  °C 


Fig.  6.  Effect  of  quenching  temperature  difference  on  strength  of 
composites  of  alumina  with  zirconia  dispersions  of  dif¬ 
ferent  sizes,  quenched  into  water  at  20°C. 
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As  an  illustration  of  the  decrease  in  the  coefficient  of  ther¬ 
mal  expansion  or  the  existence  of  a  transient  compressive  surface 
compressive  stress.  Fig.  7  compares  the  strength  behavior  of  silicon 
nitride  (lightly  dotted  curve)  and  silicon  nitride  with  20  vol.% 
Zr02,  subjected  to  a  similar  quench  as  the  alumina  of  Fig.  6.  The 
particle  size  of  the  zirconia  was  such  that  the  temperature  for  the 
phase  transformation  ranged  from  about  600  to  850°C.  Since  the  spe¬ 
cimens  were  at  room  temperature  prior  to  reheating  below  850°C  the 
zirconia  is  in  the  low  temperature  monoclinic  modification.  There¬ 
fore,  in  quenching  from  below  this  temperature,  no  advantage  is  ob¬ 
tained  from  the  phase  transformation.  In  Fig.  7  this  is  indicated 
by  the  data  points  in  the  form  of  the  crosses.  However,  on  heating 
to  above  850°C,  the  high-temperature  tetragonal  phase  is  obtained, 
which  upon  quenching  reverts  back  to  the  monoclinic  form  so  that 
thermal  stress  failure  is  avoided  as  indicated  in  Fig.  7.  Fracture 
by  quenching  from  600-850°C  can  be  avoided  only  if  the  specimens  are 


Fig.  7.  Thermal  expansion  and  retained  strength  after  thermal 
quench  of  silicon  nitride  with  zirconia  dispersions. 

For  meaning  of  different  data  points  the  reader  is  re¬ 
ferred  to  the  text.  Dotted  lines  indicate  strength  be¬ 
havior  of  silicon  nitride  without  zirconia  dispersions. 

first  heated  to  higher  temperature  to  achieve  the  tetragonal  crystal 
structure,  followed  by  slow  cooling  to  the  desired  temperature  from 
which  the  specimens  are  to  be  quenched.  In  Fig.  7  the  feasibility 
of  this  procedure  is  shown  by  the  data  points  indicated  by  circles. 
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Figure  8  shows  the  effect  of  the  ZrC>2  dispersions  on  the  ther¬ 
mal  shock  behavior  of  zirconium  silicate.  For  the  100%  zirconium 
silicate,  the  strength  loss  behavior  is  typical  of  unstahle  crack 
propagation  shown  in  Fig.  1.  The  effect  of  the  10  wt.%  zirconia 
dispersions  is  to  modify  the  fracture  mode  to  the  stable  one,  as 
indicated  by  the  monotonic  decrease  in  strength  for  AT  -  ATC.  This 
effect  must  be  atrributed  to  an  increase  in  the  size  of  the  failure 
initiating  flaws.  Of  interest  to  note  is  that  this  increase  in  the 
flaw  size  did  not  result  in  a  decrease  in  the  strength.  This  im¬ 
plies  that  the:  decrease  in  flaw  size  must  have  been  accompanied  by 
a  simultaneous  increase  in  fracture  toughness,  due  to  microcrack 
toughening.  It  is  clear  from  the  data  of  Fig.  8  that  the  zir¬ 
conium  silicate  with  zirconia  dispersions  is  much  preferred  over 
the  100%  zirconium  silicate  as  a  material  for  high-temperature 
structures  susceptible  to  thermal  stress  failure. 


Temperature  ,  °C 


Fig.  8.  Effect  of  quenching  temperature  difference  on  strength  of 
zirconium  silicate  with  and  without  zirconia  dispersions, 
quenched  into  water  at  20°C. 

CONCLUSION 

The  discussion  and  experimental  data  presented  in  this  paper 
indicate  that  the  incorporation  of  a  zirconia  dispersed  phase, 
can  lead  to  substantial  improvements  in  the  thermal  stress  resis¬ 
tance  of  brittle  structural  materials. 
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The  effect  of  compositional  and  microstructural  variables  and  pro¬ 
cessing  conditions  on  the  room  temperature  thermal  diffusivity  of  hot- 
pressed  and  reaction-sintered  silicon  nitride  was  determined.  The 
thermal  diffusivity  for  hot-pressed  silicon  nitride  increases  with  8- 
content.  Maximum  thermal  diffusivity  is  reached  at  about  3  wt.%  MgO. 
The  higher  thermal  diffusivity  of  the  8-phase  is  attributed  to  its 
higher  purity  level  and  the  less  distorted  crystal  structure  compared 
to  the  a-phase. 

In  reaction-sintered  nitride  the  thermal  diffusivity  is  strongly 
influenced  by  the  relative  amount  and  needle-like  morphology  of  the 
a-phase.  Correlations  of  the  thermal  diffusivity  with  mechanical 
properties  are  discussed. 


1 .  Introduction 

The  conduction  of  heat  in  dielectric  materials  depends  strongly 
on  the  density  and  nature  of  structural  lattice  imperfections  such  as 
dislocations,  vacancies,  lattice  vibrations  and  the  degree  of  crystal¬ 
linity  [l,  2,  3],  Microstructural  defects  such  as  pores,  grain  boundaries 
(at  low  temperature)  and  microcracks  also  have  a  significant  effect  on 
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heal  conduction  properties  [4,  5,  6].  Impurities  and  alloying  elements 
also  are  known  to  inhibit  the  rate  of  heat  conduction  in  dielectric 
materials  considerably  [4,  7]. 

Because  of  their  many  advantageous  properties,  dielectric  ceramics, 
such  as  silicon  nitride  and  others,  are  candidate  materials  for  high- 
temperature  energy  conversion  devices  such  as  the  all-ceramic  turbine 
or  in  the  form  of  replacements  for  many  metallic  components  in  internal 
combustion  engines.  The  performance  and  efficiency  of  such  systems 
in  addition  to  many  other  properties  depends  critically  on  the  thermal 
conductivity  and  thermal  diffusivity  of  the  materials  of  construction. 
High  engine  efficiency  frequently  requires  a  low  thermal  conductivity 
in  order  to  minimize  heat  losses.  On  the  other  hand  for  good  resis¬ 
tance  to  failure  due  to  thermal  stresses  high  values  of  the  thermal 
conductivity  and  diffusivity  are  preferred.  For  purposes  of  reliable 
engineering  design  and  the  achievement  of  the  proper  trade-offs  it 
is  essential  that  the  values  of  thermal  conductivity  and  thermal  dif¬ 
fusivity  of  such  dielectric  ceramics  and  the  variables  which  control 
them  are  well  understood. 

In  this  respect,  silicon  nitride  is  expected  to  be  very  complex. 
Silicon  nitride  in  the  hot-pressed  form  (HPSN)  can  exhibit  different 
densities  as  well  as  different  crystal  structures.  Hot-pressed  sili¬ 
con  nitride  can  consist  of  varying  amounts  of  the  a-phase  and  the  8- 
phase  which  are  bonded  together  by  a  grain  boundary  phase.  The  amount, 
composition,  crystallinity  of  the  grain  boundary  phase  are  expected  to 
play  a  major  role  in  establishing  the  values  for  the  heat  conduction 
properties.  Reaction-sintered  silicon  nitride  (RSSN)  can  exhibit 
variations  in  density,  free  unreacted  silicon  content,  different  crystal 
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structures  as  well  as  a  wide  range  of  pore  sizes  and  shape.  The  matrix 
of  reaction-sintered  silicon  nitride  consists  mainly  of  the  two  modi¬ 
fications  a  and  8  which  are  different  in  grain  size  and  grain  shape. 

The  a-matte  containing  very  small  pores  is  characterized  by  very  thin 
whisker-shaped  crystallites;  the  8-phase  exists  in  equiaxed  grain  [8,  9]. 
All  these  microstructural  variables  are  expected  to  effect  heat  conduction 
profoundly.  The  wide  variation  of  reported  values  of  the  thermal  con¬ 
ductivity  and  diffusivity  of  silicon  nitrides  including  CVD-Si^N^  and 
the  sialons  is  well  documented  [lO,  11,  12,  13,  14]. 

The  purpose  of  the  present  study  was  to  determine  the  thermal 
diffusivity  of  well-characterized  hot-pressed  and  reaction-sintered 
silicon  nitride  in  which  by  systematic  variations  of  grain  size  and 
composition  of  the  starting  powder  and  through  careful  process  control 
many  of  the  pertinent  compositional  and  microstructural  parameters 
were  isolated  as  much  as  possible. 

2 .  Experimental 

2.1  Sample  Preparation 

2.1.1  Hot-pressed  silicon  nitride 

The  microstructure  of  hot-pressed  silicon  nitride  can  be  controlled 
by  variation  of  phase  composition  and  morphology  of  the  starting  powder 
and  by  the  type  and  amount  of  sintering  aids  as  well  as  by  change  of 
the  processing  parameters  as  temperature,  pressing  time  and  pressure 
during  hot-pressing  process  [15,  16,  17,  18].  In  this  investigation 
MgO  was  used  as  sintering  aid.  To  study  the  effect  of  crystal  structure 
and  the  amount  of  glassy  phase  two  sets  of  specimens  were  prepared: 


3 


Set  I  : 


Variation  of  the  amount  of  0-*phase  at  constant  amount  of 
5  wt.%  MgO  sintering  aid.  These  specimens  were  hot--prpssed 
at  a  temperature  of  1973  K  for  various  lengths  of  time  to 
result  in  certain  percentage  of  the  3-phase. 

Set  II:  Variation  of  the  amount  of  MgO  content  between  1  and  10  wt.% 
for  constant  8-amount  of  100%.  Specimens  were  hot-pressed 
at  2053  K.  Hot-pressing  time  was  changed  in  such  a  way 
that  for  each  grade  with  different  MgO  content  complete  trans¬ 
formation  was  just  achieved. 

For  the  preparation  of  these  above  specimens  a  commercial  Si^Nj^ 

•k 

powder  was  used  with  a  0.04  fraction  of  8-phase  and  a  specific  surface 

area  of  8.8  m  /g.  The  impurity  contents  of  the  powder  were  (in  wt.%): 

0.02  Al;  0.04  Ca;  0.04  Fe;  0.6  C.  This  powder  was  mixed  with  different 

amounts  of  MgO  under  isopropanol  using  a  teflon  container  with  agate 

balls.  The  powders  were  hot-pressed  under  vacuum  at  a  pressure  of  31 
2 

MN/m  .  From  the  hot-pressed  discs,  75  mm  diameter,  6  mm  height,  circu¬ 
lar  (diameter  12.5  mm)  or  rectangular  plates  (5  mm  by  5  mm)  with  a  thickness 
of  2.1  mm  were  prepared.  These  samples  were  cut  in  such  a  way  that  the 
measured  thermal  diffusivity  corresponded  to  heat  flow  parallel  to 
the  hot-pressing  direction.  Both  sides  of  the  plates  were  diamond  ground. 

2.1.2  Reaction-sintered  silicon  nitride 

Three  sets  of  specimens  of  reaction-sintered  silicon  nitride  were 
prepared  following  procedures  described  elsewhere  [9,  19,  20] .  In  set  I 
artificial  pores  of  different  sizes  were  introduced  by  mixing  the  start¬ 
ing  silicon  powders  with  4  vol.%  of  an  organic  combustible  wax  of  dif- 

k 

H.  C.  Starck,  Goslar,  West-Germany . 
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ferent  particle  sizes  (0.36  pm;  63-90  pm;  125-180  pm).  These  mixtures 
were  isostatically  pressed  and  the  wax  particles  were  burnt  out  at 
653  K.  In  set  II  and  III  at  nearly  constant  total  porosity  the  pore 
structure  was  varied  by  using  different  grain  sizes  of  the  starting 
silicon  powder  (<  10  pm;  10-37  pm;  37-63  pm) .  These  specimens  were 
prepared  by  injection-molding.  The  specimens  of  all  three  sets  nitrided 
at  1673  K  for  100  hrs  in  a  mixture  of  90  vol.%  Nj  and  10  vol.%  H ^  at 
constant  pressure  of  950  mbar.  More  experimental  details  of  the  nitri- 
dation  procedure  are  given  by  [9]. 

2.2  Microstructural  and  Crystallographic  Characterization 

The  density  was  determined  by  the  Archimedes  principle  using  1^0 
at  20°C  for  HPSN  and  mercury  at  20°C  and  atmospheric  pressure  for  RSSN. 
Pore  and  grain  structure  of  HPSN  and  RSSN  were  examined  by  light  and 
scanning  electron  microscopy  on  unetched  and  on  etched  sections  and  on 
fracture  surfaces.  The  content  of  the  a-  and  8-phases  were  established 
by  quantitative  X-ray  crystallography.  The  free  Si  content  was  deter¬ 
mined  by  means  of  optical  metallography  and  by  X-ray  crystallography. 

For  HPSN  quantitative  grain  size  data  were  determined  from  polished 
and  etched  surfaces  (etching  with  HF  +  H2O2  +  HNO^  -  determination  of 
the  mean  grain  intercept  [l7,  18]).  For  RSSN  changes  in  the  distribution 
of  the  pore-opening  diameter  of  the  bottle-shaped  pores  (micropores)  were 
determined  using  mercury  pressure  porosimetry.  The  degree  of  nitridation 
in  RSSN  was  calculated  from  the  weight  grain  during  the  nitriding  pro¬ 
cess  [8],  The  morphology  of  the  a-phase  in  RSSN  was  determined  on  Ar 
ion-etched  microsections  [8]  and  on  unetched  fracture  surfaces. 
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During  the  hot-pressing  process  at  temperatures  above  1500°C 
a  liquid  phase  is  formed  and  the  a-crystals  dissolve  in  the  liquid 
and  the  3-crystals  are  recipitated.  Using  MgO  as  sintering  aid  the 
molten  magnesia  silicates  at  the  grain  boundaries  solidifies  during 
cooling  mainly  in  the  glassy  state.  The  thickness  of  the  grain  boundary 
phases  is  in  the  order  of  10-100  &  [ 17 ] .  During  the  a  to  &  transfor¬ 
mation  an  elongated  growth  of  6-crystals  takes  place  due  to  different 
growth  rates  perpendicular  and  parallel  to  the  c-axis  [l7].  The  mean 
grain  size  varies  from  0.3  to  1.5  pm  with  a  shortest  diameter  of  the 
rod-like  3-grains  between  0.2  and  2  pm  and  the  length  of  the  elongated 
grains  up  to  10  pm. 

Figure  1  shows  the  density  and  the  3/a+3  ratio  of  the  hot-pressed 
specimens  of  set  I  with  5  wt.%  MgO  as  a  function  of  hot-pressing  time 
at  1973  K.  The  densification  occurs  in  the  early  stage  of  the  hot- 
pressing  process.  As  expected  the  6-content  increases  with  increased 
duration  at  temperature.  The  kinetics  of  the  a  to  6  transformation 
are  strongly  influenced  by  the  amount  of  the  glassy  phase,  the  number 
of  the  pre-existing  6-grains  in  the  starting  powder  and  the  processing 
conditions  such  as  temperature  and  pressure  during  hot-pressing  [l5,  16,  17], 
With  regard  to  the  amount  of  the  sintering  aid  this  means  that  the  time 
required  for  complete  conversion  of  the  a-Si^N^  to  e-Si^N^  decreases  with 
increasing  MgO  content. 

Figure  2  shows  the  change  in  microstructure  during  a  to  8  trans¬ 
formation  in  set  I.  After  short  hot-pressing  time  a  fine  grained  mostly 
equiaxed  grain  structure  is  observed.  With  increasing  hot-pressing 
time  a  change  in  grain  morphology  from  equiaxed  to  prismatic  rod-like 
grains  takes  place  accompanied  by  the  growth  of  elongated  grains. 


For  the  specimens  in  Set  II,  Fig,  3  gives  the  density  as  a  function 

of  MgO  content  for  hot-pressing  times  required  for  a  100%  conversion  of 

the  a-  into  the  6-phase.  The  density  indicates  a  slight  decrease  with 

increasing  MgO  content.  The  deviations  from  theoretical  density  of 
3 

3.20  g-cm  probably  is  the  result  of  the  relatively  high  amount  of 
the  grain  boundary  glassy  phase.  For  all  MgO  contents  the  microstructure 
consists  of  elongated  6-grains.  The  values  of  the  mean  grain  intercepts 
are  in  the  range  of  0.7  to  1.2  ym.  No  systematic  changes  in  grain  size 
were  observed. 

2.2.2  Reaction  sintered  silicon  nitride 

Table  1  lists  the  data  for  the  microstructural  characterization 
for  the  RSSN  samples.  As  indicated  by  Morgan  [2l]  a  number  of  questions 
on  the  nature  of  the  formation  of  the  a-  and  6-phase  in  RSSN  still  remain 
to  be  answered.  For  the  method  of  preparation  used  for  the  specimens  of 
this  study,  the  a-phase  most  probably  is  formed  by  a  gas-phase  reaction 
between  the  Si  and/or  SiO  and  N£,  which  leads  to  whisker-shaped  crystal¬ 
lites  in  the  form  of  a  matte.  At  nitriding  temperatures  below  the  melt¬ 
ing  point  of  Si,  the  6-phase  is  formed  by  a  diffusion  process  of  the 
nitrogen  to  the  Si-Si^N^  interface  [8,  9],  As  indicated  by  the  data 
in  Table  1,  the  size  of  the  natural  pores  In  sets  II  and  III  range 
from  about  0.06  to  0.25  ym  for  the  micropores  and  from  7  to  34  ym  for 
the  macropores.  The  micropores  are  thought  to  correspond  to  the  spacings 
between  the  very  line  fibers  in  the  a-matte  [8].  The  6-grains  were 
found  to  be  equiaxed  ranging  in  size  from  about  0,5  to  3  ym.  Further 
microstructural  details  for  these  materials  were  presented  in  earlier 
publications  [8,  19], 
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In  set  I  at  nearly  constant  other  microstructural  parameters  only 
the  size  of  the  artificially  introduced  macropores  are  changed.  In 
the  sets  II  and  III  the  amount  of  a-phase  decreases  from  batches  A 
to  C,  while  the  diameter  of  the  micro-  and  macro-pores  increase.  The 
change  in  the  a/a+0  ratio  and  in  the  size  of  micro-  and  macropores  is 
due  to  the  variation  of  the  grain  size  of  the  starting  silicon  powder 
and  to  the  corresponding  change  in  size  and  arrangement  of  pores  in 
the  green  body  [9].  Figure  4  shows  the  microstructure  of  the  a-phase 
for  the  two  specimens  A  and  C  of  set  III,  which  exhibited  large  dif¬ 
ferences  in  thermal  diffusivity  for  the  same  density.  These  micro¬ 
graphs  clearly  indicate  the  difference  in  morphology,  in  particular 
the  increase  in  thickness  and  in  contact  area  between  the  fibers  in 
batch  C. 

2 . 3  Technique  for  Measuring  the  Thermal  Diffusivity 

The  thermal  diffusivity  was  measured  by  the  laser-flash  technique 
[ 1 1 ] .  Samples  of  the  silicon  nitride  in  the  form  of  thin  rectangular 
or  circular  plates  were  flashed  on  one  side  by  means  of  a  50  J  glass- 
Nd  laser.  A  liquid  Nj  cooled  infra-red  detector  coupled  to  a  storage- 
oscilloscope  was  used  to  monitor  and  record  the  transient  temperature 
of  the  opposite  side  of  the  sample.  The  recorded  trace  of  the  tempera¬ 
ture  as  a  function  of  time  permitted  calculation  of  the  thermal  dif¬ 
fusivity  [22].  Measurements  were  made  only  at  room  temperature,  at 
which  any  differences  in  the  thermal  diffusivity  are  expected  to  be  most 
pronounced  without  having  the  resort  to  specimen  cooling  and  major 
equipment  modifications. 
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2 . 4  Elastic  Properties,  Strength  and  Fracture  Toughness 

Frequently  for  porous  materials  correlation  are  found  between 
the  amount  of  pore  phase,  pore  morphology  and  the  thermal  and 
mechanical  properties.  For  this  reason,  experimental  data  for  the 
mechanical  behavior  can  assist  in  the  proper  interpretation  of  the 
thermal  property  data.  For  the  samples  of  reaction-sintered  silicon 
nitride  of  the  present  study,  the  elastic  properties,  fracture  stress 
(determined  in  4-point  bending  tests)  and  fracture  toughness  were  re¬ 
ported  in  detail  as  part  of  a  previous  study  [l9].  The  results  are 
shown  in  Table  2. 

3.  Experimental  Results  and  Discussion 

3.1  Hot-Pressed  Silicon  Nitride 

Figure  5  shows  the  thermal  diffusivity  of  the  silicon  nitride 
with  5%  MgO  as  a  function  of  hot-pressing  time.  As  indicated  in 
Fig.  1,  increasing  time  for  hot-pressing  results  in  the  transformation 
of  the  a-  into  the  8-phase.  Comparisons  of  Fig.  1  and  5  indicate 
clearly  that  the  thermal  diffusivity  shows  a  significant  increase 
with  increasing  B-content. 

Some  explanations  for  this  observation  can  be  given.  For  very 
short  times  of  hot-pressing  (t  =  10  mm)  the  increase  can  be  attributed 
to  an  increase  in  density.  For  the  samples  subjected  to  longer  hot- 
pressing  times  in  which  complete  densification  occurred  the  increase 
in  thermal  diffusivity  must  be  attributed  to  changes  in  the  relative 
amounts  of  the  a-  and  8-phase  and  to  corresponding  changes  in  grain 
morphology  and  heat  transfer  hehavior . 
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It  is  expected  that  heat  conduction  in  the  a-  and  3-silicon 
nitride  at  room  temperature  occurs  primarily  by  phonon  transport. 

For  this  reason,  differences  in  the  thermal  diffusivity  must  arise 
because  of  differences  in  the  phonon  velocity  or  phonon  mean  free 
path.  Since  the  elastic  properties  and  density  of  a-  and  B-Si^N^ 
do  not  differ  significantly  (as  far  as  the  authors  are  aware)  a  dif¬ 
ference  in  thermal  diffusivity  due  to  differences  in  phonon  velocity 
seems  unlikely.  Phonon  mean  free  path  differences  can  be  attributed 
to  at  least  two  possible  effects.  Firstly,  differences  appear  to  exist 
in  the  crystal  structure  of  the  a-  and  B-modifications. 

Henderson  and  Taylor  [23]  suggested  that  the  a-phase  has  a  more 
strained  crystal  structure  than  the  3-phase.  For  this  reason,  phonon 
scattering  is  expected  to  occur  to  a  greater  degree  in  the  a-phase 
than  in  the  3-phase,  thereby  yielding  a  lower  value  for  the  thermal 
diffusivity  of  the  former. 

A  second  reason  for  this  increase  in  thermal  diffusivity  with 
increasing  B-content  may  be  related  to  the  formation  and  growth  mechanism 
of  the  B-Si^N^.  Although  still  a  topic  of  discussion  [24],  the  forma¬ 
tion  of  the  B-phase  is  thought  to  result  from  the  dissolution  of  the 
a-phase  in  the  liquid  grain  boundary  phase  and  a  re-precipitation  of 
the  silicon  nitride  in  the  8-form  [l7].  It  is  likely  that  the  solubility 
of  impurities  in  the  B-phase  is  lower  than  in  the  liquid  grain  boundary 
grain  boundary  phase.  For  this  reason,  it  is  expected  that  during  the 
a-  to  8-  transformation  an  enrichment  of  the  purities  in  the  grain 
boundary  phase  occurs,  so  that  the  purity  of  the  B-grains  is  higher 
than  the  original  a-grains.  Impurities  or  other  atoms  are  well-known 


10 


to  be  effective  phonon  scatterers  [4,  7]  so  the  higher  purity  of  the 
6-phase  is  responsible  for  its  higher  diffuslvity. 

The  difference  in  the  size,  shape  and  orientation  of  the  6-grains 
compared  to  the  a-grains  should  also  be  examined  as  possible  contributing 
factors  to  the  observed  increase  in  thermal  diffusivity.  Since  the 
phonon  mean  free  path  in  the  silicon  nitride  at  room  temperature  is 
much  smaller  than  the  grain  size  any  effect  of  the  latter  appears  un¬ 
likely.  A  change  in  grain  shape  for  randomly  oriented  grains  contained 
in  a  matrix,  as  may  be  inferred  from  the  analytical  results  of  Bruggeman 
[25]  and  Niesel  [26]  should  not  lead  to  a  change  in  thermal  diffusivity. 
However,  as  shown  by  Powers  [27]  from  the  results  of  Fricke  [28]  an 
elongation  of  the  grains  in  a  preferred  orientation  will  result  in  an 
increase  or  decrease  in  the  direction  of  heat  flow  parallel  and  per¬ 
pendicular  to  the  direction  of  preferred  orientation,  resp.,  regardless 
of  the  relative  values  of  the  thermal  conductivities  of  the  two  phases. 
However,  since  the  6-grains  are  elongated  perpendicular  to  the  hot- 
pressing  direction  [29]  and  the  direction  of  heat  flow,  the  present 
observation  cannot  be  attributed  to  this  latter  effect.  Therefore,  it 
can  be  concluded  that  the  increase  in  thermal  diffusivity  with  in¬ 
creasing  B-content  after  complete  densification  must  be  because  the 
thermal  conductivity/dif fusivity  of  the  6-phase  is  inherently  higher 
than  of  the  a-phase. 

Figure  7  shows  the  data  for  the  thermal  diffusivity  as  a  function 
of  MgO  content  after  complete  a-to-6  transformation.  As  can  be  noted, 
with  increasing  MgO  content,  the  thermal  diffusivity  first  increases  to 
a  maximum  followed  hy  a  decrease  at  the  highest  values  of  MgO  content. 

The  initial  increase  in  thermal  diffusivity  at  low  MgO  content  is  thought 


to  be  due  primarily  to  a  reduction  of  residual  porosity  with  increasing 
MgO  content.  Maximum  thermal  diffusivity  is  reached  at  about  V/ 

MgO.  This  MgO-content  is  commonly  used  in  practice  to  achieve  optimum 
mechanical  behavior.  The  decrease  in  thermal  diffusivity  with  higher 
MgO  content  most  likely  can  be  attributed  to  the  increasing  amount  of 
glassy  grain  boundary  phase.  Glassy  materials  generally  have  values  of 
thermal  conductivity  or  diffusivity  some  order  of  magnitude  less  than 
crystalline  solids.  For  high  MgO  contents  the  higher  amount  of  magnesium 
silicates  is  expected  to  be  particularly  effective  in  controlling  the 
overall  rate  of  heat  transfer  through  hot-pressed  MgO-Si^N^. 

3 . 2  Reaction-Sintered  Silicon  Nitride 

Table  2  lists  the  experimental  data  for  the  thermal  diffusivity, 
elastic  properties,  strength  and  fracture  toughness  for  all  specimens 
sets  I,  II,  ITT.  Comparisons  with  the  microstructural  and  compositional 
data  in  Table  1  leads  to  the  following  correlations  and  conclusions: 

The  data  for  the  specimens  for  set  I,  which  were  prepared  to 
ascertain  the  effect  of  macropore  size  on  thermal  mechanical  properties 
near  constant  other  microstructural  variables  showed  that  no  such  pore 
size  effect  exists  for  the  thermal  diffusivity  and  Young's  modulus,  in 
agreement  with  theory  [30].  The  decrease  in  strength  with  increasing 
macro-pore  size  is  in  accordance  with  theoretical  considerations  based 
on  pore  size  to  flow  size  ratios  [31,  32], 

The  data  for  the  specimens  in  sets  II  and  III  indicate  a  pro¬ 
nounced  negative  correlation  between  the  thermal  diffusivity  and  a/ct+B 
ratio  and  a  positive  correlation  with  the  size  of  the  micro-  and  macro¬ 
pores.  As  indicated  by  the  data  for  set  I,  this  variation  in  thermal 
diffusivity  cannot  be  attributed  to  the  variation  in  size  of  the  macro- 
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pores.  The  negative  correlation  between  the  strength  and  the  size  of 
the  macropores  is  to  be  expected  from  pore  size  to  flow  size  considera¬ 
tions,  referred  to  earlier. 

Particularly  noteworthy  for  the  data  for  sets  IT  and  III  is  the 
strong  negative  correlation  between  the  thermal  diffusivity  and  the 
strength  at  nearly  constant  Young’s  modulus. 

The  above  observations  and  correlations  are  not  in  general  agree¬ 
ment  with  the  behavior  of  single-phase  brittle  structural  materials. 

For  these,  with  decreasing  density,  strength  shows  the  greatest  relative 
decrease  followed  by  Young's  modulus,  thermal  conductivity  and  thermal 
diffusivity  in  accordance  with  theoretical  predictions  [33-37],  Also, 
in  the  absence  of  radiation  heat  transfer,  thermal  conductivity  and 
diffusivity  should  be  unaffected  by  pore  size. 

An  answer  to  these  observations  must  be  found  in  the  unique  micro¬ 
structure  of  reaction  sintered  silicon  nitride.  Fibrous  microstructures, 
due  to  the  long  path-length  for  phonon  heat  transfer  will  exhibit  a  much 
lower  thermal  conductivity  [38]  than  a  matrix  with  the  same  relative 
density  containing  isolated  equi-dimensional  pores.  For  this  reason, 
the  thermal  conductivity/diffusivity  of  RSSN  is  expected  to  decrease 
with  increasing  a-content.  Furthermore,  differences  in  fiber  morphology 
for  a  given  a-content  are  expected  to  affect  the  thermal  diffusivity 
as  well.  Generally,  the  smaller  the  thickness  and  inter-spacing  of  the 
fibers  for  constant  bulk-density  the  lower  will  be  the  thermal  con¬ 
ductivity,  in  part  due  to  increase  phonon  scattering  at  the  fiber 
surfaces,  reduced  contact  areas  between  the  fihers  and  reduction  in 
heat  transport  by  the  gaseous  phase  and  radiation  within  the  pores. 

For  this  reason,  the  observed  variations  in  the  thermal  diffusivity  for 
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present  RSSN  samples  are  due  primarily  to  differences  in  morphology 
and  the  relative  amount  of  the  a-phase.  This  effect  is  enhanced 
further  because  as  stated  earlier,  the  lower  thermal  diffusivity  of 
the  a-phase  compared  to  the  6-phase  is  thought  to  be  due  to  the  more 
distorted  crystal  structure  of  the  former. 

The  inverse  correlation  between  strength  and  thermal  diffusivity 
at  near  constant  Young's  modulus,  particularly  noticable  for  the 
specimens  of  set  II,  probably  also  can  be  attributed  to  the  amount 
and  size  of  the  a-fibers.  An  increase  in  fiber  amount  is  expected  to 
improve  the  fracture  toughness,  as  observed,  which  improves  strength. 

The  free  Si-content  does  not  seem  to  have  a  pronounced  effect 
on  the  thermal  diffusivity,  possibly  because  its  relatively  small  amount. 

3 . 3  General 

it  should  be  noted  that  although  for  both  the  hot-pressed  and 
reaction-sintered  silicon  nitride,  the  thermal  diffusivity  increases 
with  increasing  6-content;  these  effects  occur  for  different  reasons. 

For  the  hot-pressed  silicon  nitride,  the  more  perfect  crystal  structure 
and  the  higher  purity  of  the  6-phase  leads  to  the  increase  in  the  thermal 
diffusivity.  In  the  reaction-sintered  silicon  nitride,  however,  the 
fibrous  morphology  and  amount  of  the  more  strained  a-silicon  nitride 
appears  to  be  the  deciding  factor  in  establishing  the  value  of  the 
thermal  diffusivity. 

The  results  of  the  present  study  suggest  that  by  careful  process 
modifications,  a  measure  of  control  can  be  exerted  over  the  heat  con¬ 
duction  properties  of  silicon  nitride.  For  high  thermal  stress  resistance 
which  requires  high  thermal  conductivity  the  relative  amount  of  the  6- 
phase  should  be  maximized.  For  reaction-sintered  silicon  nitride,  however. 
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this  latter  approach  involves  an  undesirable  trade-off  with  a  decreasing 


strength.  On  the  other  hand,  for  engine  components  which  require  low 
thermal  conductivity  in  order  to  maintain  high  engine  efficiency 
because  of  reduced  heat  losses,  the  relative  a-content  should  be 
kept  as  high  as  possible  for  both  the  hot-pressed  as  well  as  the  reac¬ 
tion-sintered  silicon  nitride. 
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Table  1.  Microstructural  and  Compositional  Characteristics 
of  Reaction  Sintered  Si.N. 
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Table  2.  Thermal  Diffusivity,  Tensile  Strength, 
Elastic  Properties  and  Fracture  Tough- 
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Fig,  2,  Change  in  grain  morphology  with  increasing  3/a+6  ratio  in  silicon 
nitride  with  5  wt,%  MgO  hot  pressed  at  1973K  for:  a,  10  min;  b, 

40  min  and  c,  100  min. 
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HOT-PRESSING  TIME  tmln) 


Thermal  diffusivity  at  room  temperature  of  hot- 
pressed  silicon  nitride  with  5%  MgO  (set  I)  as 
a  function  of  hot-pressing  time  at  1973°K. 
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Thermal  diffusivity  at  room  temperature  of  silicon 
nitride  hot-pressed  at  2053°K  (set  II)  with  8/a+6= 
1  as  a  function  of  MgO  content. 
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The  hot-pressing  of  silicon  nitride  requires  a  densif ication  aid 
such  as  magnesium  oxide.  During  hot-pressing  the  grains  of  the  silicon 
nitride  in  the  a-modification  are  transformed  gradually  into  elongated 
grains  of  the  3-modification  most  probably  by  a  dissolution  and  re¬ 
precipitation  process.  The  3-grains  exhibit  a  preferred  orientation 
such  that  their  direction  of  elongation,  which  coincides  with  the 

crystallographic  c-axis,  is  oriented  perpendicularly  to  the  hot-pressing 
1-4 

direction  .  This  preferred  orientation  has  been  shown  to  result  ir. 
a  pronounced  anisotropy  in  strength  at  high  temperature  as  well  as  frac¬ 
ture  toughness.  Recent  results^  indicate  that  the  a-to-3  transformation 
resulted  in  a  significant  increase  in  the  room  temperature  thermal  diffusi- 
vity  of  hot-pressed  silicon  nitride  in  a  direction  parallel  to  the  hot- 
pressing  direction.  It  was  suggested  that  this  latter  effect  arose  because 
the  thermal  conductivity  of  the  3-phase  was  higher  than  for  the  a-phase. 

Theoretical  solutions  of  Powers^  based  on  the  original  theory  of 
Fricke^,  for  the  thermal  conductivity  of  composites  consisting  of  matrices 
with  oriented  elongated  dispersions  in  the  shape  of  prolate  ellipsoids 
of  revolution  indicate  that  the  thermal  conductivity  of  the  composite 
would  be  a  maximum  in  the  direction  of  elongation.  This  effect  would 
occur  regardless  whether  the  matrix  or  dispersions  had  the  higher  thermal 
conductivity.  At  least  qualitatively  this  result  is  expected  to  occur 
for  hot-pressed  silicon  nitride,  due  to  the  elongation  of  the  3-grains, 
perpendicular  to  the  hot-pressing  direction.  For  this  reason,  the  thermal 
conductivity  and  diffusivity  perpendicular  to  the  hot-pressing  direction 
would  be  expected  to  be  higher  than  the  corresponding  values  parallel  to 
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the  hot-pressing  direction.  The  purpose  of  this  communication  is  to 
present  experimental  data  in  support  of  this  conclusion. 

For  this  purpose  specimens  were  cut  from  disks  approximately  70  mm 

wide  by  6  mm  thick  prepared  by  vacuum  hot-pressing  in  graphite  dies  at  a 
-2  5 

pressure  of  35  MN-m  for  previous  studies  on  the  effect  of  micro¬ 
structure,  composition  on  mechanical  properties  and  thermal  diffusivity"* 
measured  parallel  to  the  hot-pressing  direction.  For  the  specimens,  the 
processing  parameters,  MgO  content,  density  and  B/ct+8  ratio  are  listed 
in  Table  1.  The  aspect  ratio  of  the  8-grains  was  about  5.  Additional 
discs  containing  10  wt.%  MgO  were  hot-pressed  as  well  for  a  range  of  time 
periods  as  listed  in  Table  1.  Specimens  for  the  measurement  of  the 
thermal  diffusivity  in  the  form  of  plates  approximately  5  mm  square  by  2  mm 
thick  were  cut  from  larger  specimens  used  for  measurements  of  bend  strength. 
Most  specimens  for  thermal  diffusivity  were  cut  from  the  same  strength 
specimens  in  order  to  minimize  possible  effects  of  non-uniformity  in  density 
or  composition  within  the  hot-pressed  disks. 

9 

The  thermal  diffusivity  was  measured  by  the  laser-flash  technique 
using  equipment  described  elsewhere^.  The  transient  temperature  of  the 
specimen  during  the  measurement  was  monitored  by  a  liquid  Ng  cooled  IR- 
detector.  Data  were  obtained  at  room  temperature  at  which  any  directional 
effects  on  the  thermal  diffusivity  would  be  expected  to  be  a  maximum  without 
cooling  the  specimens  below  room  temperature. 

The  experimental  data  are  shown  in  Figs.  1,  2  and  3.  Comparison 
shows  that  for  all  compositions  and  heat  treatments  the  thermal  dif¬ 
fusivity  perpendicular  to  the  hot-pressing  direction 
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exceeds  the  corresponding  value  parallel  to  the  hot-pressing  direc» 
tlon,  in  agreement  with  the  above  hypothesis.  The  observation  that 
the  thermal  diffusivity  increases  with  increasing  8-amount  for  both 
directions  relative  to  the  hot-pressing  direction  is  confirmation 
that  the  thermal  conductivity/dif fusivity  of  the  8-phase  is  signifi¬ 
cantly  higher  than  that  of  the  a-phase.  This  was  attributed^  to  the 
fact  that  the  crystal  structure  of  the  a-phase  is  a  more  distorted 
structure  with  shorter  phonor.  mean  free  path  than  the  8-structure  and 
probably  also  to  impurity  segregation  in  the  grain  boundary  phase  so 
that  the  8-grains  have  lower  impurity  content  than  the  original  a- 
phase.  The  maximum  in  the  thermal  dif fusivity  as  a  function  of  MgO 
content  probably  occurs  because  once  complete  densi f ication  near  3  wt . % 
MgO  is  achieved,  further  additions  of  MgO  result  in  a  lowering  of  the 
thermal  diffusivity  as  the  result  of  the  lower  value  of  the  glassy 
grain-boundary  phase.  The  slight  increase  in  the  thermal  diffusivity 
with  increasing  hot-pressing  time  for  10%  MgO  content  possibly  is  due  to 
an  increase  in  crystalline  order.  X-ray  studies  of  this  effect  were 
not  positive,  possible  as  the  result  of  the  low  thickness  of  the  grain 
boundary  phase  (<  100  8) . 

In  general,  the  present  results  indicate  that  even  at  a  single 
temperature,  the  heat  transport  properties  of  hot-pressed  silicon 
nitride  are  subject  to  a  number  of  simultaneous  variables  including 
the  effect  of  composition,  microstructure,  processing,  purity,  as 
well  as  orientation.  For  purposes  of  practical  application  of  sili¬ 
con  nitride  it  is  imperative  that  such  effects  are  well  understood. 
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COMPOSITION,  PROCESSING  CONDITIONS  AND  DENSITY  OF  HOT-PRESSED  SILICON  NITRIDE 


THERMAL  OIFFUSIVITY  (cm's")  do  THERMAL  DIFFUSIVITY  (cm's 


.  1.  Effect  of  hot-pressing  direction  on  the  thermal 
diffusivity  of  silicon  nitride  subjected  to 
various  hot-pressing  time  to  increase  8-fraction , 


Fig.  2.  Effect  of  hot-pressing  direction  on  the  thermal 
diffusivity  of  silicon  nitride  with  a  range  of 
MgO  content  for  p/a+8  ratio  near  unity. 
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Fig.  3.  Effect  of  hot-pressing  direction  on  the  thermal 
diffusivity  of  silicon  nitride  with  10  wt.t  MgO 
and  8/i+B  =  1  as  a  function  of  hot-pressing  time 


CHAPTER  V 


OBSERVATIONS  ON  THE  SUB-CRITICAL  GROWTH  AND  HEALING  OF  MICRO¬ 
CRACKS  IN  BRITTLE  CERAMICS 
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Department  of  Materials  Engineering 
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Blacksburg,  Virginia  24061 


ABSTRACT 

Micro-cracking  in  brittle  composites  was  monitored  by  measur¬ 
ing  the  temperature  dependence  of  thermal  diffusivitv  by  the  laser 
flash  method.  Depending  upon  the  material  system,  micro-cracks  can 
exhibit  a  time  dependent  growth  or  healing  or  a  combination  of  both. 
A  theoretical  basis  for  these  observations  was  established  by 
analyzing  the  stability  and  nature  of  crack  propagation  of  pre¬ 
cursor  micro-cracks  in  a  spherical  inclusion  contained  in  an  in¬ 
finite  matrix  with  different  elastic  properties. 


INTRODUCTION 

Microcracked  materials,  in  view  of  their  potential  for  high 
fracture  toughness* ,  high  thermal  shock  resistance-  and  low  ther¬ 
mal  conductivity*,  are  of  considerable  technical  interest.  Opti¬ 
mization  of  these  properties  requires  a  detailed  understanding  of 
the  variables  which  control  the  formation  of  microcracks. 

Micro-cracking  Is  observed  in  composites  with  large  mismatches 
in  the  coefficients  of  thermal  expansion  of  the  components.  Ex¬ 
tensive  micro-crack ing  can  also  occur  in  polycrvstal 1 inc  aggregates 
with  high  thermal  expansion  anisotropy  of  the  individual  grains. 

Mi rro-r rack i ng  occurs  due  to  the  internal  stresses  which  arise  from 
the  non-uniform  thermal  expansion  whenever  t lie  material  is  subject¬ 
ed  to  a  temperature  change  from  the  value  at  which  the  material  is 
stress-free.  This  latter  temperature  usually  is  the  fabrication 
temperature.  For  linearly  clastic  materials  the  magnitude  of  the 
internal  stresses  is  a  function  of  the  thermal  expansion  differ- 


i 


3 


A 


lin  os,  the  range  of  temperature  change  and  the  elastic  properties 
of  the  components^  or  the  polycrystalline  array.  The  magnitude  of 
the  internal  stresses  can  also  be  affected  by  stress  relaxation 
phenomena  as  governed  by  the  kinetics  of  available  creep  mechanisms 
and  the  time  period  over  which  such  relaxation  is  permitted  to  oc¬ 
cur  . 


Microcracking  is  also  strongly  dependent  on  grain  size  ,  in 
spite  of  the  fact  that  the  magnitude  of  the  internal  stresses  is 
independent  of  the  scale  of  the  microstructure.  This  observation 
has  led  to  the  development  of  a  theory  of  microcrack  formation  bas¬ 
ed  on  the  criterion  that  microcracking  will  occur  only  for  grains 
of  sufficient  size  such  that  the  elastic  energy  of  the  residual 
stress  field  equals  the  surface  energy  required  to  create  a  crack 
equal  to  the  grain  size*.  This  hypothesis  does  not  define  the  con¬ 
ditions  for  the  onset  of  microcrack  formation. 

More  recently,  the  conditions  for  the  formation  of  micro-cracks 
have  been  defined  in  terms  of  the  size  and  geometry  of  micro-crack 
precursors  in  the  form  of  crack-like  residual  pores  at  triple  points 
and  grain  boundaries^.  From  these  precursors,  microcracks  will 
form  when  the  stress  intensity  factor  exceeds  the  critical  stress 
intensity  factor.  Since  the  size  of  such  precursors  frequently  is 
proportional  to  the  grain  size,  an  indirect  grain  size  effect  is 
expected . 

The  variables  which  control  microcracking  can  be  studied  by 
the  measurement  of  those  material  properties  which  are  affected 
by  microcracks.  These  include  strength,  fracture  toughness,  elas¬ 
tic  behavior  and  thermal  conductivity  or  diffusivity.  These  latter 
two  properties  are  preferred  since  their  measurement  does  not  in¬ 
volve  the  application  of  a  mechanical  load,  which  could  affect 
crack  stability. 

For  a  material  with  randomly  oriented,  equal  sized,  penny¬ 
shaped  cracks,  the  thermal  conductivity  (K)  In  the  absence  of  heat 
transfer  across  the  cracks  is-*: 

K  =  K  [l  +  8tU3/9]"1  (1) 

o 

where  K0  is  the  thermal  conductivity  of  the  crack-free  material, 

N  is  the  number  of  cracks  per  unit  volume  and  P.  Is  the  radius  of 
the  cracks.  Since  microcracking  has  negligible  effect  on  the 
specific  heat  per  unit  volume,  eq.  1  also  describes  the  effect  of 
micro-cracking  on  the  thermal  diffusivity. 

Equation  1  indicates  that  any  changes  in  the  thermal  conduc¬ 
tivity  or  diffusivity  can  be  attributed  to  changes  in  the  crack 
density  or  crack  size  or  combinations  thereof.  For  this  reason, 
measurements  of  the  thermal  diffusivity  of  microcrackod  material 
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can  provide  information  on  the  varini  Us  viiUli  control  microcrack 
stability  and  propagation  behavior.  The  purpose  of  this  paper  is 
to  present  a  number  of  such  experimental  data.  These  will  be  inter¬ 
preted  qualitatively  first,  followed  by  a  quantitative  analysis. 

EXPERIMENTAL 

The  experimental  data  for  the  thermal  diffusivity  to  be  pre¬ 
sented  for  a  number  of  different  materials  include  recent  observa¬ 
tions  as  well,  as  data  obtained  in  a  number  of  previous  studies. 

These  latter  studies  contain  descriptions  of  the  equipment  and  pro¬ 
cedures  of  the  laser-flash  diffusivity  method  used  for  these  mea¬ 
surements  as  well  as  details  on  material  fabrication  and  micro¬ 
structures. 

Figure  1  shows  the  thermal  diffusivity^  as  a  function  of  tem¬ 
perature  for  a  coarso-pra ined  pot ycrystal I i ne  l rnn-i i Innate.  Be¬ 
cause  of  its  relatively  low  values  of  thermal  diffusivity  compared 
to  a  fine-grained  iron-titanate,  the  material  of  fig.  1  must  have 
undergone  extensive  micro-cracking  as  it  was  cooled  from  the  an¬ 
nealing  temperature  prior  to  the  measurement  of  the  thermal  diffu¬ 
sivity.  It  is  significant  to  the  objectives  of  the  present  study 
that  the  thermal  diffusivity  displays  a  significant  hysteresis  in 
such  a  manner  that  the  thermal  diffusivity  on  cooling  considerably 
exceeds  the  corresponding  values  obtained  during  heating.  To  pro¬ 
vide  an  explanation  for  this  behavior  it  should  be  noted  that  on 
heating  the  material  the  differential  thermal  strains  responsible 
for  the  microcracking  are  reduced  or  even  eliminated.  At  suffi¬ 
ciently  high  temperature  the  excess  surface  energy  of  the  micro¬ 
cracks  can  bo  reduced  by  material  transport  bv  diffusional  process¬ 
es,  This  leads  to  crack-healing  i.e.,  a  decrease  in  crack  size  and 
possible  crack  density  and  a  corresponding  increase  in  the  thermal 
di f  f usl vitv . 

Figure  2  shows  the  thermal  diffusivity  of  a  f ine-gra ined 
magnesium  dftitanate  also  subjected  to  a  single  heating  and  cool¬ 
ing  cycle.  In  contrast  to  the  behavior  noted  in  fig.  1,  these 
data  show  that  the  thermal  diffusivity  on  cooling  has  values  below 
those  found  on  the  heating  part  of  the  cycle.  This  suggests  that 
during  the  thermal  cycle  an  increase  in  crack  densitv  and/or  crack 
size  occurred. 

It  is  critical  to  note  that  those  data  suggest  that  the  micro¬ 
cracks  formed  during  cooling  from  the  fabrication  or  annealing  tem¬ 
perature  are  not  necessarily  in  an  equilibrium  condition.  As  an 
explanation  for  the  data  of  fig.  2  it  is  suggested  that  during  the 
first  cooling,  the  precursor  microcracks  were  not  of  sufficient 
size  that  they  could  lead  to  microcrack  Formation  at  room  tempera¬ 
ture.  On  reheating  the  specimens,  however,  they  could  grow 
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fig.  1 .  Therma 1  di f  f u- 
sivity  of  thermal ly  cycled 
coarse-grained  iron-tita- 
nate. 
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Fig.  2.  Thermal  diffusi¬ 
vity  of  thermally  cycled 
fine  grained  magnesium 
dititanate. 
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snhrr  i  t  ical  1  y  to  a  latpe  enough  si7.e,  such  that  on  cooling  to  room 
temperature,  criticallitv  was  reached  and  microcracks  formed. 

Figure  3  shows  the  thermal  diffusivity  for  magnesium  ditita- 
natc®*  which  indicates  a  combination  of  the  behavior  noted  for  figs. 
1  and  2.  At  high  temperatures  the  expected  crack  healing  occurs 
but  on  cooling  to  room  temperature  the  values  for  the  thermal  dif¬ 
fusivity  fall  below  the  values  at  the  beginning  of  the  thermal 
cycle.  This  behavior  is  attributed  to  the  existence  of  a  relative- 
Iv  large  size  distribution  of  tin  microcrack  precursors.  During 
cooling  from  the  annealing  temnerature  those  precursors  which  were 
largo  enough  led  to  the  formation  of  microcracks.  However,  the 
precursors  which  were  not  of  sufficient  size,  required  a  second 
reheat  for  them  to  grow,  which  then  led  to  additional  microcracking 
during  the  second  cooling,  t.o.,  the  cooling  part  of  the  thermal 
cycle  for  the  measurement  of  the  thermal  diffusivity.  Again,  as 
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Fig.  3.  Thermal  diffusi- 
vity  of  thermally  eye  1  oil 
coarse-grained  magnesium 
dititanate  exhibiting  both 
crack  healing  and  growth. 
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stated  earlier,  microcrack  configurations  developed  during  a  single 
cycle  may  not  represent  an  equilibrium.  In  general,  the  number 
and  size  of  the  microcracks  and  their  effect  on  material  properties 
are  a  function  of  the  thermal  history  of  the  material.  Furthermore, 
theoretical  assessments  of  microcracking  should  include  various 
statistical  aspects  of  microcrack-precursors  including  not  only 
their  size,  but  also  their  orientation  and  geometry. 

Further  support  for  the  above  conclusions  is  given  in  figs. 

4 a  and  4b,  which  show  the  thermal  diffusivitv  of  composites  of 
magnesium  oxide  containing  20  wt . X  SiC  cycled  to  1300<’C  and  1400°C, 
respectively.  both  these  figures  show  qualitatively  similar  beha¬ 
vior  to  that  shown  in  fig.  3.  However,  the  specimen  cycled  to 
1400°C  shows  a  thermal  diffusivitv  on  cooling  which  is  higher  than 
for  tlie  specimen  cycled  to  1300°C.  These  differences  can  be  attri¬ 
buted  to  the  differences  in  the  temperature  dependence  of  the  pro¬ 
cesses  for  microcrack  healing  and  precursor  growth.  At  1300°C,  the 
interna]  stresses  will  be  higher  than  at  1400°C.  For  this  reason, 
at  the  lower  temperature  the  material  will  exhibit  a  higher  ten¬ 
dency  for  precursor  growth  than  at  the  higher  temperature.  The 
opposite  effect  will  exist  for  crack-healing,  which  is  more  likely 
to  occur  at  1400°C  than  at  1300°C.  For  these  reasons,  heating  to 
1400°C  is  expected  to  show  a  greater  improvement  in  thermal  dif- 
fusivity  on  cooling  than  heating  to  1300nC,  in  agreement  with  the 
observations.  This  effect  is  further  evidenced  by  the  diffusivity 
data  for  MgO  +  20  wt.7  SiC  composites  which  showed  an  8.1 7,  11.37., 
and  15.47.  decrease  in  room  temperature  thermal  diffusivity  values 
for  specimens  thermally  cycled  to  1300°C  and  held  for  0,  3,  and 
10  hours  respectively. 
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Fig.  4n.  Thermal  diffusivity  of 
magnesium  oxide  with  20  wt.%  sili¬ 
con  carbide  thermally  cycled  to 
1 100°C. 


Fig.  4b.  Thermal  diffusivity 
of  magnesium  oxide  with  20  wt.% 
silicon  carbide  thermally  cveled 
to  1400°C. 


Similar  effects  are  indicated  in  fig.  5,  which  shows  the 
values  for  the  thermal  diffusivity  at  room  temperature  before  and 
after  thermal  cycling  of  SJC-MgO  composites  as  a  function  of  SiC 
content.  The  decrease  in  the  thermal  diffusivity  after  thermal 
cycling  for  all  compositions  is  evidence  for  the  growth  and/or 
formation  of  additional  cracks.  It  Is  of  interest  to  note  that  a 
minimum  in  the  thermal  diffusivity  occurs  at  about  20  wt.%  SiC. 

This  implies  that  at  this  composition  microcracking  is  most  exten¬ 
sive.  At  first  sight,  it  may  be  expected  that  the  number  of  micro- 
cracks  should  be  proportional  to  the  amount  of  silicon  carbide, 
with  a  corresponding  monotonic  decrease  in  the  thermn’  diffusivity. 
The  fact  that  this  is  not  observed  is  thought  to  be  related,  as 
described  in  greater  detail  elsewhere",  to  the  change  in  phase 
distribution.  At  low  SiC  contents,  the  silicon  carbide  represents 
a  discontinuous  dispersed  phase.  At  the  higher  SiC  contents,  how¬ 
ever,  part i c le-to-part icle  contact  will  occur  so  that  both  the  SiC 
and  Mgfl  become  continuous  phases.  This  will  have  a  major  effect 
on  tlie  magnitude  anti  distribution  of  the  internal  stress  field, 
which  can  have  a  profound  effect  on  the  onset  of  microcracks  due 
to  orientation  effects  of  the  precursors,  as  well  the  ability  of 
the  composite  to  undergo  internal  stress  relaxation. 

Evidence  of  the  stability  of  microcracks  controlled  by  still 
another  phenomenon  is  given  in  fig.  6,  which  shows  the  thermal  dif¬ 
fusivity  of  a  relatively  fine-grained  iron-titanate  subjected  to 
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Fig.  5.  Effect  of  thermal  cycle  to 
1400°C  on  room  temperature  thermal 
diffusivity  of  MgO-SiC  composites. 

repeated  thermal  cycles^.  These  data  indicate  that  repeated  ther¬ 
mal  cycling  leads  to  permanent  recovery  of  the  thermal  diffusivity. 
This  indicates  that  following  the  first  cooling  from  the  manufac¬ 
turing  or  annealing  temperature,  microcracks  exhibit  maximum  dimen¬ 
sions,  which  are  permanently  decreased  following  subsequent  heat 
treatments.  This  effect  is  expected  for  micro-cracks  formed  from 
relatively  small  precursors.  In  analogy  to  crack  formation  under 
conditions  of  thermal  shock*0  these  cracks  can  form  by  unstable  pro¬ 
pagation  to  reach  a  size  such  that  they  are  subcritical  with  res¬ 
pect  to  the  temperature  condition  required  for  their  creation. 
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Fig.  6.  Effect  of  temperature  and  thermal 
cycling  on  the  thermal  diffusivity  of  nnlv- 
crystallin.e  FcTTiOr  annealed  for  4  hrs.  at 
1100°. 
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Such  cracks  should  exhibit  permanent  healing  especially  if  partially 
healed,  so  that  on  subsequent  e\  I es  repropagation  occurs  in  a 
stable  mode  only. 

The  collective  evidence  o!  Lhese  experimental  data  is  clearly 
indicative  of  the  non-equilibrium  nature  of  microcracks.  In  parti¬ 
cular,  their  size  and  numbt r  can  be  strongly  influenced  by  sub- 
critical  growth  and  healing  effects.  A  quantitative  basis  for 
these  effects  can  be  provided  by  a  simple  analysis  of  microcrack 
stability  and  propagation  behavior. 


ANALYSIS 

The  mechanical  model  chosen  for  the  analysis  consists  of  a 
single  isotropic  inclusion  contained  within  an  infinite  isotropic 
matrix.  Microcracking  is  assumed  to  occur  within  the  inclusion. 
Since  the  stress  (strain)  state  within  the  inclusion  is  uniform,  by 
treating  the  relevant  properties  of  the  inclusion  with  cracks  as 
pseudo-cont inuum  properties,  analytical  solutions  are  easily  ob¬ 
tained.  Estimates  of  micro-crack  formation  at  the  interface  or 
within  c lie  matrix  where  the  stress  state  is  non-uniform  require 
numerical  methods.  However,  the  analytical  results  to  be  obtained 
for  microcracking  within  the  inclusion,  also  qualitatively  describe 
behavior  of  microcracking  at  the  interface  or  within  the  matrix. 

The  analysis,  similar  to  the  one  by  Griffith,  is  based  on  the 
hypothesis  that  micro-cracks  will  be  unstable  if  for  propagation 
the  decrease  in  the  elastic  energy  of  the  internal  stress  field 
associated  with  the  inclusion  equals  or  exceeds  the  corresponding 
increase  in  surface  energy.  In  contrast  to  the  Griffith  theory, 
this  analysis  assumes  the  existence  of  two  surface  energy  levels. 
The  first  one  corresponds  to  the  minimum  energy  level  required  for 
suhcritical  crack  growth.  The  second  (higher)  one  corresponds  to 
the  value  of  surface  energy  required  for  fast  fracture. 

For  an  isostatic  stress  o  in  the  inclusion,  the  total  elastic 
energy  (W)  within  the  inclusion  and  matrix  isll: 

W  =  2tto2(1-2v  )R3/F,  +  tto2(1+v  ) R3/E  (2) 

p  p  m  m 

where  v  is  Poisson's  ratio,  E  is  Young's  modulus,  R  is  the  inclu¬ 
sion  radius  and  p  and  m  refer  to  the  inclusion  and  matrix  respec¬ 
tively. 

In  terms  of  the  misfit-strain  due  to  the  mismatch  of  the  co¬ 
efficients  of  thermal  expansion  of  the  inclusion  and  matrix,  W  can 
be  written: 

W  =  2ffA«2AT2R3/((l-2vp)/Ep  +  (l+vm)/2F.m)  (3) 
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where  An  is  the  difference  in  coefficients  of  thermal  expansion 
between  matrix  and  inclusion  and  AT  is  the  cooling  range. 

For  randomly  oriented  penny-shaped  microcracks  of  equal  size. 
Young's  modulus  Ep  of  the  inclusion  can  be  expressed: 

F.  =  F  [l  +  16(l-\>2  )N?.3/9(l-2v  )]_1  (4) 

P  po  P  P  J 

where  N  is  the  number  of  cracks  per  unit  volume  of  inclusion,  Fp0 
is  the  Young's  modulus  of  the  inclusion  without  cracks  and  l  is  the 
radius  of  the  cracks.  The  condition  2.<R  is  assumed  to  exist  always, 
for  the  validity  of  the  mechanical  model. 

The  total  surface  energy  of  the  cracks  within  the  particle 
is: 

T  =  4  tt  R  3  Ny  tt  C, 2  /  3  (5) 

where  y  is  the  surface  energy.  * 

Crack  instability  occurs  when 

d/d£(W-T)  >  0  (6) 

Substitution  of  eqs.  3  and  5  into  eq .  6  followed  by  differen¬ 
tiation,  leads  to  the  critical  temperature  difference  ATC  based  on 
the  fracture  surface  energy  Yf  for  unstable  crack  propagation: 

f7rYf ( l-2v  )  ,  cE  (1+v  )  16(l-v  2)NjA 

at  *  MU _ E _ J  V'5  JE2 _ DL_  +  I  +  _ 0 _ _ 

c  (Aaj2E  (i-v2T?J  2E  (l-2v  )  1  9(l-2v  ) 

pop  mp  P 

(7) 

The  expression  for  the  critical  temperature  difference  ATSC 
based  on  the  fracture  surface  energy  ys  for  subcritical  crack  growth 
can  be  obtained  from  equation  7  by  using  ys  in  place  of  yf. 

Small  cracks  of  length  1.Q  which  on  instability  at  AT  =  ATC 
propagate  in  an  unstable  mode  when  d2(W+T)/dS.2  >  0,  will  arrest 
when  the  total  elastic  strain  energy  released  is  converted  to  frac¬ 
ture  surface  energy.  This  will  occur  at  a  final  crack  length 

l.  =  (3  E  (1+v  ) / 16  E  ( 1— v  )Nft  Y'2  (8) 

t  po  m  m  p  o 

Equations  7  and  8  are  not  conveniently  expressed  in  dimension¬ 
less  form.  For  this  reason,  ATc,  ATSC  and  If  were  calculated  for 
an  MgO  inclusion  within  a  SiC  matrix  for  a  value  of  crack  density 
N  =  10fl/unit  volume  and  a  ratio  Yf/yst.  =  4.  The  solid  curves  In 
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fig.  7  denote  AT  and  AT  .  Similarly,  is  shown  by  the  dotted 
sc  c  J  '  f  1 

curve . 

Figure  7  shows  three  regimes  of  different  crack  behavior.  In 
the  regime  labeLled  (A),  the  crack  will  propagate  critically.  Sub- 
critical  crack  growth  occurs  in  region  B,  whereas  region  C  indicates 
crack  healing.  A  number  of  conclusions  can  be  made  from  fig.  7 
related  to  the  experimental  data  presented  earlier. 

It  is  critical  to  note  that  a  change  in  AT  also  involves  a 
change  in  the  absolute  value  of  temperature.  This  in  turn,  can  have 
a  large  effect  on  the  rate  of  the  diffusional  processes  responsible 
for  crack  growth  and  crack  healing.  For  this  reason,  on  cooling 
to  low  temperatures  microcracks  are  not  expected  to  exhibit  their 
equilibrium  length,  because  crack  growth  and  healing  is  suppressed. 
In  other  words,  microcracks  can  be  "super-cooled".  Only  at  tem¬ 
peratures  where  rates  of  crack  growth  or  healing  are  sufficiently 
high,  can  microcracks  reach  their  equi librium  length  corresponding 
to  these  particular  values  of  AT  and  temperature.  This  implies 
that  microcracked  materials  which  in  service  are  subjected  to 
variations  in  temperatures,  generally  will  not  be  at  equilibrium 
and  are  expected  to  show  temperature  and  time-dependent  properties. 

A  number  of  characteristics  of  microcrack  stability  and  pro¬ 
pagation  can  be  discussed  on  the  basis  of  the  numerical  data  in 
fig.  7.  For  small  values  of  crack  length  such  that  in  eq.  7  the 
term  involving  the  crack  size  is  small  compared  to  unity,  ATC  is 
independent  of  crack  density  and  inclusion  size.  In  addition,  eq . 

8  clearly  indicates  that  for  a  given  crack  density,  the  final 
crack  length  which  results  from  unstable  crack  propagation  at 
AT  =  ATC  also  is  independent  of  the  inclusion  size.  This  suggests 
that  theories  for  the  existence  of  a  critical  grain  size  based  on 
energy  considerations  will  require  reexamination. 

Inclusion  size  will  have  an  effect  on  ATC  when  for  a  given 
number  of  cracks,  the  cracks  are  sufficiently  large  such  that  the 
effective  Young's  modulus  of  the  inclusion  is  affected.  This  be¬ 
havior  is  responsible  for  the  increase  in  ATSC  and  ATC  with  in¬ 
creasing  crack  length  to  the  right  of  the  minimum  of  the  solid 
curves  in  fig.  7. 

A  microcracked  material,  reheated  to  a  temperature  close  or 
equal  to  its  fabrication  or  annealing  temperature  (AT  -*•  0),  will 
undergo  crack  healing.  This  should  manifest  itself  in  an  increase 
in  the  thermal  diffusivity,  in  agreement  with  the  data  shown  in 
figs.  1,  3  and  4. 

The  value  of  final  crack  length  (denoted  by  if"  in  fig.  7) 
which  results  from  instability  of  a  short  crack  of  length  H0"  at 


Fir,.  7.  Microcrack  stability  diagram  in¬ 
dicating  critical  and  sub-critical  crack 
propagation  and  crack  healing  for  three 
cooling  histories  (numerical  data  for 
MgO  inclusion  in  SiC  matrix). 

ATC,  is  subcritical  with  respect  to  ATSC.  Such  a  crack  for  any 
value  of  AT  £  ATSC  will  undergo  permanent  crack  healing.  For  this 
type  of  crack  behavior,  the  crack  size  will  be  largest  after  cool¬ 
ing  from  the  preparation  temperature,  when  the  thermal  diffusivitv 
also  will  show  its  minimum.  Due  to  crack  healing,  heating  will 
cause  a  permanent  increase  in  the  thermal  diffusivity,  as  shown  by 
the  data  in  fig.  6. 

Materials  with  a  range  of  sizes  of  precursors  will  undergo 
partial  microcracking  due  to  those  precursors  which  are  of  critical 
size.  On  reheating,  the  microcracks  will  undergo  healing  whereas 
the  precursors  which  were  too  small  for  microcrackinp  will  undergo 
crack  growth.  On  a  second  cooling,  those  precursors  which  grew  to 
critical  size  will  propagate  into  microcracks.  This  behavior  mani¬ 
fests  itself  in  a  recovery  in  the  thermal  diffusivity  at  higher 
temperatures  and  a  room  temperature  diffusivity  value  which  is 
lower  at  the  end  of  the  thermal  cycle  than  at  the  beginning.  This 
behavior  is  in  agreement  with  the  data  shown  in  figs.  3  and  A.  The 
relative  rates  of  crack  healing  and  precursor  growth  are  expected 
to  be  very  sensitive  to  temperature  in  view  of  the  competitive  ef¬ 
fects  of  the  magnitude  of  temperature  and  residual  stress  on  the 
kinetics  of  crack  healing  and  growth.  If  on  the  first  cooling, 
none  of  the  precursors  were  of  sufficient  size  to  create  micro¬ 
cracks,  but  grew  on  reheating,  thermal  diffusivity  as  a  function 
of  temperature  as  shown  in  fig.  2  would  be  expected. 

The  above  discussion  suggests  that  many  of  the  experimental 
observations  can  be  explained  by  the  analytical  results  for  a 
simple  mechanical  model  for  microcrack  stability  and  propagation 
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CHAPTER  VI 
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ABSTRACT 


The  stresses  due  to  "thermal  trapping"  of  Internally  absorbed  radiation 
in  solid  materials  are  analyzed  for  a  flat  plate  symmetrically  heated  by 
radiation  and  cooled  by  convection  with  finite  values  of  the  heat  transfer 
coefficient.  The  thermal  stress  state  consists  of  an  initial  transient 
followed  by  a  steady  state  stress  at  t  +  «  and  involves  a  reversal  in  sign 
of  the  stresses.  The  maximum  value  of  the  transient  stress  increases  with 
optical  thickness  whereas  the  maximum  value  of  the  steady-state  stress 
occurs  at  an  optical  thickness,  pa  -  1.3  with  zero  stresses  at  pa  =  0  and  °°. 


INTRODUCTION 


\ 

The  current  interest  in  solar  radiation  as  an  important  energy  source 
has  focused  attention  on  the  problem  of  thermal  stress  failure  of  materials 
of  construction  of  solar  installations.  Analyses  of  thermal  stresses  due 
to  radiative  heat  transfer  involve  non-linear  integral  equations.  For  this 
reason  the  usual  approach  of  obtaining  complete  transient  or  steady-state 
solutions  for  the  temperature  equations  from  which  the  thermal  stresses 
are  evaluated,  becomes  rather  complex.  To  circumvent  this  problem  at 
least  for  some  instances  it  should  be  noted  that  for  a  body  initially  at 
low  temperature  suddenly  subjected  to  radiation  heating,  the  temperature 
rise  at  the  time  when  the  thermal  stresses  are  at  or  near  their  maximum 
value,  is  sufficiently  small  that  any  re-emitted  radiation,  to  a  first 
approximation,  can  be  neglected.  For  practical  purposes  then,  the  body 
can  be  considered  to  be  subjected  to  a  constant  heat  flux. 

This  latter  approach  was  followed  in  an  analysis  of  thermal  stresses 
due  to  radiation  heating  of  materials  essentially  opaque  to  the  incident 
radiation,  which  was  absorbed  in  the  immediate  surface  regions  of  the 
material  [1].  This  analysis  was  extended  to  include  the  spectral  depen¬ 
dence  of  the  absorption  coefficient  typical  for  dielectric  materials, 
which  to  a  good  approximation  can  be  described  to  be  transparent  below 
a  given  wave-length  and  opaque  above  this  value  of  wave-length  [2], 

In  addition  to  complete  opacity  or  complete  transparency,  thermal 
radiation  can  also  interact  with  materials  such  that  it  is  absorbed  as 
it  is  being  transmitted  through  the  material.  The  rate  of  absorption  is 
governed  by  the  absorption  coefficient.  Internal  absorption  of  radiation, 
depending  on  the  external  cooling  condition,  can  lead  to  the  phenomenon 
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referred  to  as  "thermal  trapping"  in  which  the  interior  of  the  plate  reaches 
temperatures  higher  than  the  surface  or  surrounding  medium  13].  The  non- 
uniform  non-linear  temperature  distributions  associated  with  thermal  trap¬ 
ping  should  result  in  thermal  stresses.  This  effect  was  analyzed  recently 
for  a  flat  plate  symmetrically  heated  by  normally  incident  thermal  radiation 
and  cooled  by  convection  for  heat  transfer  coefficients  (h)  for  simplicity 
taken  as  h  =  0  and  °°  [4].  The  present  writers  became  aware  of  the  general 
solutions  of  Viskanta  et.  al  [3]  for  the  thermal  trapping  effect  involving 
finite  values  of  the  heat  transfer  coefficient.  Using  these  solutions, 
the  purpose  of  the  present  study  is  to  analyze  the  rhermal  stresses  in  a 
flat  plate  symmetrically  heated  by  radiation  and  cooled  by  convection  for 
finite  values  of  the  heat  transfer  coefficient.  Thermal  stresses  in 
assymmetrically  heated  plates  will  be  presented  in  future  reports. 


ANALYSIS 

A.  Boundary  and  Initial  Conditions 

The  flat  plate,  considered  infinite  in  extent,  is  located  in  the  yz- 
plane  and  -a  <  x  <  a.  At  x  =  a  and  -a  the  plate  is  subjected  at  t  =  0  to  a 
constant  normally  incident  uniform  thermal  radiation.  At  x  =  — a  and  a  the 
plate  is  cooled  by  Newtonian  convection  with  a  heat  transfer  coefficient  h. 
The  optical  properties  of  the  material  are  assumed  to  be  "grey"  i.e,, 
independent  of  wave-length.  The  reflectivity  of  the  plate  will  be  assumed 
sufficiently  low  that  the  effect  of  multiple  reflections  within  the  plate 
can  be  neglected.  The  material  properties  which  affect  the  temperature 
distribution  and  thermal  stresses,  such  as  the  eraissivity  (absorptivity), 
absorption  coefficient,  thermal  conductivity  and  thermal  diffusivity,  the 
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coefficient  of  thermal  expansion.  Young's  modulus  and  Poisson's  ratio  are 
assumed  to  he  independent  of  temperature.  As  stated  earlier,  it  will  be 
assumed  that  over  the  time  period  for  the  stresses  to  reach  a  value  equal 
to  or  nearly  equal  to  their  maximum,  the  temperatures  in  the  plate  remain 
sufficiently  low  that  the  intensity  of  the  re-emitted  radiation  to  a  first 
approximation  is  negligible  in  comparison  to  the  incident  radiation. 

The  intensity  of  the  heat  flux,  q  within  the  plate  due  to  the  thermal 
radiation  entering  at  x  =  a  is: 


q 


£q. 


-P(a-x) 

e 


(1) 


where  c  is  the  emissivity  (or  c  =  1  -  r  where  r  is  the  reflectivity)  q^ 
is  the  intensity  of  the  incident  heat  flux  and  y  is  the  absorption  coefficient. 
Similarly,  the  radiation  within  the  plate  due  to  the  incident  radiation 


at  x  =  -a  is: 


q  =  cqoe 


-y (a+x) 


(2) 


From 


eqs.  1  and  2  the  total  radiative  heat  flux  (q^)  within  the  plate  is: 


q .  =  2eq  e 
t  o 


■pa 


cosh  (P  x) 


(3) 


The  rate  of  internal  heat  absorption  (g"')  within  the  plate  is: 


g"'  =  2yEq^e  lJacosh(Px) 


(4) 


Solutions  for  the  transient  temperature  (T)  are  obtained  by  solving  [5]: 


32T  g"'(x)  =  1  3T 
3x’z  k  k  9x 


where  k  is  the  thermal  conductivity  and  <  is  the  thermal  diffusivity. 

For  the  present  heating  and  cooling  conditions  the  initial  and  boundary 


conditions  are: 


(6) 


T(x,0)  =  T 


I  «.*> 


3T  ( 

3^  (a’° 


-  £  (T(a,t) 


3T  .  . 

3x  (“a>C)  = 


■klT^«  -  V 


where  h  is  the  heat  transfer  coefficient. 

For  h  =  0,  the  initial  and  boundary  conditions  are  given  by  Eqs.  (6) 
and  (7)  and 


f  (-a,t)  -fl  (a,t)  =0 


B.  Solutions  for  transient  temperatures. 


Solution  of  Eq.  (5)  for  the  initial  and  boundary  conditions  given  by 


Eqs.  (6),  (7),  (8)  and  (9)  yields: 


T(x,t)  =  T  +  E  B  (t)  cos  (X  x) 
o  n=o  n  n 


where  A  are  the  roots  of  the  transcendental  equation 
n 


A  tan  (A  a)  =  h/k 
n  n 

G  2 

and  B  (t)  =  -~T~  (1  "  e  n  ) 

n 

2<pcq  e  y  cos  (X  a)  (ewa-e  ya)  +  X  sin(X  a)  (e'l3+e  pa) 

, . ,  or  n  n  n  , 

with  G  “  — -  [ - ] 

n  N  k  2  .  ,2 

n  y  +  X^ 

n 


and  N  =  a  +  cos  (^_a) 

n  Xzk  u 

n 
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C.  Thermal  stresses. 


Kxnressions  for  the  thermal  stresses  o  (x,t)  are  obtained  by  sub- 

y,z 

stitution  of  Eq.  (11)  for  T  in  [6]: 


u 

y.z 


aE 

1-v 


[-T 


a 


Tdx  + 

2a3 


a 

f  Txdx] 
-a 


(1A) 


which  yields: 


a 

y*z 


-aE 


z“  B  (t)  cos  (*nx)  +  - 

n=o  a (1-v) 


-  Bn<0 

Z  — - -  sin  (A  a) 

A  n 


For  convenience  the  numerical  results  are  reported  in  terms  of  the 
non-dimensional  temperature: 

T*  =  (T  -  T  )k/eq  a 
o  o 

the  non-dimensional  stress: 

■k 

a  =  o(l-v)k/aEqoea 


the  non-dimensional  time: 


and  the  non-dimensional  distance,  x/a. 


RESULTS  AND  DISCUSSION 

Figure  1  shows  the  stresses  at  the  center  (x=o)  and  surface  (x  -  -a  and  a) 
of  the  plate  (with  a  =  1  cm)  as  a  function  of  time  for  a  range  of  values 

of  the  optical  tickness,  and  h/k  arbitrarily  chosen  equal  to  0.09  cm 

Figure  2  shows  the  stresses  in  the  center  and  surface  of  the  plate 

(a  =  1  cm)  as  a  function  of  time  for  a  range  of  values  of  h/k  (cm  ^)  for 

arbitrary  value  of  the  optical  thickness,  pa  =  3. 
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Figures  3  and  A  give  the  maximum  values  of  the  transient  and  steady- 
state  thermal  stresses  as  a  function  of  optical  thickness,  |ia. 

From  the  above  figures  a  number  of  effects  may  be  noted.  Firstly, 
it  may  be  observed  that  after  an  initial  transient  the  thermal  stresses 
reach  a  steady-state  value  at  t  The  maximum  value  of  the  transient 

stresses  increases  with  increasing  optical  thickness,  pa  as  shown  in  Fig.  3. 
In  contrast,  the  steady-state  value  of  the  thermal  stresses  as  shown  in 
Fig.  A  exhibits  a  maximum  at  pa  =  1.3  and  approaches  zero  as  pa  +  0  or  °°. 
Secondly,  the  magnitude  of  the  transient  thermal  stresses  is  a  function 
of  h/k,  whereas  the  steady-state  stresses  are  independent  of  h/k.  It  is 
of  interest  to  note  that  the  maximum  transient  stresses  decrease  with  in¬ 
creasing  h/k,  in  contrast  to  many  literature  solutions  for  transient  thermal 
stresses  in  convective  heat  transfer  which  indicate  that  the  magnitude  of 
maximum  thermal  stress  increases  with  increasing  heat  transfer  coefficient. 
The  independence  of  the  steady-state  stress  on  h/k  occurs  because  the 
profile  of  the  temperature  distribution  within  the  plate  is  independent  of 
h/k,  in  spite  of  the  fact  that  the  magnitude  of  the  temperatures  for  a  given 
profile  is  an  inverse  function  of  the  heat  transfer  coefficient.  A  steady- 
state  temperature  distribution  will  be  attained  regardless  of  the  magnitude 
of  the  heat  transfer  coefficient,  with  the  exeception  of  h  =  0,  For  this 
latter  value  of  h,  although  the  thermal  stresses  after  an  initial  transient 
will  attain  a  constant  value  [A],  the  temperatures  will  continue  increasing 
linearly  with  time  (at  least  within  the  assumption  that  the  re-emitted 
radiation  can  be  neglected) . 

The  "thermal  trapping"  and  associated  thermal  stresses  are  most  ef¬ 
fectively  illustrated  for  the  limiting  value  of  the  heat  transfer  coef¬ 
ficient  h  =  ”.  For  this  value  no  increase  in  surface  temperature  occurs, 
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and  the  temperature  rise  within  the  plate  is  due  to  the  internal  absorption 
of  the  incident  radiation.  As  indicated  by  Fig.  4,  this  effect  is  most 
pronounced  at  an  optical  thickness  pa  =  1.3. 

Probably  the  most  interesting  observation  as  indicated  by  Figs.  1  and 
2  is  the  reversal  in  sign  of  the  stresses  in  going  from  the  transient  to 
the  steady-state  condition.  This  phenomenon  is  the  direct  result  of  the 
thermal  trapping  and  the  spatial  non-uniformity  of  the  heat  generation  within 
the  plate,  which  is  a  maximum  in  the  surface  regions  of  the  plate.  During 
the  very  initial  time  periods  after  initiation  of  radiation  at  t  =0,  the 
regions  of  the  plate  in  the  immediate  surface  area  will  be  at  a  higher  tem¬ 
perature  than  the  interior  of  the  plate.  During  this  period,  the  heat 
losses  by  convection  from  the  plate  are  relatively  small  for  small  values  of 
h  with  the  heat  flowing  primarily  from  the  surface  regions  towards  the 
interior  of  the  plate.  This  results  in  a  temperature  distribution  such 
that  the  surface  is  placed  in  compression,  whereas  the  center  of  the  plate 
is  in  tension.  However,  as  the  plate  approaches  the  steady-state  condition 
(t  -»  <*>)  ,  the  surface  temperatures  will  be  sufficiently  high  that  the  heat  losses 
at  the  surface  by  convection  become  appreciable.  In  fact,  the  maximum  surface 
temperature  will  correspond  to  the  condition  that  the  heat  lost  by  convection 
equals  the  total  radiation  energy  absorbed.  Under  this  condition,  all  heat 
generated  within  the  plate  must  flow  outward.  This  implies  that  the  maximum 
temperature  occurs  at  the  center  of  the  plate  (for  this  symmetric  case) ,  which 
will  result  in  a  tensile  stress  in  the  plate  surface  and  a  compressive  stress 
in  the  plate  interior.  However,  as  h/k  -*  <*>,  no  stress  reversal  in  the  sur¬ 
face  takes  place  because  the  surface  temperature  remains  constant  re¬ 
sulting  in  the  tensile  stress  for  t  >  0.  This  effect  of  the  reversal 


of  stress  is  illustrated  in  Figs.  5  and  6  which  give  the  spatial  variation 

ic 

of  the  temperature  and  stress  for  increasing  values  of  t  for  the  specific 
case  of  h/k  =  1  cm  ^  and  pa  =  10.  As  far  as  these  authors  are  aware  the 
occurrence  of  such  a  stress  reversal  has  not  been  reported  in  the  literature 
to  date.  It  appears  that  this  phenomenon  requires  a  spatially  non-uniform 
internal  heat  generation  and  may  not  only  be  encountered  in  the  absorption 
of  thermal  radiation  but  possibly  nuclear  radiation  as  well,  such  as  in 
the  form  of  neutrons  or  other  nuclear  particles. 
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Figure  1.  Maximum  thermal  stress  in  a  partially  absorbing  flat  plate 

symmetrically  heated  by  normally  incident  radiation  and  cooled 
by  convection  for  a  range  of  values  for  optical  thickness,  pa 
with  h/k  =  0.03  cm.-^  and  a  =  1  cm. 


Figure  2.  Time  dependence  of  maximum  thermal  stress  in  a  partially  absorbing 
flat  plate  symmetrically  heated  by  normally  incident  radiation 
and  cooled  by  convection  for  various  values  of  h/k  (cm."l)  with 
optical  thickness,  pa  =  3  and  a  =  1  cm. 

Figure  3.  Maximum  values  of  transient  thermal  stresses  as  a  function  of 

optical  thickness  pa  in  a  partially  absorbing  flat  plate  symmetri¬ 
cally  heated  by  normally  incident  radiation  and  cooled  by  convection 
for  various  values  of  h/k  (cm.-^-)  with  a  =  1  cm. 


Figure  4.  Maximum  values  of  steady  state  stress  as  a  function  of  optical 
thickness  pa  in  a  partially  absorbing  flat  plate  symmetrically 
heated  by  normally  incident  radiation  and  cooled  by  convection 
for  various  values  of  h/k  (cm.'l)  with  a  =  1  cm. 


Figure  5.  Temperature  distribution  in  a  partially  absorbing  flat  plate 

symmetrically  heated  by  normally  incident  radiation  and  cooled 
by  convection  for  various  values  of  time  with  h/k  =  1  cm.~l, 
pa  =  10  and  a  =  1  cm. 


Figure  6.  Stress  distribution  in  a  partially  absorbing  flat  plate  symmetrically 
heated  by  normally  incident  radiation  and  cooled  by  convection  for 
various  values  of  time  with  h/k  =  1  cm.“l,  pa  =  10  and  a  =  1  cm. 
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Fig.  A.  Maximum  values  of  steady  state  stress  as  a  function  of  optical  thickness  pa  in  a  partially  absorbing 
flat  plate  symmetrically  heated  by  normally  incident  radiation  and  cooled  by  convection  for  various 
values  of  h/k  (cm.~i)  with  a  =  1  cm. 
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Fig.  5.  Temperature  distribution  in  a  partially  absorbing  flat  plate  symmetrically  heated  by  normally  incident 
radiation  and  cooled  bv  convection  for  various  values  of  time  with  h/k  =  1  cm."^-,  pa  =  10  and  a  “  1  cm 
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ABSTRACT 


Thermal  stresses  in  a  partially  absorbing  brittle  ceramic  plate 
assymmetrically  heated  by  radiation  and  cooled  by  convection  on  the 
same  surface  with  finite  and  infinite  heat  transfer  coefficient  are 
analyzed. 

Comparison  of  the  results  from  this  study  with  those  obtained  in 
Part  I,  where  the  plate  is  subjected  to  radiation  heating  in  the  front 
surface  and  cooled  by  convection  at  the  rear  surface,  indicates  that 
the  magnitude  of  the  maximum  steady  state  tensile  thermal  stresses 
are  nearly  identical.  Significant  differences  are  found  in  the  transient 
thermal  stresses  and  the  temperature  distribution.  The  relative  tempera¬ 
ture  levels  in  the  present  case  are  found  to  be  significantly  lower  than 
those  obtained  in  Part  I.  Implications  of  these  results  to  the  design 
and  operation  of  engineering  structures  such  as  concentrated  solar 
receivers  are  discussed. 


r. 


INTRODUCTION 


A  number  of  analyses  on  the  role  of  the  properties  which  affect  the 

thermal  stress  resistance  of  semi-absorbing  ceramics  were  presented 
12  3 

previously  ’  ’  .  For  these  studies  the  geometry  consisted  of  a  flat 
plate  subjected  to  symmetric  radiation  heating  and  convective  cooling  as 
well  as  assymmetric  radiation  heating  on  one  side  of  the  plate  with  cool¬ 
ing  on  the  opposite  face.  For  these  two  heating  and  cooling  conditions, 
the  magnitude  of  the  maximum  tensile  thermal  stress  was  found  to  be  a 
complex  function  of  the  optical  thickness  (absorption  coefficient  tiroes 
plate  thickness).  At  steady-state  (t  =  °°) ,  for  non-zero  heat  transfer 
coefficient  the  thermal  stresses  were  equal  to  zero  for  values  of  the 
optical  thickness  equal  to  0  and  00 .  For  a  non-zero  heat  transfer  co¬ 
efficient  the  maximum  value  of  steady-state  tensile  thermal  stress  for 
the  symmetric  and  assymmetric  radiation  heating  occurred  at  values  of 
the  optical  thickness  of  1.334  and  2,  respectively. 

Comparison  of  the  numerical  results  shows  that  for  the  symmetrically 
heated  and  cooled  plate  the  magnitude  of  the  thermal  stresses  considerably 
exceed  those  for  assymmetric  heating  and  cooling.  The  reason  for  this  is 
that  in  the  assymmetric  case,  the  non-uniformity  of  the  thermal  expansions 
in  part  can  be  accommodated  by  bending.  In  view  of  the  symmetric  tem¬ 
perature  distributions,  this  mode  of  deformation  cannot  occur  for  symmetric 
heating  and  cooling. 

Comparison  of  the  temperature  distributions,  on  the  other  hand,  shows 
that  for  high  values  of  optical  thickness  (i.e.  opaque  materials)  the 
temperature  differences  encountered  at  steady-state  for  assymmetric  heating 
and  cooling  significantly  exceed  those  for  symmetric  heating  and  cooling. 
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This  occurs  because  in  the  symmetric  heating  and  cooling  case,  the  heat 
is  removed  from  the  same  surface  subjected  to  radiation,  with  little  or 
no  heat  absorption  in  the  interior  of  the  material.  In  fact,  for  an 
optical  thickness  which  approaches  infinity,  the  temperature  within 
the  plate  at  steady-state  (time  t  ->  00 )  is  uniform.  Under  this  condition 
the  rise  in  temperature  above  ambient  is  inversely  proportional  to  the 
convective  heat  transfer  coefficient  (with  exception  of  h  =  0).  In  con¬ 
trast,  for  the  plate  heated  by  radiation  at  the  front  face  and  cooled 
by  convection  at  the  back  face,  all  the  heat  must  be  conducted  through 
the  thickness  of  the  plate.  This  will  result  in  a  highly  non-unifortn 
temperature  distribution  even  at  an  optical  thickness  of  infinity  and 
steady-state  conditions  for  which  the  magnitude  of  thermal  stresses  equals 
zero.  For  high  radiation  fluxes  this  can  result  in  excessive  front  face 
temperatures  even  for  moderate  or  low  rear  face  temperature  especially 
for  ceramic  materials  with  moderate  or  low  values  of  thermal  conductivity. 
This  latter  effect  can  have  important  implications  for  the  design  and 
operation  of  concentrated  solar  energy  receivers  which  generally  will 
operate  in  an  assymmetric  heating  mode  and  for  which  ceramic  materials 
are  serious  contenders  for  the  materials  of  construction.  In  order  to 
avoid  excessively  high  temperatures  at  the  front  face  subjected  to  radia¬ 
tion,  may  well  require  limiting  the  temperature  of  the  back-face.  This 
in  turn  would  reduce  the  temperature  of  the  heat  transfer  fluid  with  an 
associated  decrease  in  overall  efficiency  of  the  receiver.  If  so,  a  de¬ 
sign  modification  appears  in  order.  One  such  modification  is  suggested 
by  the  previous  analyses,  which  indicate  that  assymmetric  heating  and 
convective  cooling  at  the  same  side  of  the  plate  can  result  in  the  lower 
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thermal  stresses  associated  with  the  assymmetric  heating  and  the  greater 
temperature  uniformity  observed  for  the  symmetric  heating  and  cooling. 

Accordingly,  the  objective  of  this  study  is  to  analyze  the  role  of 
the  material  properties  in  the  thermal  stress  resistance  and  temperature 
distribution  in  a  ceramic  plate  subjected  to  radiative  heating  and  con¬ 
vective  cooling  at  the  same  side,  with  no  heat  transfer  (i.e.,  completely 
insulated)  at  the  opposite  (back)  face  of  the  plate. 

II.  ANALYSIS 

The  flat  plate  of  thickness  2a,  located  parallel  to  y-z  plane  and 
bounded  by  x  =  a  and  -a,  is  initially  at  a  uniform  temperature,  T  .  It 
is  heated  by  an  incident  radiation,  qQ,  and  cooled  by  Newtonian  convection 
with  heat  transfer  coefficient,  h,  at  x  =  -a.  Similar  to  the  previous  study^ 
it  was  assumed  that  the  reflectivity  of  the  plate  is  small  enough  to 
neglect  the  effect  of  multiple  internal  reflection  and  the  physical 
properties  of  the  material  are  independent  of  temperature.  It  was  also 
assumed  that  the  temperature  of  the  plate,  over  the  time  period  required 
for  stresses  to  reach  their  maximum,  was  low  enough  that  the  re-emission 
of  the  absorbed  energy  is  negligible  in  comparison  to  the  incident  rad¬ 
iation.  In  addition,  the  optical  properties  of  the  plate  was  taken  to 
be  independent  of  the  wavelength  of  the  incident  radiation. 

A.  Temperature  Distribution 

Due  to  the  normally  incident  radiation,  q^,  on  the  front  face  (x  =  -a) 
of  the  plate,  the  flux  transmitted  through  the  surface  is  EqQ  where  e=l-r 
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with  r  being  the  reflectivity.  The  rate  of  internal  heat  generation 
(g’"  )  within  the  plate  is  given  by': 


(x)  =  peqoe 


-p (a+x) 


where  11  is  the  absorption  coefficient.  The  temperature  solutions  (T) 
for  the  specific  initial  and  boundary  conditions  were  obtained  by  solving 


the  differential  equation  : 


32T  g'”  (x)  =  f  1.  3T 
3x2  k  V  3x 


where  k  is  the  thermal  conductivity  and  <  is  the  thermal  diffusivity. 

Equation  (1)  was  solved  for  three  values  of  heat  transfer  coefficient, 
h,  viz.  h  =  0,  finite  and  infinity.  The  initial  and  boundary  conditions 
along  with  the  temperature  solutions  are  as  follows: 


a.  Finite  heat  transfer  coefficient,  h 


The  initial  and  boundary  conditions  are: 

-s  T 

T(x,0)  =  Tq;  (a ,  t)  =  0  (2) 

k  9--(-a,t)  =  h(T(-a,t)  -  Tq}  (3) 


The  solution  of  equation  (1)  subjected  to  the  initial  and  boundary  con¬ 
ditions  (2)  and  (3)  is  given  by: 


where 


T(x,t)  =  T  +  l  B  (t)  Cos  X  (a-x) 
°  n  n  n 

n=0 


c  >2 
B  (t)  =  — n- —  [l  -  e~'cAn  t] 

n  ^n2 


KqQpc  X^Sin(2X  a)  -  ue 

Gn  =  mr  ^  p2  +  X  2 

n  n 


-2pa 


(4) 

(5a) 


(5b) 
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a  + 


(5c) 


Nn  =  3  +  IFk  cos2(2Xna> 
n 


and  X  is  the  root  of  the  transcendental  equation: 
n 


X  tan  ( 2 X  a)  =  — 
n  n  k 


(6) 


b.  Heat  transfer  coefficient,  h  =  00 


The  initial  and  boundary  conditions  are: 


T(x,0)  =  T(-a,t)  =  T 


E  (a’t}  =  0 


(7) 

(8) 


The  temperature  solution  for  conditions  (7)  and  (8)  is  given  by: 


where 


and 


T(x,t)  =  T  +  E°°  B  (t)  sin{X  (x+a)}  (9) 

o  „  n  n 

n=0 

G  l2. 

Bn(t)  =  t1  "  e"  n  1  (10a) 

xq  ye  X  -  ye  2,Jasin(2Xna) 

G„  *  Hsr  - ]  <10b) 

n 

X  =  ( 2.^~1-)-7L  with  n=Q,  1,  2 -  (10c) 

n  4a 


The  thermal  environment  with  h=0  is  identical  to  the  one  discussed 
in  an  earlier  study  and,  therefore,  the  expressions  for  temperature  and 
thermal  stresses  for  this  case  can  be  found  in  reference  3. 


B.  Thermal  Stresses 


The  expressions  for  thermal  stresses  were  obtained  from 


a  =  f^-[-T  +  —  /3 Tdx  +  fa Txdx] 
y,z  1-v  2a  -a  2a3  -a 


(11) 
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where  a  is  the  coefficient  of  thermal  expansion,  E  is  the  Young's  modulus 
of  elasticity,  v  is  the  Poisson's  ratio  and  T  is  the  temperature, 

fhe  stresses  were  obtained  for  two  cases  with  finite  and  infinite 
heat  transfer  coefficients: 

a.  Finite  heat  transfer  coefficient,  h 

For  tliis  case  the  temperature  distribution  T  in  equation  (4)  was 
substituted  in  equation  (11)  to  yield  the  thermal  stress  distribution: 


°y.z  =  ~  T^v  Bn(t)cos(An(a-x)} 
n=0 


t  aE  r-  B„<0  Sln(2»na) 
„-0  »„ 


~  „  l-cos(2A  a) 

+  3»E  S-B  (t)[_ 


2a3(1-v)^_0  n 


A  2 
n 


a  sin(2A  a) 
- 

n 

(12) 


where  Bn(t)  and  A^  are  defined  in  equations  (5)  and  (6) , 


b.  Heat  transfer  coefficient  h  = 

Thermal  stress  distribution  for  this  case  was  obtained  from  equations 
(9)  and  (11)  : 


3  =  -7— - r  £  B  (t)  sin  f  A  (x+a)  } 

y,z  (1-v)  n  n 


aE 


2a(l-v) 


r°°  B  ( t ) 

l  n 

n=l  A 


3xaE 


2a3(l-v)  n 

n=l 


sin(2A  a) 
B  (t)[- 


V 


(13) 
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where  B  (t)  and  A  are  defined  in  equation  (10). 
n  n 

The  maximum  steady  state  tensile  stress  is  found  to  occur  in  the 
front  face  of  the  plate  (x  =  -a)  and  is  given  by: 


•'  (-a,00) 

max 


ctF. 

2a( 1-v) 


t 

n=l 


3a  E _ °°  ^n 

2a2(l-v)  xA  2 
n=l  n 


[■ 


sin(2A 

n 


a) 


(14) 


where  is  defined  in  equation  (10b). 

Equation  (14)  can  be  simplified  for  optically  very  thick  (ua  -+  “) 
and  optically  very  thin  (ua  <<  1)  plates  to: 


=  0  (pa  -*■  »)  (15n) 

TTla  X 

aEq  cpa2 

°max(~a  ,a>)  =  3(l-vjk  (va<<1)  (15b) 

Equating  a  with  the  tensile  strength  S  ,  the  expressions  for  the 
max  c 

maximum  permissible  heat  flux  q  ^  to  which  the  plate  can  be  subjected 
without  the  risk  of  failure  under  steady  state  condition  is  given  by: 


q  -*■  00  (pa  -*■  “>)  (16a) 

max 


3S  (l-v)k 

q  = — £ - - -  (pa  <<1) 

wax  aEcya2 


(16b) 


III.  NUMERICAL  RESULTS  AND  DISCUSSION 

1-4 

Following  previous  practice  ,  the  temperatures,  stresses  and  time 
will  be  presented  in  the  non-dimensional  form: 
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Cl  7a) 


T*  "  (T  '  To)k/cV 
* 

a  =  ci(l-y)  k/aEeq^a  (17b) 

and  a  =  <t/a?  (17c) 

for  convenience,  hereafter,  referred  to  as  temperature,  stress  and  time, 
respectively. 

Figures  la  and  2a,  for  a  range  of  values  of  time  and  h  =  0.06  watts, 
cm.  2°C  1.  show  the  spatial  variation  of  the  temperature  within  the  plate 
for  values  of  optical  thickness  pa  =  0.5  and  30,  respectively.  These 
results  indicate  that  during  the  initial  stages  of  heating,  the  regions 
of  the  plate  near  the  front  surface  are  at  a  higher  temperature  than  those 
regions  near  the  back  surface.  This  is  because  the  rate  of  heat  generation 
within  the  plate  near  the  front  is  higher  than  near  the  back.  This  effect 
is  more  pronounced  for  the  higher  value  of  optical  thickness,  since  the 
incident  radiation  does  not  penetrate  as  far  into  the  plate  as  for  the  lower 
values  of  optical  thickness. 

At  higher  values  of  time  (t  -*■«>),  however,  the  temperature  near  the 
rear  of  the  plate  exceeds  those  near  the  front.  This  effect,  is  expected 
for  any  non-zero  value  of  the  heat  transfer  coefficient  since  all  the  ab¬ 
sorbed  heat  must  be  removed  in  the  front,  resulting  in  a  heat  flow  in  a 
direction  opposite  to  that  of  the  incoming  radiation. 

Figures  lb  and  2b  show  the  distribution  of  stresses  corresponding 
to  the  temperatures  shown  in  figs,  la  and  2a.  For  pa  =  0.5,  the  stresses 
in  the  front  surface  of  the  plate  are  always  tensile  even  at  the  shortest 
values  of  time  at  which  the  front  of  the  plate  is  at  higher  temperature. 

As  pointed  out  in  Part  I  of  this  study,  tensile  thermal  stresses  can  occur 
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in  the  hottest  part  of  a  structure,  when  the  temperature  distribution 
is  assymmetric  and  exhibits  a  downward  concavity.  This  effect  is  absent 
for  pa  =  30  for  which  the  temperature  distributions  exhibit  an  upward 
concavity  over  most  of  its  range  with  the  exception  of  the  region  immediate¬ 
ly  at  the  surface  of  the  plate  at  the  lower  values  of  time. 

For  purposes  of  comparison  with  the  results  obtained  in  Part  I  of 
this  study,  fig.  3  shows  the  spatial  distribution  of  temperature  at  t  =  » 
and  h  =  0.06  watts. cm.  ^°C.  ^  for  a  range  of  values  of  pa.  For  other 
values  of  the  heat  transfer  coefficient,  the  temperature  distributions 
are  identical  to  those  shown  in  fig.  3  with  the  exception  that  the  tempera¬ 
ture  at  the  front  surface  of  the  plate  is  inversely  proportional  to  the 
heat  transfer  coefficient.  As  expected,  for  all  values  of  pa,  the  tem¬ 
perature  distribution  shows  a  positive  gradient  since  all  the  absorbed 
heat  is  removed  at  the  front  of  the  plate.  For  pa  =  the  plate  is  at 
uniform  temperature  throughout,  since  at  thermal  equilibrium  all  the  heat 
absorbed  by  radiation  in  the  immediate  surface  is  removed  by  convention 
without  flowing  into  the  plate  interior. 

Figure  4  shows  the  effect  of  the  magnitude  of  the  heat  transfer 
coefficient  on  the  magnitude  of  maximum  tensile  stress  for  pa  =  3.0.  The 
magnitude  of  the  peak  transient  tensile  thermal  stress  at  the  beginning 
of  the  thermal  history  is  an  inverse  function  of  the  heat  transfer  co¬ 
efficient.  At  thermal  equilibrium  (t  ,  the  magnitude  of  stress  is 
independent  of  the  value  of  the  heat  transfer  coefficient. 

Figures  5a  and  5b  show  the  transient  behavior  of  the  maximum  tensile 

^20 

stress  for  a  range  of  values  of  pa  and  values  of  h  =  0.06  watts. Cm.  C, 
and  °°,  respectively.  For  h  =  0.06  watts. Cm.  ^°C.  ^  and  pa  2  the  stresses 
go  through  a  transient  peak,  before  reaching  their  equilibrium  value. 


For  h  =  0.06  watts.cm.  ^°C.  ^  and  pa  <  2  and  all  values  of  pa  for  h  =  00 
no  such  peak  transient  thermal  stress  exists. 

Figure  6  summarizes  the  values  for  maximum  tensile  stress  as  a  function  of 
optical  thickness,  pa.  Above  an  optical  thickness,  pa  i  11,  the  peak  transient 
stresses  exceed  the  steady-state  stresses.  Conversely,  for  pa  <  11,  the  steady- 
state  stresses  exceed  the  transient  stresses.  For  pa  <  2  no  peak  transient 
stresses  are  found.  Maximum  steady-state  stresses  are  found  at  pa  ;  2. 

A  comparison  of  the  present  results  with  those  obtained  in  Part  I 
of  this  study  indicates  that  for  assymmetric  radiation  heating  with  con¬ 
vective  cooling  at  the  front  or  rear  surface,  the  magnitude  of  the  maximum 
steady-state  tensile  thermal  stresses  are  nearly  identical, but  significant 
differences  arc  found  in  the  transient  thermal  stresses  and  the  temperature 
distributions. 

For  purposes  of  discussion  of  this  latter  effect,  the  value  of  optical 
thickness,  pa  =  can  be  taken  as  an  extreme  case.  At  steady-state  (t  =  °°) 
the  stresses  for  both  front  and  rear  cooling  are  identically  equal  to  zero. 

For  front  cooling  as  found  in  the  present  study,  the  zero  thermal  stress 
state  results  from  the  temperature  state  which  is  completely  uniform  through¬ 
out  the  plate  thickness  as  shown  in  fig.  3.  In  contrast,  for  cooling  at  the 
rear  surface,  the  heat  absorbed  in  the  front  is  transferred  by  solid  con¬ 
duction  throughout  the  thickness  of  the  plate.  As  shown  in  fig.  4a  of  Part 
I  of  this  study,  this  results  in  a  linear  temperature  distribution.  Such 
a  linear  temperature  distribution  in  an  unconstrained  plate,  also  results 
in  a  zero-thermal  stress  state.  The  temperature  gradient  in  this  case 
equals,  dT/dx  =  -cq  /k. 

The  above  observation  is  critical  to  the  design  and  operation  of 
structures  assymmetrically  heated  by  radiation  and  cooled  by  convection 
such  as  receivers  of  concentrated  solar  energy.  For  high  collector 
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efficiency  it  is  critical  that  the  temperature  of  the  working  fluid  be 
as  high  as  possible.  This  requires  the  high  temperature  of  the  surface 
from  which  the  heat  is  transferred  to  the  working  fluid.  For  cooling  at 
the  rear  surface,  this  implies  that  the  front  surface  will  be  at  a  higher 
temperature  than  the  rear  surface,  as  for  pa  »  all  absorbed  radiation 
must  be  transported  through  the  plate  by  solid  conduction. 

This  effect  can  be  illustrated  by  a  numerical  example.  For  a  peak 
3  2 

solar  intensity  of  10  watts/m  concentrated  by  a  factor  of  200  at  the 
surface  of  the  plate,  with  a  value  of  thermal  conductivity  of  0.05  watts/cm. °C 
(representative  of  a  ceramic  at  higher  temperature)  and  a  value  of  e  =  1-r  = 
0.9,  the  temperature  gradient  in  the  plate  dT/dx  -  -300°C.cm  For  a 
plate  of  unit  thickness,  this  value  of  temperature  gradient,  yields  a  value 
for  the  temperature  difference  between  the  front  and  rear  of  the  plate  of 
300°C.  For  rear  surface  temperatures  as  high  as  possible  such  a  tempera¬ 
ture  difference,  will  result  in  front  surface  temperatures,  which  could 
result  in  oxidation,  creep,  evaporation  and  greatly  reduced  service-life, 
depending  on  the  ceramic  chosen  as  structural  material.  Furthermore,  even 
in  the  absence  of  thermal  stress,  high  temperature  gradients  within  a 
material  can  cause  excessive  warpage  which  also  must  be  taken  into  account. 

For  a  temperature  gradient  of  300°C.cm  ^ ,  the  radius  of  curvature  R,  in¬ 
duced  in  a  plate  of  aluminum  oxide  (a  -  8  x  10  t  becomes  R  =  400  cm. 

Not  only  will  this  cause  large  undesirable  displacement,  but  the  plate 
may  even  start  acting  as  a  lens. 

No  such  disadvantages  are  found  for  a  plate  heated  by  radiation  and 
cooled  by  convection  on  the  same  side.  For  pa  =  00  and  t  =  »,  the  tempera¬ 
tures  are  uniforn  and  of  a  magnitude  considerably  below  those  for  a  plate 
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heated  by  radiation  and  cooled  by  convection  at  opposite  surfaces. 

As  indicated  by  the  combined  results  of  Part  I  and  II  of  this 
study,  from  the  point  of  view  of  minimizing  temperatures  and  thermal  stress¬ 
es,  a  solar  receiver  consisting  of  a  dual-plate  design  can  be  used  to  ad¬ 
vantage.  The  front  plate  should  be  highly  transparent  (pa  -*■  0)  ,  to  trans¬ 
mit  as  much  solar  radiation  as  possible.  On  the  other  hand  the  second  plate 
should  be  highly  opaque  (pa  °°)  .  The  working  fluid  should  be  placed 
between  the  plates  to  be  heated  by  the  second  plate  which  absorbs  the  radia¬ 
tion.  This  design  reduces  problems  with  temperature  related  material 
degradation  to  a  minimum.  The  above  is  a  further  illustration  of  the  well 
known  concept  that  structures  made  of  structural  ceramics  must  be  designed 
in  such  a  manner  to  fully  benefit  from  the  advantages  ceramic  materials  can 
offer. 
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Fig.  lb.  Spatial  variation  of  thermal  stresses  for  various  values  of  time  in  a  partially  absorbing  flat 

plate  subjected  to  radiation  heating  and  convection  cooling  at  the  front  face  with  heat  transfer 
coefficient,  h  =  0.06  (watts.  cm.-2oc-l-^  K  =  9.3  watts>  Cm.‘ll,C-l,  ua  =  0.5  and  a  =  1  cm. 
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Spatial  variation  of  thermal  stresses  for  various  values  of  time  in  a  partially  absorbing  flat 
plate  subjected  to  radiation  heating  and  convection  cooling  at  the  front  surface  with  heat  rra 
fer  coefficient,  h  =  0.06  (watts,  cm.-^oc--*-) ,  K  =  0.3  watts.  cm.-l°C-l,  ua  =  30  and  a  =  1  cm. 
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^8*  Spatial  distribution  of  steady  state  (t  -*■  °°)  temperature  for  various 
vaiucs  of  pa  in  a  partially  absorbing  flat  plate  subjected  to  radia¬ 
tion  heating  and  convection  cooling  at  the  front  face  with  heat 
transfer  coefficient,  h  =  0.06  (wat ts . cm. -2*c~l) ,  k  =  0.3  watts. cm.' 
°C-1  and  a  =  1 
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CHAPTER  VIII 


EFFECT  OF  NATURE  OF  CONCAVITY  OF  TEMPERATURE  DISTRIBUTION 
ON  POSITION  AND  SIGN  OF  MAXIMUM  THERMAL  STRESS1 2 3 

2  3  2 
K.  Satyamurthy  ,  D.  P.  H.  Hasselman  and  J.  P.  Singh 


INTRODUCTION 


Engineering  solutions  for  the  thermal  stresses  in  structures  with 
positive  coefficient  of  thermal  expansion,  generally  indicate  that  the 
stresses  are  compressive  (negative)  at  or  near  the  position  in  the  struc¬ 
ture  where  the  temperature  has  its  highest  value.  Similarly,  the  tensile 
stresses  generally  are  found  at  or  near  the  position  of  lower  temperature. 

An  exception  to  the  above  observation  was  found  in  a  recent  analysis  [l] 
of  the  thermal  stresses  in  a  flat  plate  assymmetrically  heated  by  thermal 
radiation  absorbed  in  the  plate  interior.  The  results  of  this  study  showed 
that  the  stresses  exhibited  their  maximum  tensile  value  at  the  surface  of 
the  plate  being  irradiated  which  is  at  the  highest  value  of  temperature. 

In  order  to  provide  an  explanation  for  this  observation  it  should  be 
noted  that  thermo-elastic  theory  does  not  predict  a  unique  temperature  pro¬ 
file  to  achieve  a  thermal  stress  distribution.  The  validity  of  this  statement 

1  The  results  of  study  were  obtained  as  part  of  a  research  program  on  the 
thermo-mechanical  behavior  of  brittle  structural  materials  supported  by 
the  Office  of  Naval  Research  under  contract:  N00014-78-C-0431. 

2 

Research  Associate,  Department  of  Materials  Engineering,  Virginia  Polytechnic 
Institute  and  State  University,  Blacksburg,  Virginia  24061. 

3  Whittemore  Professor  of  Materials  Engineering,  Virginia  Polytechnic  Institute 
and  State  University,  Blacksburg,  Virginia  24061. 


can  be  examined  by  superposing  a  zero-stress  linear  temperature  distribution 
with  arbitrary  gradient(s)  [2]  on  any  other  temperature  distribution  regard¬ 
less  of  its  nature  or  source.  This  shows  that  any  arbitrary  relationship 
between  the  stress  and  temperature  at  any  point  in  a  structure  can  be  ob¬ 
tained  by  the  appropriate  choice  of  the  temperature  gradients. 

Such  a  superposition  of  linear  temperature  distribution  was  not  done 
in  the  plate  subjected  to  radiation  heating  referred  to  above.  For  this 
reason,  the  existence  of  the  maximum  tensile  thermal  stress  at  the  position 
of  highest  temperature  must  be  attributed  to  other  reasons.  An  examination 
of  the  general  results  showed  that  the  temperature  distribution  exhibited  a 
downward  cancavity  due  to  the  internally  absorbed  radiation.  Such  tempera¬ 
ture  distribution  contrasts  with  the  generally  concave  upward  temperature 
distribution  encountered  in  transient  heat  transfer  for  which  the  exchange 
of  heat  between  the  structure  and  environment  occurs  at  the  surface.  It 
was  postulated  that  the  nature  of  the  concavity  of  the  temperature  dis¬ 
tribution  as  the  result  of  the  internal  absorption  of  the  incident  radia¬ 
tion,  was  directly  responsible  for  the  finding  that  the  tensile  thermal 
stresses  occurred  in  the  hottest  part  of  the  plate.  The  purpose  of  this 
communication  is  to  present  a  simple  analysis,  which  demonstrated  the 
validity  of  this  hypothesis. 


ANALYSIS 


A  thin  plate  located  in  the  y-z  plane  and  -a<x<a  in  a  rectangular 
coordinate  system  will  be  considered.  The  plate  is  subjected  to  a  mono¬ 
tonic  distribution  of  temperature  T  of  the  form: 

T  =  T  [1+C(f^)n]  (1) 

o  za 


2 
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where  is  a  reference  temperature  and  C  and  n  are  constants.  For  C  >  0, 
the  temperature  distribution  Is  concave  downward  or  upward  for  0  <  n  <■  1  and 
n  >  1,  respectively.  Similarly,  for  C  <0,  the  temperature  distribution  is 
concave  downward  or  upward  for  n  >  1  and  0  <  n  <  1,  respectively. 

The  thermal  stress  distribution  can  be  calculated  from  [2] 


F  -i  3-  o  & 

a  =  ~~  [-T  +  ~  f  Tdx  +  /  Txdx]  (2) 

y,z  1-v  2a  „  3 
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Substitution  of  eq.  1  into  2  and  integration  yields: 
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Examination  of  eqs.  1  and  3  shows  that  for  a  concave  upward  temperature 
distribution  (C  >  0,  n  >  1  or  C  <  0,  0  <  n  <  1) ,  the  position  of  maximum  ten¬ 
sile  stress  coincides  with  the  position  of  minimum  temperature.  In  contrast 
for  temperature  distribution  with  downward  concavity  (C  <  0,  n  >  1  or  C  >  0, 

0  <  n  *■  1)  ,  the  maximum  tensile  stress  occurs  at  the  position  of  maximum 
temperature.  For  n  =  1,  the  temperature  distribution  is  linear  with 
zero  thermal  stress  as  expected  [2],  regardless  of  the  value  of  C.  These 
results  indicate  that  the  direction  of  the  deviation  from  linearity  (i.e. 
the  nature  of  the  concavity)  of  the  temperature  distribution  controls  the 
relative  position  of  maximum  thermal  stress  and  maximum  temperature,  in 
agreement  with  the  hypothesis  presented  earlier. 


DISCUSSION 


It  should  be  noted  that  if  a  symmetric  temperature  distribution  had 
been  chosen  the  position  of  maximum  tensile  thermal  stress  would  have 


occurred  at  the  position  of  lowest  temperature,  regardless  of  the  nature 


of  the  concavity  of  the  temperature  distribution.  It  appears  then  that  for 
the  position  of  the  maximum  tensile  stress  to  coincide  with  the  position  of 
the  highest  temperature,  requires  an  assymmetric  concave  downward  tempera- 
ture  distribution.  Such  a  temperature  distribution  will  result  from  an 
assymmetrically  distributed  non-uniform  internal  heat  generation.  This  is 
the  case  for  the  aforementioned  analysis  Cl]  of  the  thermal  stresses  in  a 
plate  subjected  to  assymmetric  thermal  radiation  absorbed  within  the  interior 
of  the  plate.  Similar  effects  can  be  expected  for  assymmetrically  absorbed 
nuclear  radiation  or  non-uni formly  distributed  radio-activity. 
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ABSTRACT 

A  fracture-mechanical  analysis  is  presented  for  the  effect  of 
crack  interaction  on  the  stability  and  propagation  of  cracks  in 
brittle  ceramic  materials  subjected  to  severe  thermal  stress, 
based  on  a  model  consisting  of  a  flat  plate  with  a  rectangular  ar¬ 
ray  of  columnar  and  co-linear  cracks.  It  is  shown  that  the  com¬ 
bined  effect  of  the  interaction  between  the  columnar  and  co-linear 
cracks  on  thermal  stress  intensity  factors  and  effective  elastic 
behavior  of  the  plate  results  in  a  major  improvement  in  thermal 
stress  fracture  initiation.  However,  extensive  unstable  crack 
propagation  and  interaction  of  closely  spaced  co-linear  cracks  will 


result  in  complete  failure  due  to  crack  coalescence. 


I.  INTRODUCTION 


A.  General 

Fngineering  structures  or  components  for  use  at  high  temperature 
must  rely  on  ceramics  as  the  only  suitable  class  of  materials  of 
construction.  Ceramic  materials  due  to  their  strong  ionic  or 
covalent  bonding  exhibit  high  melting  points,  high  chemical  sta¬ 
bility,  good  creep  resistance  and  other  favorable  properties. 

High-temperature  structures  or  components  usually  are  re¬ 
quired  to  operate  under  conditions  of  high  thermal  stress  during 
transient  heating/cooling  or  steady-state  operation.  Practical 
experience  has  shown  that  ceramic  materials  under  these  condi¬ 
tions  are  prone  to  catastrophic  failure.  One  of  the  principal 
reasons  for  this  is  that  in  these  materials,  in  contrast  to  many 
metals  and  allovs,  thermal  stress  relaxation  by  plastic  flow  can¬ 
not  readily  occur.  Furthermore,  the  corresponding  low  'values  of 
fracture  toughness  make  ceramic  materials  prone  to  failure  even 
at  moderate  thermal  stress  levels.  This  low  fracture  toughness 
also  tends  to  make  such  failure  completely  catastrophic  with  the 
crack  traversing  through  the  component,  rendering  it  unsuitable  for 
continued  satisfactory  service.  This  latter  pehnomenon  is  criti¬ 
cal  since  frequently  in  practice  thermal  environments  can  be  of 
such  severity  that  failure  even  in  ceramic  materials  with  opti¬ 
mum  thermal  stress  resistance  cannot  be  avoided.  In  fact,  mate- 


riaL  development  in  refractory  technology  usually  is  more  directed 
toward  promoting  rapid  crack  arrest  in  improving  thermal  stress 
fracture,,  rather  than  avoiding  such  fracture  altogether. 

For  this  reason,  it  is  imperative  that  for  the  purpose  of 
reliable  engineering  design  and  selection  of  materials,  the  factors 
and  mechanisms  which  affect  thermal  stress  failure  are  well  under¬ 
stood.  Such  understanding  should  cover  the  initiation  of  frac¬ 
ture  as  well  as  tlie  propagation  of  cracks  following  thermal  stress 
fracture.  Since  thermal  stress  fracture  frequently  involves 
multiple  crack  formation,  the  effects  of  crack  interaction  should 
be  understood  as  well.  This  latter  topic  represents  the  objective 
of  the  present  study.  In  order  to  provide  a  proper  framework  for 
the  analysis  and  results  to  be  presented,  a  brief  review  of  the 
major  variables  which  control  the  initiation  of  thermal  stress 
fracture  and  crack  propagation  in  brittle  ceramics  will  be  given 
f i rst . 


B.  Background 

Thermal  stress  fracture  of  ceramic  materials  can  arise  from 
a  wide  variety  of  thermal  and  mechanical  conditions  including  ex¬ 
ternal  constraints,  steady-state  or  transient  heat  flow,  convec¬ 
tive  or  radiative  heat  transfer,  thermal  buckling  and  many  others 
For  the  purpose  of  material  selection  and  the  development  of  re¬ 
liable  design  criteria,  the  role  of  the  pertinent  external  vari- 


nblos  and  mart-rial  properties  which  affect  thermal  stress  failure 
should  he  carefully  assessed.  For  this  reason  it  is  conventional 
practice  in  ceramic  technology  to  obtain  from  the  literature  or 
derive  expressions  for  the  thermal  stresses  for  a  given  mode  of 
heat  transfer  and  mechanical  boundary  conditions.  By  defining 
a  failure  condition,  these  expressions  can  be  rearranged  to  yield 
the  maximum  thermal  environment  to  which  a  material  can  be  sub¬ 
jected  without  incurring  the  risk  of  thermal  stress  failure  in 
terms  of  the  appropriate  environmental,  geometric  and  material 
variables. 

As  an  example,  the  maximum  allowable  temperature  difference 
across  the  thickness  of  a  concentric  hollow  cylinder  undergoing 
steady-state  heat  flow  can  be  written  (Kingery,  1955;  Hasselman, 
1978)  : 

AT  =  Co ,( 1 -v) /aE  (1) 

max  r 

where  C  is  a  geometric  constant  related  to  the  ratio  of  outer  to 
inner  radii  of  the  cylinder,  o  is  the  tensile  fracture  stress, 
v  is  Poisson’s  ratio,  a  is  the  coefficient  of  thermal  expansion 
and  E  is  Young's  modulus  of  elasticity. 

As  a  further  example,  the  maximum  sudden  change  in  ambient 
temperature  (AT  )  to  which  a  solid  cylinder  under  conditions 
of  convective  heat  transfer  can  be  subjected  (.Taeger,  1945)  for 
Bolt  number,  8  -  bh/k  <<  1  is: 
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AT  -  4.95  .Tc(l-u)k/bhaE  (2) 

max  r 

where  k  Is  the  thermal  conductivity,  h  is  the  heat  transfer 

coefficient  and  b  is  radius  of  cylinder. 

Also,  the  maximum  black-body  radiation  temperature  (T  ) 

max 

to  which  a  spherical  body  at  low  initial  temperature  can  be  sub¬ 
jected  is  (Hasselman,  1963): 

1/4 

T  =  (5  a,(l-v)k/aEpeb)  '  (3) 

max  t 

where  r  is  the  absorptivity  and  p  is  Stef an-Boltzmann  constant. 

Eqs.  1,  2  and  3  suggest  that  for  a  given  geometry  and  body 
size,  ceramic  materials  with  high  thermal  stress  resistance  for 
the  three  heat  transfer  conditions  referred  to  above  should  have 
high  values  of  the  "thermal  stress  resistance  parameters": 

o^.(l-v)/aE  ;  o^(l-v)k/aF.  ;  {a  ^  ( l-v)k/aEE  } 1  ^ 

(4) 

These  three  parameters  suggest  that  ceramic  materials  with 
high  thermal  stress  resistance  should  have  high  values  of  tensile 
strength  and  thermal  conductivity  in  combination  with  low  values 
of  the  coefficient  of  thermal  expansion.  Young's  modulus  of  elas¬ 
ticity,  absorptivity  and  Poisson's  ratio.  A  review  of  the  above 
and  many  other  thermal  stress  resistance  parameters  for  other 
thermal  environments  and  failure  phenomena  was  presented  recently 
by  one  of  the  present  writers  (Hasselman,  1978).  Important  to 
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note  is  that  these  thermal  stress  resistance  parameters  are  ap¬ 
propriate  to  the  initiation  of  thermal  stress  fracture. 

The  nature  of  crack  propagation  in  thermally  stressed  brittle 
ceramics  was  analyzed  by  Hasselman  (1969,  1970).  A  simple  mechani¬ 
cal  model  in  the  form  of  a  flat  plate  was  considered.  The  plate 
was  uniaxially  constrained  from  free  thermal  expansion.  The  plate 
had  N  cracks  per  unit  area  of  equal  size,  oriented  perpendicularly 
to  the  direction  of  constraint.  In  the  absence  of  crack  interaction 
the  effective  Young's  modulus  (E^^.)  of  the  plate  was  described  by: 

E  ,,  =  E  (1  +  2TtNa2)_1  (5) 

ef  f  o 

where  Eq  is  Young's  modulus  of  the  crack-free  plate  and  a  is  the 
half-length  of  the  crack. 

An  energy  balance  analysis,  following  the  well-known  Griffith 
(1920)  theory,  showed  that  the  critical  temperature  range  (AT  ) 
of  cooling  the  plate  required  for  crack  instability  is: 

ATc  =  (2Y/a27iEa)1,/2(l  +  2nNa2)  (6) 

For  two  values  of  crack  density,  fig.  1  shows  AT^  as  a  func¬ 
tion  of  crack  size  indicated  by  the  solid  curve.  For  short  crack 

2 

length  such  that  27rNa  <<  l,  AT^  is  independent  of  crack  density, 

1/2 

and  is  inversely  proportional  to  a  .  For  high  crack  length 

2  3/2 

such  that  2uNa  >>  1,  AT^  is  proportional  to  a  ,  due  to  the 

increased  compliance  of  the  plate.  For  a  given  crack  density  this 
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implies  that  high  value.*;  of  AT  (i.e.,  high  thermal  stress  re¬ 
sistance)  can  be  achieved  by  increasing  crack  length.  At  first 
sight  this  conclusion  appears  contradictory.  However,  this  effect 
arises  because  the  failure  stress  is  controlled  by  a  single  crack 
only,  whereas  the  thermal  stress  is  governed  by  the  elastic  be¬ 
havior  of  the  plate,  which  in  turn  is  governed  by  the  collective 

effect  of  all  cracks.  A  minimum  in  AT  occurs  at  a  crack  size, 

c 

a  which  represents  the  minimum  crack  size  for  which  the  change 
m 

in  effective  Young's  modulus  becomes  significant.  At  a  given 
value  of  AT  the  cracks  are  unstable  between  two  values  of  crack 
length. 

At  AT  for  a  crack  with  initial  size,  a  >  a  crack  propaga- 
c  o  m 

tion  occurs  in  a  stable  mode.  For  increase  in  AT  to  values  greater 
than  AT^,  the  relation  between  AT^  and  crack  length  is  given  by  eq . 

For  an  initial  crack  size,  aQ<  a  ,  the  elastic  energy  release 
rate  exceeds  the  energy  required  to  create  the  new  crack  surfaces. 
This  will  cause  the  crack  to  propagate  in  an  unstable  (dynamic) 
mode.  Following  the  instability  at  ATc  the  crack  will  accelerate 
and  reach  maximum  velocity  when  it  reaches  a  size  equal  to  the  high 
er  value  of  critical  crack  size.  For  this  reason  the  crack  will 
continue  to  propagate  until  the  kinetic  energy  and  additionally 
released  strain  energy  is  converted  to  fracture  energy.  For  a  <<a 
the  final  crack  size  (a^)  which  results  from  this  mode  of  crack 
propagation  can  be  derived  to  be  (Hasselman,  1971)  : 
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a,  =  (47rNa  ) 
r  o 


(7) 


The  clotted  lines  in  fig.  1  show  the  value  of  Oj.  which  results 
from  the  unstable  crack  propagation  for  cracks  with  initial  size, 
a  <  a  .  As  indicated  by  fig.  1  and  eq.  7  the  final  crack  length 

o  m  7  r,  i  c 

is  an  inverse  function  of  the  initial  crack  length. 

The  results  shown  in  fig.  1  are  quantitatively  valid  only  for 
the  particular  mechanical  model  chosen,  and  assume  the  absence  of 
crack  interaction.  Nevertheless  on  the  basis  of  dimensional 
grounds  the  predicted  fracture  behavior  should  agree  qualitatively 
with  experimental  observations.  Data  in  support  of  this  latter 
conclusion  will  be  reviewed  briefly. 

Crack  sizes  in  brittle  materials  are  not  easily  established 
experimentally.  For  this  reason,  it  is  common  practice  to  monitor 
the  extent  of  crack  propagation  indirectly  by  the  measurement  of 
load-bearing  ability  (tensile  strength)  after  subjecting  the  cera¬ 
mic  to  thermal  shock.  In  order  to  provide  a  basis  on  which  to 
interpret  the  changes  in  load-bearing  ability,  fig.  2  shows  the 
crack  length  and  corresponding  changes  in  strength  following  ther¬ 
mal  shock  for  unstable  and  stable  crack  propagation. 

For  high-strength  materials  (i.e,  initial  crack  size,  a^  <  a  ) , 

the  unstable  mode  of  crack  propagation  results  in  a  discontinuity 

in  strength  at  AT^.  Following  arrest,  the  crack  is  sub-critical 

so  that  strength  for  AT  >  AT  remains  invariant  until  at  AT^  the 

c  c 
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crack  will  resume  propagation  in  a  stable  mode.  Experimental 
data  which  support  this  predicted  strength  behavior  is  presented 
in  a  number  of  studies  (Hasselman,  1970;  Gupta,  1972,  Mai  and 
Atkins ,  1975)  . 

Since  the  crack  length  following  thermal  stress  fracture  is 
inversely  proportional  to  the  crack  length  prior  to  fracture  ini¬ 
tiation,  the  strength  retained  is  expected  to  be  an  inverse  func- 
of  the  strength  prior  to  thermal  shock,  as  shown  experimentally 
by  Hasselman  (1970) . 

Low  strength  materials  with  a  >  a  will  show  stable  crack 

o  m 

propagation  and  a  monotonic  decrease  in  strength  as  shown  by 
Larson  et  al.  (1974). 

The  number  of  cracks  participating  in  the  fracture  process 
will  govern  whether  crack  propagation  will  occur  in  a  stable  or 
unstable  mode  confirmed  by  Larson  and  Hasselman  (1975). 

For  unstable  crack  propagation,  the  strength  loss  at  AT^  will 
be  an  inverse  function  of  the  number  of  cracks  which  form,  as  shown 
experimentally  by  Bertsch  et  al.  (1974). 

High  densities  of  small  cracks  should  result  in  very  high 
thermal  stress  resistance  as  confirmed  by  Rossi  (1969)  by  the 
comparative  thermal  shock  behavior  of  hot-pressed  magnesium  oxide 
and  composites  of  magnesium  oxide  with  a  dispersed  phase  of  tungsten. 

The  above  prediction  confirmed  by  experimental  evidence  sug¬ 
gests  that  at  least  qualitatively  the  thermal  stress  fracture  and 
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crack  propagation  behavior  of  brittle  ceramics  is  well  understood. 


The  possible  role  of  crack  interaction  will  be  examined  next. 

II.  EFFECT  OF  CRACK  INTERACTION 

Multiple  cracks  can  occur  in  an  almost  infinite  number  of  re¬ 
lative  orientation.  Clearly  an  exhaustive  analysis  of  all  possible 
interaction  effects  is  beyond  the  scope  of  the  present  study.  The 
specific  configuration  chosen  for  an  analysis  of  the  effect  of 
crack  interaction  consists  of  a  flat  plate  with  a  rectangular  array 
of  co-linear  and  columnar  cracks,  shown  in  fig.  3.  For  this  crack 
configuration,  Delameter  et  al.  (1975)  obtained  results  for  the 
stress  intensity  factors  and  effective  Young's  modulus  as  a  func¬ 
tion  of  co-linear  and  colimnar  crack  spacings.  By  uniaxially  con¬ 
straining  the  plate  from  free  thermal  expansion,  analogous  to  the 
mechanical  model  used  in  the  previous  study,  the  results  of  Dela¬ 
meter  et  al.  (1975)  can  be  used  to  analyze  the  effct  of  crack  inter¬ 
action  on  crack  stability  and  propagation  behavior  in  a  thermal 
stress  field.  For  this  purpose,  the  numerical  data  of  Delameter 
et  al.  were  curve-fitted  to  a  polynomial  which  gave  an  excellent 
fit  with  the  numerical  data  to  within  one  percent  for  the  range  of 
values  of  0  £  2a/d  £  0.8  and  0.2  £  2a/b  £  2.0.  All  the  figures  in 
this  section  are  limited  by  the  extent  of  the  numerical  results  of 


Delameter . 


The  effect  of  interaction  between  co-linear  cracks  should  be 


examined  first.  For  a  single  row  of  cracks,  fig.  4  shows  the 
critical  temperature,  ATc  required  for  crack  instability  as  a 
function  of  crack  size  (in  arbitrary  units)  for  three  values  of 
inter-crack  distance.  Since  for  a  single  row  of  co-linear  cracks, 
the  compliance  of  an  infinite  plate  is  not  affected,  AT^  always 
decreases  with  increasing  crack  length.  For  short  cracks  (a/d  < 

0.1)  such  that  interaction  affects  are  absent  the  log  AT^  -  log  a 
plot  has  a  slope  of  -1/2,  as  expected.  However,  as  a/d  approaches 
0.5,  the  value  of  the  slope  decreases  very  rapidly  until  at  a/d  = 

0.5  the  value  of  AT^  goes  to  zero.  This  condition  occurs  because 
at  a/d  =  0.5  all  crack  coalesce. 

Figure  5  shows  AT^  for  a  single  column  of  cracks.  For  wide 
crack  spacings  such  that  b/2a  >>  1,  the  value  for  AT£  corresponds 
to  the  value  for  a  single  crack.  The  log  ATC  vs.  log  a  curve 
has  a  slope  which  is  less  than  -1/2.  In  fact,  as  b/a  -*■  1,  the 
slope  becomes  positive,  strongly  indicative  of  the  positive 
effect  on  thermal  stress  resistance  due  to  the  interaction  bet¬ 
ween  columnar  cracks.  This  effect  is  not  related  to  the  decrease 
in  Young's  modulus  which  for  a  single  column  of  cracks  in  an  infinite 
plate  is  negligible,  but  is  due  to  the  "shielding"  between  cracks 
which  in  effect  decreases  the  stress  intensity  factor. 

To  show  the  combined  effect  of  co-linear  and  columnar  crack 
interaction,  fig.  6  shows  AT^  for  constant  co-linear  crack  spacing  d. 
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as  a  function  of  a/d  for  a  range  of  values  of  b/d.  In  effect, 
for  a  given  value  of  crack  size,  fig.  6  examines  the  effect  of 
columnar  inter-crack  spacing  on  AT^  for  a  plate  with  a  rectangular 
array  of  cracks.  The  data  show  that  AT^  is  well  in  excess  of  the 
value  for  a  single  row  of  cracks.  This  effect  is  due  to  the  com¬ 
bined  effect  of  the  "shielding"  between  columnar  cracks  and  the 
decrease  in  effective  Young's  modulus  with  decreasing  b/d  ratio. 

It  should  be  note*d,  however,  that  for  a  given  b/d  ratio,  AT^  goes 
througli  a  maximum  with  increasing  value  of  a/d,  and  approaches  zero 
as  a/d  -*  0.5.  This  behavior  is  in  contrast  with  the  results  shown 
in  fig.  1,  which  indicates  that  in  the  absence  of  crack  interaction 
ATc  increases  monotonical ly  with  increasing  crack  length. 

Figure  7  shows  similar  results  for  AT^  for  constant  inter- 
columnar  spacing  b  as  a  function  of  a/b  for  a  range  of  values  of 
d/b.  In  effect  for  a  given  crack  size,  this  shows  the  effect  of 
co-linear  crack  spacing  on  AT^.  The  results  indicate  that  with 
decreasing  value  of  d/b,  AT^  increases  over  the  value  of  the  single 
column.  This  effect  occurs  in  spite  of  the  fact  that  as  indicated 
in  fig.  4,  AT^  decreases  with  decreasing  co-linear  crack  spacing. 
The  actual  increase  in  AT^  is  due  to  the  fact  that  the  positive 
effect  of  decreasing  Young's  modulus  with  decreasing  value  of  d 
(i.e.,  increasing  crack  density)  on  ATC  is  greater  than  the 
negative  effect  on  ATC  due  to  decrease  in  co-linear  crack  spacing, 
d. 
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Conclusions  similar  to  those  which  can  be  drawn  from  the  data 

in  figs.  6  and  7,  are  also  evident  in  fig.  8,  which  shows  AT^  for 

-1  1/2 

constant  crack  density,  N  =  (bd)  ,  as  a  function  of  aN  for  a 

1/2  1/2 

range  in  values  of  dN  .As  long  as  aN  remains  small,  AT^ 
with  crack  interaction  exceeds  the  corresponding  value  for  the  ab¬ 
sence  of  crack  interaction  as  the  direct  result  of  the  decrease  in 

1/2 

effective  Young's  modulus.  However,  at  the  higher  values  of  aN 

*  1/2 

and  lower  values  of  dN  ,  as  the  result  of  increased  co-linear 

crack  interaction  AT  falls  below  the  value  of  AT  in  the  absence 

c  c 

of  crack  interaction. 

The  data  presented  so  far  concentrated  on  the  onset  of  crack 
propagation.  As  related  earlier,  the  extent  of  crack  propagation 
also  is  a  decisive  variable  in  assessing  the  thermal  shock  resistance 
of  brittle  ceramics.  The  interaction  between  cracks  also  is  ex¬ 
pected  to  affect  the  final  crack  size  ( a^) .  Values  for  a^-  were 
calculated  on  the  assumption  that  all  the  elastic  energy  released 
during  crack  propagation  is  converted  to  surface  fracture  energy. 

The  results  obtained  are  shown  in  fig.  9  for  a  number  of  values  of 
dN  of  fig.  8.  Primarily  as  the  result  of  the  columnar  crack 

interaction,  the  final  crack  length  is  less  than  the  corresponding 
one  in  the  absence  of  crack  interaction.  This  latter  conclusion 
is  valid  only  as  long  as  the  final  crack  length  remains  small  re¬ 
lative  to  the  co-linear  crack  spacing.  If,  however,  the  final 
crack  length  approaches  the  value  of  the  co-linear  crack  spacing. 
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an  interesting  effect  should  occur,  indicated  by  the  curve  for 
1/2 

dN  =2.  The  dotted  part  of  this  curve  which  falls  below  the 
crack  stability  curve,  AT^,  denotes  the  position  of  crack  arrest 

for  a  crack  with  initial  size  a  <  a  .  The  crack  is  sub-critical 

o  m 

only  as  long  as  it  falls  below  the  corresponding  solid  curve 
(AT^)  for  continued  propagation  in  a  stable  mode.  This  sub¬ 
criticality  disappears  at  a  value  of  final  critical  crack  size, 

a,  at  which  solid  curve  for  AT  and  the  dotted  curve  for  crack 
f  c 

arrest  intersect.  For  this  reason,  if  the  final  crack  length 
would  exceed  the  values  indicated  by  the  point  of  intersection 
of  the  dotted  and  the  solid  curves,  the  crack  will  again  be  un¬ 
stable  and  will  continue  to  propagate  until  complete  crack 
coalescence  occurs.  In  other  words  the  material  will  completely 
disintegrate.  The  same  effect  will  occur  for  stable  crack  propa¬ 
gation  whenever  crack  length  exceeds  the  value  corresponding  to 
the  peak  in  AT^  vs.  a  curve.  These  latter  results  contrast  with 
those  obtained  for  the  absence  of  crack  interaction,  which  do 
not  predict  this  fracture  mode. 

These  above  results  have  important  implication  for  the  reten¬ 
tion  of  load-bearing  ability  of  brittle  ceramics  subjected  to 
thermal  shock  and  in  which  strong  crack  interaction  effects  are 
present.  Such  strength  behavior  is  shown  schematically  in  fig. 

10.  Curve  (a)  results  from  the  extensive  propagation  of  an 
initially  very  short  crack  (in  a  high-strength  material)  which 
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due  to  co-linear  crack  Interaction  will  not  exhibit  crack  arrest. 


Curve  (b)  indicates  the  strength  behavior  of  a  material  in  which 
the  crack  arrests  to  be  subcritical  followed  by  stable  crack 
propagation  and  crack  coalescence  resulting  in  complete  disin¬ 
tegration  of  the  material.  Curve  (c)  corresponds  to  a  crack 
which  propagated  in  a  stable  mode  until  crack  coalescence  and 
complete  failure.  These  writers  are  not  aware  of  experimental 
data  which  demonstrate  this  behavior.  A  search  of  the  literature 
for  such  data  or  an  experimental  research  program  to  establish 
such  data  is  recommended  as  highly  valuable. 

III.  SUMMARY  AND  CONCLUSIONS 

The  effect  of  crack  interaction  on  the  initiation  of  thermal 
stress  fracture  and  resulting  crack  propagation  in  brittle  mate¬ 
rials  was  investigated  using  mechanical  model  consisting  of  a  flat 
plate  with  a  rectangular  array  of  co-linear  and  columnar  cracks. 
The  general  results  indicate  that  the  crack  interaction  in  com¬ 
bination  with  changes  in  the  effective  elastic  behavior  of  the 
plate,  has  a  beneficial  effect  on  thermal  stress  resistance. 

Only  when  the  crack  sizes  approach  the  co-linear  crack  spacing, 
thermal  stress  resistance  is  decreased.  Implications  of  these 
results  on  the  retention  of  the  load-bearing  ability  of  brittle 
ceramics  fractured  by  thermal  shock  are  discussed. 
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l'ig.  1.  Crack  stability  and  propagation  behavior  for  a  rigidly 
held  thermally  stressed  flat  plate  for  two  values  of 
density  of  non- interacting  cracks (Hasselman ,  1971) 

Fig.  2.  Crack  propagation  and  strength  behavior  for  (A)  and  (B) : 

unstable  and  (C)  and  (D) :  stable  mode  of  propagation  of 
non-interacting  cracks  under  conditions  of  thermal  shock 
(Larson  et.  al.,  1974) 

Fig.  3.  Flat  plate  with  rectangular  array  of  cracks  for  the  analy¬ 
sis  of  the  effect  of  crack  interaction  on  the  thermal 
fracture  of  brittle  ceramics 

Fig.  4.  Effect  of  colinear  crack-interaction  on  the  critical  tem¬ 
perature  difference  for  crack  instability 

Fig.  3.  Effect  of  columnar  crack-interaction  on  the  critical  tem¬ 
perature  difference  for  crack  instability 

Fig.  6.  Combined  effect  of  colinear  and  columnar  crack-interaction 
on  the  critical  temperature  difference  for  constant  co¬ 
linear  crack  spacing 

Fig.  7.  Combined  effect  of  colinear  and  columnar  crack-interaction 
on  the  critical  temperature  difference  for  constant  colum¬ 
nar  crack  spacing 

Fig.  8.  Combined  effect  of  colinear  and  columnar  crack-interaction 

on  the  critical  temperature  difference  for  crack  instability 
for  a  range  of  value  of  colinear  and  columnar  crack  spacing 
with  a  constant  value  of  crack  density,  N 

Fig.  9.  Effect  of  crack  interaction  on  the  extent  of  crack  propa¬ 
gation  following  unstable  crack  propagation 

Fig.  10.  Predicted  strength  behavior  of  brittle  ceramics  as  a  func¬ 
tion  of  the  severity  of  thermal  shock  for  stable  and  un¬ 
stable  propagation  of  interacting  cracks 
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Fig.  2.  Crack  propnRntion  and  strength  behavior  for  (A)  and-  (P):  unstable  and 
(C)  and  (D) :  stable  mode  of  propagation  of  non-internet ing  cracks 
under  conditions  of  thermal  shock  (Larson  et  nl.,  1974) 
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Fig.  4.  F.ffect  of  co-llncar  crack-interaction  on  the  critical  temporal  n  re  <11  f 
ference  for  crack  instability 
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Fig.  5.  Effect  of  columnar  crack-interaction  on  the  crilic.nl  temperature 
difference  for  crack  instability 


Fig.  6.  Combined  effect  of  colinear  and  columnar  crack-interaction  on  the  critical  temperature 
difference  for  constant  colinear  crack  spacing 


Fig.  8.  Combined  effect  of  colinear  and  columnar  crack-interaction  on  the  critical  tempera 
difference  for  crack  instability  for  a  range  of  values  of  colinear  and  columnar  cr 
spacing  with  a  constant  value  of  crack  density,  N 
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Fig.  9.  Effect  of  crack  Interaction  on  the  extent  of  crack  propagation 
following  unstable  crack  propagation 


unstable  propagation  of  interacting  cracks 
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STRESSES  DUE  TO  THERMAL  TRAPPING  IN  SEMI-ABSORBING  MATERIALS 
SUBJECTED  TO  INTENSE  RADIATION 


J.  P.  Singh,  J.  R.  Thomas,  Jr.,  and  D.  P.  H.  Hasselman 

Departments  of  Materials  and  Mechanical  Engineering 
Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  Virginia  24061 


ABSTRACT 

Analytical  results  are  presented  for  the  thermal  stresses  re¬ 
sulting  from  "thermal  trapping"  in  semi-absorbing  materials  sub¬ 
jected  to  symmetric  and  assymmetric  radiation  heating  and  convec¬ 
tive  cooling  with  finite  heat  transfer  coefficient,  h.  The  tran¬ 
sient  stresses  during  heat-up  in  both  cases  were  found  to  be  an 
inverse  function  of  the  heat  transfer  coefficient,  h,  and  increase 
monotonically  with  the  optical  thickness,  pa.  In  contrast,  the 
steady  state  stresses  were  independent  of  h  and  exhibited  a  maximum 
at  pa  =  1.3  for  symmetric  heating  and  at  pa  =  2  for  assymmetric 
heating  with  zero  stresses  at  pa  =  0  and  ».  For  the  symmetrically 
heated  plate,  the  transition  from  the  transient  to  the  steady-state 
condition  involved  a  reversal  in  the  sign  of  thermal  stress  at  any 
position.  For  the  assymmetrically  heated  plate  the  steady  state 
maximum  tensile  thermal  stresses  occur  at  the  position  of  the  high¬ 
est  temperature,  ie.  the  front  surface. 


INTRODUCTION 

Many  engineering  structures  such  as  aerospace  vehicles  and 
concentrating  solar  collectors  are  subjected  to  high  intensity  radi¬ 
ant  heat  fluxes.  In  order  to  withstand  the  high  temperature  in¬ 
volved,  such  structures  must  be  constructed  of  materials  with  high 
melting  ooint  and  good  mechanical  properties.  Unfortunately,  such 
materials  tend  to  be  brittle  and  are  susceptible  to  catastrophic 
failure  under  the  influence  of  transient  and  steady  state  thermal 
stresses  resulting  from  non-linear  temperature  distribution! >2. 
Therefore,  for  the  selection  of  the  optimum  material  and  reliable 


engineering  design,  it  is  necessary  to  analyze  the  thermal  stresses 
arising  from  radiation  heating  and  understand  the  variables  which 
control  their  magnitude. 

The  analysis  of  thermal  stress  failure  of  an  opaque  material 
subjected  to  black  body  radiation  was  performed  by  Hasselman^  some 
time  ago.  Subsequently,  the  effect  of  the  spectral  dependence  of 
absorption  coefficient  was  taken  into  account  to  include  materials 
completely  transparent  below  and  opaque  above  a  given  wavelength^* 

For  many  materials,  however,  the  incident  radiation  may  be  absorbed 
throughout  the  thickness  and  depending  upon  the  absorption  propert¬ 
ies,  it  may  result  in  so  called  "thermal  trapping"  effect  in  which 
the  interior  temperature  of  the  body  may  exceed  the  ambient  tempera¬ 
ture.  Thermal  stresses  in  a  flat  plate  due  to  this  type  of  heat 
transfer  were  analyzed  recently  by  Hasselman  et.  al.5  for  symmetric 
radiation  heating  with  convective  cooling  for  two  limiting  values  of 
heat  transfer  coefficient,  h , equal  to  zero  and  infinity.  Subsequent 
ly  two  follow-up  studies  by  Thomas  et.  al.^  and  Singh  et.  al.7 
included  the  case  of  symmetric  and  assymmetric  radiation  heating 
and  convective  cooling  with  finite  values  of  h. 

The  purpose  of  this  paper  is  to  review  and  report  the  main 
results  for  stresses  in  the  latter  two  studies.  For  further  details 
the  reader  is  referred  to  the  original  papers. 6, 7 

ANALYSIS 

The  analysis  of  thermal  stresses  in  the  flat  plate  first  re¬ 
quired  obtaining  the  temperature  distribution,  which  was  subsequent¬ 
ly  utilized  to  calculate  the  thermal  stresses.  In  this  analysis  the 
reflectivity  of  the  material  was  assumed  to  be  low  and  hence  the 
effect  of  multiple  internal  reflections  was  neglected.  It  was  also 
assumed  that  the  maximum  value  of  the  stress  in  the  plate  will  be 
reached  before  the  temperature  becomes  sufficiently  high  that  the 
effect  of  re-emission  must  be  taken  into  account^.  The  validity  of 
this  latter  assumption  can  easily  be  substantiated  with  numerical 
examples5«6.  Finally  it  was  assumed  that  the  emissivity  and  absorp¬ 
tivity  of  the  plate  are  independent  of  wavelength.  However,  the 
spectral  dependence  of  these  properties  can  easily  be  incorporated 
in  the  analysis,  if  needed. 

The  general  analytical  approach  along  with  the  heating  and 
cooling  conditions  are  as  follows: 

A.  Symmetric  Radiative  Heating  and  Convective  Cooling 

The  rate  of  internal  heat  generation  (g"')  within  the  plate 
subjected  to  symmetric  normally  incident  radiation  q0  is  given  by: 


cosh(yx) 


(1) 


r  » 


in  o  -pa 
g  =  2ueo  e 
o 

where  \i  is  the  absorption  coefficient,  c  is  the  emissivity  (e=l-r, 
where  r  is  the  reflectivity),  a  is  the  half  thickness  of  the  plate 
and  x  is  the  through-the-thickness  coordinate  with  x=0  at  the  center 
of  the  plate. 

The  temperature  solutions  were  obtained  by  solving  the  differ¬ 
ential  equation^: 

-g+8"'(x)/k=  (1/K)|f  (2) 

where  k  is  the  thermal  conductivity,  t  is  the  time  and  k  is  the 
thermal  diffusivity. 

The  initial  and  boundary  conditions  were: 

T(x,0)  =  Tq;  (0,t)  =  0 

f(a,t|  -  -  £  {T(a, t)  -  To) 

3T  h 

^(-a,t)  =~  {T(^,t)  -  To} 

where  h  is  the  convective  heat  transfer  coefficient. 

The  solution  of  temperature  T,  thus  obtained  from  differential 
equation  (2),  was  utilized  to  obtain  the  expressions  for  thermal 
stress  fromlO: 
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which  yields: 
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where  Xn  are  the  roots  of  the  transcendental  equation 
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B.  Assymmetric  radiation  Heating  and  Convective  Cooling 


In  this  case,  the  radiation  a  is  normally  incident  to  only 
one  face  of  the  elate  at  x  =  -a  an8  the  plate  is  cooled  by  con¬ 
vection  at  x  =  a.  The  intensity  of  internal  heat  generation  (g"') 
within  the  plate  can  be  shown  to  be: 


g’"(x) 


Vcqoe 


-lj(a+x) 


(8) 


The  initial  and  boundary  conditions  are: 

’  V  I  ■ 0 


(9) 


-kJ(a,t)=h[T(a,t)-To]  (10) 

The  temperature  solutions  were  obtained  by  solving  the  differential 
equation  (2)  subjected  to  the  conditions  (9)  and  (10).  Similar  to 
case  A,  the  temperature  solution,  thus  obtained,  was  utilized  to 
find  out  the  expression  for  thermal  stress.  The  thermal  stress 
in  this  case  is  given  by: 
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where  X  are  the  roots  of  the  transcendental  equation 
n 


X  tan  ( 2 X  a)  =  h/k 
n  n 


and 


B  (t)  =  G  [1  -  e_KXn"t]/icX  2 
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NUMERICAL  RESULTS  AND  DISCUSSION 

For  the  convenience  of  data  presentation,  the  results  will  be 
presented  in  non-dimensional  form  as  follows: 

The  non-dimensional  time: 

t*  =  Kt/a2  (12) 

the  non-dimensional  temperature: 

T*  =  (T  -  Tq)  k/eqQa  (13) 

and  the  non-dimensional  stress: 

A 

a  -  o(l-v)k/aEeqQa  (14) 

For  simplicity,  the  quantities  t'*‘,  T  and  a*  hpreafter  will  be  re¬ 
ferred  to  as  time,  temperature  and  stress,  respectively. 

Figure  1  shows  the  thermal  stresses  in  the  surface  and  the 
center  of  a  symmetrically  heated  plate  as  a  function  of  time  for  a 
range  of  values  for  h/k  with  an  arbitrarily  chosen  value  of  the 
optical  thickness  pa  =  3.  Stresses  reach  the  steady  state  value 
after  an  initial  transient.  For  a  given  optical  thickness  pa,  the 
value  of  the  steady  state  stress  is  independent  of  h.  This  is  so 
because  the  steady  state  temperature  profile  within  the  plate  is 
independent  of  the  heat  transfer  coefficient  h,  even  though  the 
magnitude  of  the  temperature  for  a  given  profile  is  an  inverse  func¬ 
tion  of  the  heat  transfer  coefficient.  The  values  of  the  maximum 
transient  stress,  however,  decrease  with  increasing  h  in  contrast 
to  many  literature  solutions  for  transient  thermal  stresses  in  con¬ 
vection  heat  transfer  environment  where  the  thermal  stresses  in¬ 
crease  with  increasing  value  of  h. 

Another  noteworthy  point  in  figure  1  is  the  change  in  the  sign 
of  the  thermal  stresses  in  the  plate  with  time  for  finite  values  of 
h/k.  At  small  values  of  time  the  surface  stress  is  compressive 
whereas  the  stress  in  the  center  is  tensile.  However,  at  longer 
times  the  stresses  reverse  their  sign  and  become  tensile  and  com¬ 
pressive  in  the  surface  and  the  center,  respectively.  This 
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Fig.  1.  Maximum  thermal  stress  as  a  function  of  time  in  a  partially 
absorbing  flat  plate  symmetrically  heated  by  thermal  radia¬ 
tion  and  cooled  by  convection  for  various  values  of  h/k 
(cm.-l)  with  optical  thickness  pa=*3  and  a=l  cn. 

phenomenon  is  a  result  of  the  non-uniform  heat  generation  inside 
the  plate  and  the  associated  "thermal  trapping"  effect.  Due  to  a 
larger  portion  of  the  incident  radiation  being  absorbed  in  the 
surface,  initially  the  surface  is  at  higher  temperature  than  the 
interior.  During  this  time  period,  the  heat  loss  from  the  surface 
by  convection  is  very  small  because  of  the  low  surface  temperature 
in  combination  with  finite  value  of  h  and  the  heat  flows  primarily 
from  the  surface  to  the  interior.  This  results  in  a  temperature 
distribution  which  causes  a  compressive  stress  in  the  surface  and 
tensile  stress  in  the  center.  However,  after  a  long  period  of 
time  (t  -*■  »)  the  surface  gets  hot  with  an  appreciable  amount  of  heat 
loss  by  convection  which  eventually  will  be  equal  to  the  incident 
radiative  heat  flux  at  steady  state.  Under  this  condition,  due  to 
the  internal  heat  generation,  the  interior  of  the  plate  will  be  at 
higher  temperature  than  the  surface  resulting  in  a  tensile  stress 
in  the  surface  and  a  compressive  stress  in  the  center  of  the  plate. 
However,  as  h/k  ->  »,  no  stress  reversal  in  the  surface  takes  place 
because  the  surface  temperature  remains  constant  (equal  to  the  am¬ 
bient  temperature  )  at  all  times  resulting  in  tensile  surface  stress 
for  t  >  0. 

For  the  symmetrically  heated  plate,  figures  2  and  3  show  the 
maximum  values  of  transient  and  steady  state  stress  as  a  function 
of  optical  thickness,  pa.  It  is  interesting  to  note  that  the 
maximum  value  of  the  transient  stress  increase  with  optical  thick¬ 
ness  whereas  the  steady  state  stress  exhibits  a  maximum  at  pa*1.3 
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OPTICAL  THICKNESS,  p.0 

Fla,.  2.  T'agnitude  of  maximum  transient  thermal  stresses  as  a  func¬ 
tion  of  ootical  thickness  pa  in  partially  absorbing  flat 
plate  symmetrically  heated  by  thermal  radiation  and  cooled 
by  convection  for  various  values  of  h/k(cm.-l)  with  a=l  cm 
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Fig.  3.  Magnitude  of  maximum  steady  state  thermal  stresses  as  a 

function  of  optical  thickness  na  in  a  partially  absorbing 
flat  plate  sy  v  etrically  heated  by  thermal  radiation  and 
cooled  by  convection  for  various  values  of  h/k(cm."l)  with 
a~l  cm. 


and  approaches  zero  at  pa=Q  and  For  pa=Qv  no  heat  is  absorbed, 
resulting  in  zero  stress.  However,  for  pa-*-00,  although  the  heat  is 
transmitted  through  the  plate,  the  steady  state  temperature  profile 
is  linear  resulting  in  zero  stress  at  t=°°.  At  pa=1.3,  the  non¬ 
uniformity  in  the  internal  heat  generation  and  the  associated  ther¬ 
mal  trapping  effects  are  maximum  resulting  in  the  maximum  stress. 

Tor  an  assymmetrically  heated  plate,  figure  A  shows  the  values 
of  maximum  tensile  thermal  stress  as  a  function  of  time  for  various 
values  of  h/k  with  an  arbitrarily  chosen  value  of  pa=3.  Similar 
to  figure  1,  the  stress  reaches  a  constant  steady  state  value  in¬ 
dependent  of  h  and  the  value  of  maximum  transient  stress  is  an  in¬ 
verse  function  of  h. 


Fig.  A.  Ifaxinum  tensile  thermal  stress  as  a  function  of  time  in  a 
partially  absorbing  flat  plate  assymmetrically  heated  in 
front  by  thermal  radiation  and  cooled  at  the  rear  surface 
by  convection  for  various  values  of  h/k(cm.~l)  with  pa=3 
and  a=l  cm. 

Figure  5  shows  the  spatial  distribution  of  the  steady  state 
tensile  stress  (t  ->  °°)  in  the  assymmetrically  heated  plate.  It 
is  important  to  note  that  the  maximum  tensile  stress  occurs  in  the 
front  face  (x  =  -a)  of  the  plate  which  is  obviously  at  the  highest 
temperature.  This  initially  unexpected  result  was  shown  to  be  due 
to  a  concave  downward  temperature  distribution^  as  the  direct  re¬ 
sult  of  an  assymmetrically  non-uniform  internal  heat  generation. 
This  temperature  distribution  contrasts  with  the  usual  concave  up¬ 
ward  temperature  distribution  found  in,  many  cases  where  the  heat 
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Fig.  5.  Spatial  distribution  of  steady  state  thermal  stresses  in  a 
partially  absorbing  flat  plate  assymmetrically  heated  in 
front  by  thermal  radiation  and  cooled  at  the  rear  surface 
by  convection  for  various  values  of  pa  with  h/k(cm.-^) 
finite  and  a  =  1  cm. 

is  exchanged  only  at  the  outside  boundaries  of  the  structure.  For 
this  latter  case,  the  maximum  tensile  stresses  freouently  occur  at 
or  near  the  position  of  the  lowest  temperature. 

Figure  6  summarizes  the  values  of  maximum  tensile  transient 
and  steady  state  thermal  stress  as  a  function  of  pa  in  the  assym- 
metrically  heated  plate.  Similar  to  the  observation  made  for 
symmetric  heating  case.  (Figures  2  and  3)  the  steady  state  stress 
exhibits  a  maximum,  which  occurs  at  pa=2  because  of  the  maximum 
thermal  trapping  effect.  However,  the  maximum  transient  stress  for 
finite  value  of  h  is  a  monotonically  increasing  function  of  pa  and 
will  fall  between  the  stress  olots  for  h=0  and  h= °°  (Figure  6). 

From  Figures  (4)  and  (6),  it  is  obvious  that  for  10.7  <  pa  <  <», 
the  maximum  transient  stress  is  greater  than  the  corresponding 
maximum  steady  state  stress  and  for  0  <  pa  <  10.7,  the  converse  is 
true. 

A  comparison  of  the  results  for  thermal  stresses  in  the  symmet¬ 
rically  and  assymmetrically  heated  flat  plates  indicates  that  the 
nature  of  the  stress  profiles  in  both  cases  are  oualita tively  similar 
(Figures  1-6) .  However,  the  magnitude  of  the  maximum  stress  in  case 
of  symmetric  heating  is  higher  than  that  in  case  of  assymmetric 
heating.  The  lower  value  of  stress  in  case  of  assymmetric  heating 
is  due  to  the  fact  that  the  non-uniform  thermal  expansion  can  be 


Fig.  6.  '  agnitude  of  maximum  tensile  thermal  stresses  as  a  function 

of  pa  in  a  partially  absorbing  flat  plate  assynmetrically 
heated  in  front  by  thermal  radiation  and  cooled  at  the  rear 
surface  by  convection  with  h/k(cm.~l)  finite  and  a  =  1  cm. 

partly  accommodated  by  bending  v  ich  is  not  the  case  in  symmetric 
heating. 

The  above  discussions  have  the  following  implications:  First, 
it  is  important  to  note  that  although  the  steady  state  stress  is 
zero  at  pa=0  and  °°,  the  magnitude  of  stress  could  be  significant 
for  the  intermediate  values  of  pa.  Therefore,  it  is  preferable  to 
have  a  transparent  material  to  avoid  thermal  stress  fialure.  How¬ 
ever,  this  criterion  for  material  selection  may  not  be  compatible 
with  other  performance  criteria.  For  example,  for  high  efficiency 
of  the  solar  collectors,  it  is  desirable  to  have  high  values  of  the 
absorption  coefficient  and  emissivity  to  keep  the  amount  of  trans¬ 
mitted  and  reflected  radiation  to  a  minimum.  Clearly  a  compromise 
between  high  collector  efficiency  and  high  thermal  stress  resistance 
will  have  to  be  made. 

Finally,  the  observation  of  the  maximum  tensile  thermal  stress 
in  the  hottest  part  (front  face)  of  the  assymmetrically  heated  plate 
could  be  very  critical  for  failure  by  the  stress  corrosion  mechanism. 
Stress  corrosion,  which  could  be  promoted  by  environmental  condi¬ 
tions,  is  generally  a  thermally  activated  process  and  is  expected 
to  increase  with  stress  and  temperature.  Therefore,  in  addition 
to  having  low  thermal  stress,  it  is  important  to  maintain  low  tem¬ 
perature  levels  in  the  plate  to  ensure  high  thermal  stress  resis¬ 
tance  and  long  life.  Thus,  care  nust  be  taken  in  selection  of 
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material  and  heat  transfer  conditions  in  order  to  keep  the  critical 
regimes  at  low  temperature.  It  is  quite  obvious  then,  that  the  in¬ 
cidence  of  thermal  stress  failure  in  a  plate  subjected  to  radiation 
heating  can  be  minimized  by  proper  selection  of  material  in  combina¬ 
tion  with  the  optimum  choice  of  the  optical  thickness  and  heat  trans¬ 
fer  conditions. 
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ABSTRACT 

The  results  are  presented  for  the  effect  of  spatially  varying 
thermal  conductivity  on  the  tensile  thermal  stress  developed  in  a 
solid  and  a  hollow  circular  cylinder  subjected  to  different  heating 
conditions.  It  is  shown  that  the  maximum  tensile  thermal  stress  in 
brittle  ceramics  can  be  reduced  significantly  by  redistributing  the 
temperature  profile  using  (a)  a  spatial  variation  in  thermal  conduc¬ 
tivity,  (b)  a  spatial  variation  in  pore  content  which  in  turn  changes 
the  density,  thermal  conductivity  and  modulus  of  elasticity  and  (c) 
by  considering  the  effect  of  temperature  on  the  thermal  conductivity 
and  specific  heat.  Possible  methods  for  creating  such  variations  in 
the  material  properties  are  discussed. 


INTRODUCTION 

The  design  of  engineering  structures  or  components  operating 
under  non-iso thermal  conditions  must  take  into  account  the  presence 
of  thermal  stresses  which  arise  from,  constrained  non-uniform  ther¬ 
mal  expansion  in  addition  to  the  normal  operating  stresses.  For 
many  structures  or  components  operating  at  high-temperatures  the 
transient  thermal  stresses  during  heat-up  or  cool-down  can  exceed 
the  normal  operating  stresses  by  a  considerable  margin.  Weakening 
of  the  structure  (or  complete  failure)  due  to  the  thermal  stresses 
should  be  avoided  to  assure  long-term  reliable  operation.  This  is 
particularly  critical  for  brittle  ceramic  materials  in  which  because 
of  the  absence  of  non-linear  stress-relieving  mechanisms  and  low 
fracture  toughness,  thermal  stress  failure  can  be  highly  catastrophic 
with  corresponding  disastrous  consequences  for  continued  satisfac- 


tory  performance. 

A  number  of  options  are  available  to  the  designer,  working  in 
conjunction  with  the  materials  technologist  for  reducing  the  pro¬ 
bability  of  failure  due  to  thermal  stresses.  First,  for  many  heat 
transfer  conditions,  a  redesign  involving  a  change  in  the  geometry 
or  a  reduction  in  the  dimensions  of  the  structure  or  component  will 
reduce  the  magnitude  of  the  thermal  stresses.  This  approach,  how¬ 
ever,  usually  is  incompatible  with  other  performance  criteria  im¬ 
posed  on  the  structure  or  component.  Secondly,  the  magnitude  of 
the  thermal  stresses  can  be  reduced  by  reducing  the  rate  of  heat 
transfer.  Again,  with  few  exceptions,  this  approach  is  impractical, 
as  the  rates  of  heat  transfer  are  beyond  the  control  of  the  designer. 
A  third  approach  is  to  select  the  materials  of  construction  of  the 
structure  or  component  on  the  basis  of  figures-of-merit  for  thermal 
stress  resistance  [1],  which  reflect  the  optimum  combination  of  the 
relevant  material  properties  such  as  the  coefficient  of  thermal  ex¬ 
pansion,  Young's  modulus,  thermal  conductivity,  etc.  It  should  be 
noted,  however,  that  thermal  stress  failure  in  even  the  optimum 
materials  selected  in  this  manner,  is  not  necessarily  avoided,  as 
these  figures-of-merit  merely  compare  the  relative  thermal  stress 
resistance  of  potential  materials  of  construction.  Furthermore, 
values  of  material  properties  which  maximize  thermal  stress  resis¬ 
tance  may  not  be  compatible  with  the  property  values  required  for 
optimum  performance  of  the  structure  or  component  in  its  intended 
application. 

A  fourth  approach  to  improving  thermal  stress  resistance  is 
based  on  the  fact  that  the  magnitude  of  the  thermal  stresses  in  a 
structure  or  component  not  only  is  a  function  of  the  magnitude  of 
the  temperature  differences  encountered  within  the  structure  or  com- 
oonent,  but  also  is  a  function  of  the  temperature  distribution.  In 
principle,  therefore,  it  should  be  possible  to  decrease  the  magnitude 
of  thermal  stress  for  a  given  heat  transfer  condition  by  modifica¬ 
tions  of  the  temperature  distribution.  Since  the  distribution  of 
the  temperature  within  a  structure  or  component  is  a  function  of 
the  thermal  conductivity,  modifications  in  those  distributions  can 
be  achieved  by  incorporating  a  spatially  varying  thermal  conducti¬ 
vity.  A  number  of  studies  carried  out  by  the  present  writers  and 
co-workers  [2, 3, 4, 5]  have  shown  that  significant  reductions  in 
the  magnitude  of  the  tensile  thermal  stresses  in  brittle  ceramic 
components  can  be  achieved  by  decreasing  the  thermal  conductivity 
in  the  hotter  parts  of  the  structure.  These  results  were  obtained 
for  steady-state  heat  flow  in  a  hollow-cylinder  as  well  as  in  the 
transient  heating  of  solid  circular  cylinders.  The  reductions  in 
the  tensile  stresses  were  obtained  not  only  by  a  direct  spatial 
variation  of  the  thermal  conductivity  but  also  indirectly  by  a 
spatial  variation  in  density  (which  affects  Young's  modulus  also), 
and  by  the  negative  temperature  dependence  of  the  thermal  conduc¬ 
tivity  of  dielectric  materials  in  which  the  conduction  of  heat 
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occurs  primarily  by  phonon  transport.  It  is  the  purpose  of  this 
presentation  to  give  an  overview  of  these  results. 


CALCULATIONS,  RESULTS  AND  DISCUSSION 


A.  Constant  Heat  Flow 


Spatially  Varying  Thermal  Conductivity.  The  geometry  selected 
for  the  analysis  of  the  effect  of  a  spatially  varying  thermal  con¬ 
ductivity  under  conditions  of  steady-state  heat  flow  consisted  of 
an  infinitely  long  hollow  concentric  cylinder  with  inner  and  outer 
radius  a  and  b,  resp.  The  spatial  variation  of  the  thermal  con¬ 
ductivity  (K) ,  for  analytical  convenience  was  selected  ass 

K  -  K  U  +  f  (a/ r)n}_1  '  (1) 

o  o 

for  radially  outward  heat  flow  and 

K  =  K  {1  +  f  (r/b)n}-1  (2) 

o  o 

for  radially  inward  heat  flow. 

In  eqs.  1  and  2,  K0  represents  the  spatially  uniform  thernal 
conductivity  and  f0  and  n  are  constants.  The  general  form  of  eqs. 

1  and  2  is  such  that  the  thermal  conductivity  is  lowest  in  the  re¬ 
gions  of  the  cylinder  which  are  at  the  highest  temperature. 

From  eqs.  1  and  2  and  the  general  condition  for  steady-state 
heat  flow  in  a  radial  coordinate  system  [6],  the  temperature  dis¬ 
tribution  within  the  cylinder  was  derived  in  terms  of  the  tempera¬ 
ture  difference  (AT)  across  the  wall. 

For  brittle  materials  the  maximum  tensile  stress  is  of  primary 
interest.  In  a  hollow  cylinder  undergoing  steady-state  heat  flow, 
the  radial  and  hoop  stresses  are .exceeded  by  the  longitudinal  stress, 
which  can  be  derived  from  [7]: 

°z  =  (1-v)  {b2-a2  {  rTdr  "  (3) 

where  n  is  the  coefficient  of  thermal  expansion,  E  is  Young's 
modulus,  v  is  Poisson's  ratio  and  T  is  the  temperature.  For  the 
purpose  of  this  analysis,  a,  E  and  v  were  assumed  to  be  independent 
of  position  or  temperature. 

Figures  1  and  2  show  the  maximum  value  of  the  tensile  longi¬ 
tudinal  thermal  stress  at  r  =  b  for  radially  outward  heat  flow  and 
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Fig.  1.  Maximum  tencile  iJiemal  stress  in  a  hollow  cylinder  under 
steady-state  outward  radial  heat  flow  and  spatially  vary¬ 
ing  thermal  conductivity  (equation  1) .  Dotted  lines  indi¬ 
cate  values  for  unifom  spatial  distribution. 

r  =  a  for  the  radially  inward  flow  for  three  values  of  the  ratio  of 
outer  to  inner  radius  (b/a) .  Included  in  the  figure  are  the  values 
of  stress  for  the  soatially  constant  thermal  conductivity.  Compari¬ 
son  of  these  stress  values  shows  that  a  spatially  varying  thermal 
conductivity  can  lead  to  a  significant  reduction  in  the  magnitude 
of  thermal  stress.  Since  an  order  of  magnitude  variation  in  thermal 
conductivity  of  materials  can  be  easily  found,  these  results  should 
be  of  practical  significance. 

Spatially  Varying  Density.  The  presence  of  a  oore  phase  in 
solid  materials  has  a  pronounced  effect  on  many  material  properties 
including  the  thermal  conductivity.  In  Principle,  therefore  a 
spatially  varying  thermal  conductivity  can  be  achieved  by  incorpo¬ 
rating  a  spatially  varying  pore  concentration  in  the  material.  A 
further  advantage  is  achieved  by  the  corresponsing  decrease  in 
Young's  modulus.  Unfortunately,  the  tensile  strength  is  also  af¬ 
fected  adversely  by  the  pore  phase  which  would  lead  to  a  decrease 
in  the  thermal  stress  resistance.  However,  this  effect  can  be 
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Fig.  2.  I!aximum  tensile  thermal  stress  in  a  hollow  cylinder  under 
steady-state  inward  radial  heat  flow  and  spatially  varying 
thermal  conductivity  (equation  2) .  Dotted  lines  indicate 
values  for  uniform  spatial  distribution. 

avoided  by  distributing  the  pore  phase  such  that  full  density  is 
maintained  at  the  position  of  maximum  tensile  stress  and  having  as 
high  a  porosity  as  practically  possible  in  the  region  at  higher 
temperature,  which  generally  are  in  a  state  of  compressive  stress. 

Accordingly,  the  porosity  distribution  for  radially  outward 
heat  flow  was  arbitrarily  taken  as: 

P  =  ?Q{(b-  r)  / (b  -  a)  }n  (4) 

and  for  radially  inward  heat  flow  as: 

?  =  PQ{(r  -  a) / (b  -  a)}°  (5) 
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where  P0  is  the  porosity  at  r  =  a  and  r  =  b  for  radially  outward  and 
inward  heat  flow,  resp.  and  n  the  porosity  distribution  coefficient. 
Eqs.  4  and  5  indicate  that  P  «=  0  at  the  position  of  maximum  tensile 
stress,  i.e.,  r  =  b  and  r  =  a  for  the  radially  outward  and  inward 
heat  flow,  respectively. 

The  effect  of  the  pore  phase  on  the  thermal  conductivity,  K 
and  Young's  modulus,  E  was  chosen  to  be  given  by  the  well-known 
empirical  equation  [8,9]: 

K,  E  =  K  ,  E  e"bp  (6) 

o  o 

where  Kq  and  E0  represent  the  values  of  the  thermal  conductivity 
and  Young's  modulus  for  the  non-porous  material,  p  is  the  fractional 
porosity  and  b  is  a  constant  taken  as  3  and  4  for  K  and  E,  respec¬ 
tively.  Poisson's  ratio  was  assumed  to  remain  unaffected  by  the 
pore  phase.  In  view  of  the  analytical  complexity  of  the  variation 
of  E  and  K  as  expressed  by  eqs.  4,  5  and  6,  the  thermal  stresses 
were  obtained  by  a  finite  element  computer  program  compiled  for 
this  purpose. 

For  cadially  inward  heat  flow.  Fig.  3  shows  the  maximum  value 
of  longitudinal  tensile  stress  as  a  function  of  the  constant  n  in 
eqs.  4  and  5  for  three  values  of  the  ratio  of  the  outer  to  inner 
radius  with  PQ  =  0.3.  Included  in  the  figure  is  the  value  of  stress 
corresponding  to  P(r)  =  0.  Similar  results  were  obtained  for  radial¬ 
ly  outward  heat  flow.  These  results  indicate  that  since  the  strength 
of  the  material  remains  unaffected  in  the  critical  regions,  signi¬ 
ficant  increase  in  thermal  stress  resistance  is  achieved.  Additional 
calculations  showed  that  the  reduction  in  the  magnitude  of  the 
maximum  tensile  stress  to  a  good  approximation  was  a  linear  function 
of  the  value  of  P0- 

Of  interest  to  note  in  Fig.  3  is  that  the  value  of  stress  is 
relatively  independent  of  the  value  of  n.  This  occurs  because  the 
maximum  reduction  in  stress  due  to  the  individual  effects  of  the 
spatially  varying  porosity  on  Young's  modulus  and  thermal  conduc¬ 
tivity  occur  at  different  values  of  n.  Minimizing  the  stresses  by 
reducing  Young's  modulus  requires  a  value  of  n  -*•  0.  This  implies 
a  uniform  pore  phase,  which  results  in  a  uniform  thermal  conducti¬ 
vity  with  on  reduction  in  the  thermal  stresses.  As  far  as  a  spatial 
variation  in  thermal  conductivity  is  concerned,  optimum  reduction 
in  the  thermal  stresses  is  achieved  by  values  of  n  >>  1  such  that 
a  reduction  in  thermal  conductivity  is  achieved  in  the  immediate 
regions  of  the  tube  surface  at  the  higher  temperature. 
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Fig.  3.  T'taxinun  tensile  thermal  stress  as  a  function  of  the  ^orosity 
distribution  coefficient  n  in  a  hollow  cylinder  under  stea¬ 
dy-state  radially  inward  heat  flow  with  spatially  varying 
pore  content  (equation  5)  with  P  (r  =  b)  =  0.3. 


B.  Transient  Heating 

An  infinitely  long  solid  circular  cylinder  of  radius  a  was 
selected  as  an  appropriate  geometry  for  the  analysis  of  the  effect 
of  a  spatially  varying  thermal  conductivity  on  the  magnitude  of 
thermal  stress  under  conditions  of  transient  heating.  The  specific 
heating  conditions  consisted  of  subjecting  the  cylinder  initially 
at  thermal  equilibrium  at  a  temperature  T0,  instantaneously  to  an 
ambient  temperature  of  Too.  The  transfer  of  heat  was  considered  to 
be  Newtonian  according  to: 

K(a) (3T/3r(a))  +  h  (T(a)  -  T  )  =  0  (7) 

03 

where  h  is  the  heat  transfer  coefficient  and  T(a)  is  the  instanta¬ 
neous  surface  temperature. 

Calculations  of  the  thermal  stresses  were  carried  out  with  a 
finite  element  program  which  gave  excellent  agreement  with  stress 
values  for  spatially  uniform  thermal  conductivity  obtained  analy¬ 
tically  by  Jaeger  [IQ]. 
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The  maximum  tensile  stresses  were  of  primary  interest  which  for 
the  present  heating  conditions  are  the  hoop  stresses  at  r  =  0. 

For  convenience,  the  numerical  results  are  reported  in  terms 
of  the  non-dimensional  thermal  stress,  a *  =  o(l-v) /aEAT,  the  non- 
dimensional  time,  t*  =  Kt/a2,  the  non-dimensional  radial  coordinate, 
r/a  and  the  Biot  number,  B0  =  ah/!<0  in  which  AT  =  T®  -  T0,  ic  is  the 
thermal  dif fusivity  and  a  is  the  coefficient  of  expansion.  The  Biot 
number  is  expressed  in  terms  of  the  uniform  thermal  conductivity, 
prior  to  spatial  modification. 

Spatially  Varying  Thermal  Conductivity.  The  radial  dependence 
of  the  thermal  conductivity  was  chosen  to  be  of  the  form: 

K(r)  =  K0[l  -  C(r/a)n]  (8) 

where  n  and  C(0  <  C  <  1)  are  constants.  Equation  8  indicates  that 
at  r  =  9,  K  =  K0  for  any  value  of  C  and  n.  The  lowest  value  of 
thermal  conductivity  occurs  at  r  =  a. 

Fig.  A  shows  tyoical  results  for  the  time  dependence  of  the 
maximum  tensile  hoop  stress  at  r  =  0  for  Bo  =  5  and  a  range  of 
values  of  C  and  n.  The  peak  values  of  the  stress  for  a  spatially 
varying  thermal  conductivity  are  considerably  less  than  the  cor¬ 
responding  values  for  the  constant  thermal  conductivity,  included 
in  Fin.  k  for  comparison. 


Fig.  A.  Maximum  tensile  thermal  stress  as  a  function  of  tire  at  the 
center  of  a  solid  circular  cylinder  subjected  to  sudden 
convective  heating  for  various  soatial  distributions  of 
the  thermal  conductivity  (equation  8) . 
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Figures  5a  and  5b  show  the  spatial  variation  of  the  tempera¬ 
tures  and  stresses,  resp.,  at  the  time  of  maximum  tensile  stress 
for  the  same  values  of  C  and  n  as  shown  in  Fig.  4.  These  results 
indicate  that  the  decrease  in  tensile  stress  at  r  =  0  is  accompanied 
by  a  corresponding  increase  in  the  temperatures  and  comnressive 
stresses  in  the  surface  'r  =  a). 


NON-DIMENSIONAL  DISTANCE,  r/a 


Fig.  5a.  Spatial  distribution  of  the  transient  temperature  at  the 
time  of  maximum  tensile  stress,  in  a  solid  circular  cy¬ 
linder  subjected  to  sudden  convective  heating  for  various 
spatial  distributions  of  the  thermal  conductivity  (eoua- 
tion  8) . 

Figure  6  shows  the  maximum  values  of  tensile  hoop  stress  at 
r  =  0  as  a  function  of  n  for  a  range  of  values  of  C.  Significant 
reductions  in  the  magnitude  of  stress  can  be  achieved,  especially 
at  the  higher  values  of  C. 

Figure  7  compares  the  maximum  values  of  tensile  hoop  stress 
for  constant  thermal  conductivity  and  spatially  varying  thermal 
conductivity  (C  =  0.9,  n  =  3)  as  a  function  of  3iot  number.  The 
reduction  in  stress  for  varying  thermal  conductivity  is  significant, 
especially  at  the  higher  values  of  the  Biot  number,  at  which  mate¬ 
rials  are  most  susceotible  to  thermal  stress  failure. 


NON-DIMENSIONAL  OISTANCE  ,  r/o 


Fig.  5b.  Spatial  distribution  of  the  transient  thermal  stress  at 

the  tine  of  maximum  tensile  stress,  in  solid  circular  cy¬ 
linder  subjected  to  sudden  convective  heating  for  various 
spatial  distributions  of  the  thermal  conductivity  (equa¬ 
tion  S) . 


Spatially  Varying  Porosity.  The  distribution  of  the  pore 
phase  was  chosen  to  be: 

P  =  PQ(r/a)n  (9) 

The  distribution  for  the  thermal  conductivity  and  Young's  modulus 
was  assumed  to  be  identical  to  the  one  chosen  for  the  steady-state 
heat  flow  given  by  eq.  6.  Again,  Poisson's  ratio  was  assumed  to 
be  independent  of  the  pore  phase.  The  thermal  stresses  were  cal¬ 
culated  by  the  finite  element  program  referred  to  earlier. 

Figures  8a  and  8b  show  the  spatial  variation  in  the  transient 
temperatures  and  stresses  for  80  =  1,  5  and  20,  n  =  1  and  P0  =  0.3 
at  the  instant  of  time  the  tensile  hoop  stresses  at  the  center  of 
the  cylinder  reach  their  maximum  value.  A  value  of  P0  =  0.3  was 
chosen  as  it  corresponds  to  a  value  realizable  in  practice.  For 
purpose  of  comparison,  the  temperatures  and  stresses  for  the  fully 
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Fig.  6.  i'axinura  tensile  therrial  stress  as  a  function  of  coefficient, 
n  in  a  solid  circular  cylinder  subjected  to  convective  heat¬ 
ing  for  various  snatial  distributions  of  the  thermal  con¬ 
ductivity  (equation  -3) . 


BIOT  NUMBER,  fi9 


Fig.  7.  I'a::inum  tensile  thermal  stress  as  a  function  of  Eiot  number 
in  a  solid  circular  cylinder  with  constant  and  spatially 
varying  thermal  conductivity  (n=3,  00.9)  subjected  to 
sudden  convective  heating. 


RADIAL  DISTANCE,  r/o 


Fig.  3a.  Spatial  distribution  of  transient  tenoerature  at  the  tine 
of  maximum  tensile  stress  in  a  solid  circular  cylinder 
subjected  to  sudden  convective  heating  for  spatially  vary¬ 
ing  porosity  described  by  equation  9  with  n=l  and  ?o=0.3. 

dense  material  (?o=0)  are  included  in  Figs.  8a  and  8b.  The  numerical 
results  indicate  that  the  spatially  varying  pore  chase  results  in  a 
decrease  in  both  the  tensile  and  compressive  thermal  stresses  in  the 
cylinder,  accompanied  by  a  simultaneous  increase  in  the  surface  tem¬ 
perature  , 

Figure  9  gives  the  maximum  tensile  hoop  stress  as  a  function  of 
n  for  80  *  1,  5  and  20.  The  stress  levels  for  the  fully  dense  mate¬ 
rial  are  included  for  comparison.  A  major  reduction  in  the  magnitude 
of  the  stresses  as  the  result  of  the  spatially  varying  pore  phase  is 
evident.  Again,  this  reduction  is  greatest  at  the  higher  values  of 
the  Biot  number,  i>>r  which  thermal  stress  fracture  is  most  likely  to 
occur.  The  trend  of  the  curves  in  Fig.  9  suggests  that  it  may  be 
advantageous  to  have  a  porosity  distribution  corresponding  to  values 
of  n  as  low  as  possible.  It  should  be  noted,  however,  that  porosity 
also  has  a  strong  adverse  effect  on  tensile  strength.  This  implies 
that  for  low  values  of  n,  a  considerable  decrease  in  strength  occurs 
very  close  to  the  center  of  the  rod.  In  this  case  it  is  conceivable 
that  failure  can  occur  at  position  in  the  rod  away  from  the  center, 
so  that  the  advantage  of  the  spatial  variation  of  the  porosity  can 
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Fig.  8b.  Soatial  distribution  of  transient  thermal  stress  at  the 

time  of  maximum  tensile  stress  in  a  solid  circular  cylin¬ 
der  subjected  to  sudden  convective  heating  for  spatially 
varying  porosity  described  by  equation  9  with  n=l  and 

P  =  0.3. 
o 

be  lost.  Describing  strength  by  ea.  6  with  b  =  6,  values  of  n 
should  be  less  than  unity,  if  failure  at  positions  other  than 
r  =  0  is  to  be  avoided. 

Temperature  Dependent  Thermal  Conductivity  and  Specific  Heat. 

The  discussion  so  far  has  concentrated  on  a  spatial  variation  in 
thermal  conductivity,  either  introduced  directly  or  indirectly  by 
incorporating  a  spatially  varying  pore  phase. 

A  spatial  variation  of  thermal  conductivity  can  be  caused  by 
another  mechanism,  namely  the  temperature  dependence  of  the  thermal 
conductivity.  It  is  well  known,  that  for  many  dielectric  materials 
in  which  the  conduction  of  heat  occurs  primarily  by  phonon  transport, 
the  thermal  conductivity  exhibits  a  strongly  negative  temperature 
dependence.  For  this  reason,  during  transient  heating,  the  sections 
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Fig.  9.  T'aximum  tensile  thermal  stress  as  a  function  of  n  in  a 
solid  circular  cylinder  subjected  to  sudden  convective 
heating  for  spatially  varying  porosity  described  by  equa¬ 
tion  9  with  PQ  =  0.3  for  various  values  of  Biot  number. 

in  a  structure  or  component  at  the  higher  temperature  will  have  a 
lower  thermal  conductivity  than  those  regions  still  at  a  lower  tem¬ 
perature.  For  a  solid  circular  cylinder  this  should  result  in  a 
type  of  a  transient  soatially  varying  thermal  conductivity,  with 
a  distribution  such  that  in  analogy  to  the  previous  results,  a 
major  reduction  in  the  magnitude  of  the  tensile  thermal  stresses 
should  be  achieved. 

The  validity  of  this  hypothesis  was  verified  by  calculating 
the  maximum  tensile  thermal  stresses  in  a  solid  circular  cylinder 
of  radius  a  =  1  cm  of  high  density  aluminum  oxide  at  an  initial 
uniform  temperature  of  300°K,  subjected  to  convective  thermal  shock 
by  an  instantaneous  change  in  ambient  temperature  to  1000°C.  The 
transfer  of  heat  was  considered  to  be  Newtonian.  The  temperature 
dependence  of  the  thermal  conductivity  of  dense  aluminum  oxide  to  a 
good  approximation  is  described  by  [11] 

K  =  28/T  -  0.01  cal.°K_1 .cm.-1sec-1. 

where  T  is  the  absolute  temperature  in  °K.  The  temperature  depen¬ 
dence  of  the  specific  heat  was  also  taken  into  account  [12].  All 
other  material  properties  were  assumed  to  be  independent  of  tem¬ 
perature. 

The  thermal  stresses  were  calculated  with  the  finite  element 
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program  referred  to  earlier,  with  the  required  modifications  to 
include  the  necessary  iterative  procedure  to  take  the  temperature 
dependence  of  the  thermal  conductivity  and  specific  heat  into  ac¬ 
count  . 

Figures  10a  and  10b  compare  the  time  dependence  of  the  thermal 
stresses  for  temperature  dependent  and  independent  thermal  conduc¬ 
tivity  and  specific  heat,  resp. ,  for  a  range  of  values  of  the  Biot 
number.  The  temperature  dependence  of  the  thermal  conductivity  and 
specific  heat  causes  a  significant  decrease  in  the  magnitude  of  the 
maximum  tensile  thermal  stresses. 

Figure  11  shows  the  peak  value  of  the  transient  tensile  thermal 
stress  as  a  function  of  the  Biot  number.  The  relative  reduction  in 
the  magnitude  of  these  stresses  is  significant,  especially  at  the 
higher  values  of  the  Biot  number. 


Fig.  a.  '’arimnm  tensile  thermal  stress  as  a  function  of  time 
at  the  center  of  solid  circular  cylinder  subjected  to 
sudden  convective  heating  for  a  range  of  Biot  numbers 
with  temperature  dependent  thermal  conductivity  and 
specific  heat. 


C.  Discussion 


The  resuLts  of  the  present  study  suggest  that  at  least  for 
the  heat  flow  conditions  considered,  a  soatially  varying  thermal 
conductivity  can  lead  to  significant  reductions  in  the  tensile 
thermal  stresses.  For  all  cases  considered,  this  was  achieved  by 
decreasing  the  thermal  conductivity  in  the  sections  of  the  geometries 
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Tig.  10b .  '.'axinui;  tensile  thernal  stress  as  a  function  of  time  at 
the  center  of  solid  circular  cylinder  subjected  to  sud¬ 
den  convective  heating  for  a  range  of  3iot  numbers  with 
constant  thermal  conductivity  and  specific  heat. 


''•i't.  11.  tcoendence  of  the  maximum  tensile  thermal  stress  on  3iot 
number  in  a  solid  circular  cylinder  with  temperature 
dependent  X,  Cp  and  constant  IC,  Cp  subjected  to  sudden 
convective  heating. 


studied  which  are  at  the  higher  temperature.  At  first  sight,  this 
latter  fact  contradicts  the  results  of  thermal  stress  theory,  based 
on  spatially  uniform  properties,  that  thermal  stresses  can  be  decreas¬ 
ed  by  increasing  the  thermal  conductivity.  The  results  of  the  present 
study  substantiates  the  significant  effect  of  the  nature  of  the  tem¬ 
perature  distribution  on  the  magnitude  of  thermal  stress.  For  the 
steady-state  heat  flow  reported  in  the  earlier  part  of  this  paper, 
based  on  a  constant  AT  across  the  cylinder  wall,  the  magnitude  of 
the  stresses  is  governed  only  by  the  temperature  distribution  which 
is  controlled  by  the  relative  distribution  of  the  thermal  conducti¬ 
vity  and  not  by  its  absolute  magnitude. 

For  the  transient  heat  transfer,  a  decrease  in  thermal  stress 
is  achieved  by  the  spatial  variation  in  thermal  conductivity,  in¬ 
spite  of  the  fact  that  decreasing  the  thermal  conductivity  in  the 
surface,  in  effect,  raises  the  Biot  number,  which  should  result  in 
an  increase  in  thermal  stress.  Again,  the  strong  influence  of  the 
temperature  distribution  on  the  magnitude  of  thermal  stress  is  in¬ 
dicated.  In  this  respect,  it  should  be  noted  that  comparing  the 
numerical  data  for  the  thermal  stresses  in  terms  of  the  3iot  number, 
8o»  based  on  the  original  unmodified  value  of  thermal  conductivity, 

K0,  in  fact,  is  conservative. 

These  sv  neral  results  of  this  study  possibly  can  be  of  practi¬ 
cal  interest  Improvements  in  thermal  stress  resistance  by  a  SDati- 
ally  varying  pore  phase  already  is  common  practice  in  the  ceramic 
and  refractory  industry.  Low  thermal  conductivity  coatings  can 
also  be  used,  in  which  spatial  variations  in  other  properties  should 
also  be  taken  into  account.  Spatial  variation  in  alloy  content  can 
be  used  very  effectively  to  introduce  a  spatially  varying  thermal 
conductivity.  Spatial  variations  in  the  degree  of  crystallinity  for 
glass-ceramics  can  be  used  effectively  for  the  same  purpose,  as  well 
as  other  techniques  appropriate  for  specific  materials  required  for 
specific  purposes. 
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